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ERRATA. 



DIPPBBENTIAL BQUATIONS. 
P. 37, line 3 from bottom, for d,4f(p) read d^(p), 

46, 6 from bottom, for y=/(jF, o) read if = F(4fy o), 
66, 5 from bottom, for yd,y read yd^y. 
66, 6 /or ( 1 +p«)* read (1 +|)«)*. 

73, 2, 9, 16 from bottom, /or » 1 read = (- 1 )*-> ; and line 1, 2, /or (i»^^ 

read rf,=-.j. 

74, 1, /or dL=« readd^„=^\ and line 10, /or d,g=uM read -d;».«3f ; alio 

line 12, /or y r««rf — y. 
76, 2, /or rf„=. Af read d,^^M; and line 4, for y read — y. 

FINITE DIFFBBENCB& 

P. 23, line 3, for being 1 read being /(I). 
61, 16, for «• rtfa</ a,. 
64, 1, /or these values read this value. 
60," 3 from bottom, for = JT, read ^a'X. 
19, 9, instead of the paragraph beginning A(tt,v«) = Au» . v, + &c. 
read A(tf.o«) = Aw. . «^i + A». . «. = ( A + Atr» ) w^Vh-h 
supposing A to affect u, only, and using Av^ for an operation affecting v, only, 
such that (Atr-i)'v,= A'©^,; 

then A«(w.«,) = (A + Av-»)A(w.«H.i) = (A + Aw-»)*t*.i;^a; 
and, generally, 

A»(tt,«.) = ( A + Aw-» )»tt.«^., 

= A»tt, . w^, + n A— »tt, At;,4 ^i + ** ^ ^ ^•'"»' A'»*+— 2 + &c . 
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INTEGRATION 



DIFFERENTIAL EQUATIONS 

BETWEEN TWO OR MORE VARIABLES. 




SECTION I. 



DIFFERENTIAL EQUATIONS OF THE FIRST ORDER AND 
DEGREE. 

1. In that part of the Integral Calculus which relates 
to the integration of explicit functions of one variable, we 
have to determine the relation between y and of from the 
equation 

in the present portion, we have to determine it from the 
equation 

or to assign the relation between Xj y, z (where ^ is a function 
of the independent variables a^ and y}, or between a greater 
number of variables and their functions^ from the equation 

f(d,Xy dy«y z^ a?, y) = 0, 

or from other equations in which a greater number of variables 
and differential coefficients of higher orders are involved. 

2. A differential equation is said to be of the n^^ order,, 
when the differential coefficient of the highest order which it 
involves is the n***. 

1 
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A differential equation of any order is said, moreover, to 
be of the first, second, &c., degree, when the differential co- 
efficient which marks its order, is raised to the first, second, &c., 
power. 

To integrate a difi^erential equation of any order, is to 
pass to the primitive equation between the variables and the 
constants, from which the proposed may have been derived 
by the process of difierentiation. 

3. We shall begin with the simplest case, viz. that of 
differential equations of the first order and degree, which will 
be of the form 

M and N being functions of x and y. 

Every differential equation of the first order and degree is 
either the direct derived equation of a primitive ; or it results 
from the combination of the derived equation with its pri- 
mitive, so as to eliminate a constant which enters in each only 
to the first power ; the former sort are called exact, the latter 
inexact. 

4. First, let t^ =/(<r, y) = be an equation between i» 
and y, by virtue of which y is a function of w ; then, as is 
proved in the Differential Calculus, d^y is given by the equa- 
tion 

or, since d(*)W, d^^^ii^ are functions of a and y which we may 
represent by M and Ny 

M + Nd^y = 0, 

a dififerential equation of the first order and degree, of which 
/(^» y) = is the primitive or integral. Now if /(a?, y), 
besides other constants which are afi*ected with oo and y, 
contain a term + C independent of w and y ; this will not 
enter into M and A^, having disappeared in difi*erentiating ; 
and if there be no such term, we may add it, and/(.r, y)+C=0 
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is still a relation between ns and y which satisfies the equa- 
tion 

M+Nd^y^O; 

under this form it is called the complete integral; and the 
constant C, which does not appear in the differential equation, 
is called the arbitrary constant ; if the integral did not contain 
such a term as + C, it would not be suflSciently general, and 
would be only a particular case of the complete integral. 

We shall presently give the test which every equation of 
this sort must satisfy, and the mode of integrating it. It is 
evident that no equation of the first order which is not of the 
first degree can be exact. 

6. Next, let Ci be another constant which enters to the 
first power in the equation 

/(a?, y) + C = 0, 

then Ci will be affected with a^ and y, and will consequently 
appear to the first power in 

M + Nd^y = ; 

and if a value of Ci be obtained from either of these equations 
and substituted in the other, the result will be 

an equation of the first order and degree, involving all the 
constants which enter into /(a?, y) + C = 0, except Ci. Hence 
whilst the direct derived equation of 

/(^j y) + C = 0, viz. if + JNrd,y = 0, 

does not involve the term C which is independent of w and y, 
there will be as many other differential equations of the first 
order and degree that have /(a?, y) + C = for their primi- 
tive, as it has independent constants entering only in the 
first power; if any constant enter in a dimension above the 
first, the difierential equation obtained by eliminating it, will 
evidently not be of the first degree. 
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There are two principal methods of integrating equations 
of this sort, which consist either in separating the variables, 
by substitution, or some algebraical process; or in restoring 
the factor which makes them exact. 



Exact Differential Equations of the First Order. 
6. Let 9 be a function of ce. determined by the equation 

then the equation which gives the value of d^y is 

. rf(,)t* + rf(y)t*.d,y = 0, or Jlf+iVd,y = 0; 

the notation d^^^u^ ^{tf)'^^ implying that these differential co- 
efficients are formed on the hypothesis that y is not a function 
of 07, i. e. that a? and y are independent ; then, as proved in 
the differential Calculus, 

Conversely, an equation of the form M + Nd^y s o being 
proposed in which M and N are functions of w and y^ if the 
condition 

(which is called the criterion of integrability) be satisfied, the 
equation results from the immediate differentiation of an 
equation of the form /(a?, y) = 0; if this condition be not 
satisfied, there exists no equation by the simple differentiation 
of which, the given equation can be produced. 

7* To integrate the exact differential equation 

M -f Nd;,y = 0. 
Let the equation from which it is derived be 
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then d^g^u==^M, d^^^u^Nj and d^^^M^d^^N-y 

denoting by F a function of y which may have disappeared, 
since M is the differential coefficient of u relative to Wj on the 
hypothesis that cb and y are independent ; 

.-. dyY^N-d^,^{f^,^M\, e^ndY^f,{N-d,^if^.,M)], 

.% u^{,^M^f,{N-d^,,M)]-^C^O, 
the complete integral involving one arMtrary constant. 

a Obs. The equation Y^fy{N - d^y(f^^^M)\ will be 
absurd, unless the expression N - d^)(J{^M) be independent 
of zr; therefore its differential coefficient with respect to w 
must vanish; 

••• d^M)N- d^^d^y)(S^M) = d(,) JV- d(y)d(,)C(i)-Sf ) « d^^N -d^M 

must equal zero, which it does, since the criterion of inte- 
grability is supposed to be satisfied. 

9. As the simplest case of exact equations, we may first 
notice those in which the variables are separated; they will 
be of the form 

where X denotes a function of w only, and Y a function of 
y only; here the criterion of integrability is manifestly sa- 
tisfied, for 

d^y)jr=rf(„F=0; 

and the complete integral is 

/,jr + /,rd,y;=C, or !,X^f,Y^C. 

To this case may likewise be reduced the equation 
XY,+ YX,d,y^O, 
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which becomes, when divided by ^|F,, 

X r _, 

.-. sm * - + sin * - s sm — , 
a a a 

^ or dn-(fVTl4 + -Vl^)-«'^-'-. 

or ofy/a^ -ry* ^yy/ af - ^s aC; 

at which we may also arrive, by multiplying the proposed 
equation by ooy^ and integrating by parts, which gives^ 

or, since the part affected by the sign j!,, viz. 

is equal to zero by the proposed, 

y y/a^ - «r* + ooy/cf - y* = C. 

Ex. 2. 1 + y + y« + (1 + a? + af^ d,y = 
C (^ + y + 1) « ^ooy + iT + y - 1. 

10. The following are instances of the integration of 
exact differential equations by the method of Art. 7. 

Ex. 1. aa? + fey + c + (6a? + my + «) d,y = 0. 
rf^,^ t* = aa? 4- 6y + c ; 
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.-. u = — -f (fey + c) a? + F, 

d^yj tt = fej? + dy F== fed? + my + w ; 
.'. dyF=my + w; 

.-. F= J^wy^ + wy + C; 
.-. \ {aa^ + my*) + (fey + c) a? + ny + C = 0. 

Ex. 2. .-^ + f- 7^—1 d*9 = 0> 

V/r* + y* Vy ,yva?* + yV 



log (af + \/^ + y®) + C = 0. 



/'-Kggy-f (1 ''wy)d,y 
y^ + Sayit? - a + a* a?* 



Ex.4. ?L±£^±(1j:M^ . 0, 



Homogeneous Equations. 

11. We come next to the case of inexact equations, iii 
which the variables are separable by substitution; of these 
the most important class is homogeneous equations. 

Let M + Nd^y = be a homogeneous equation, that is, 
one in which each of the functions jlf and JV is or can be 
expressed by series of the form 

M = afaf'* + fey'a?'"-" + cy^aT-^f + &c., 
N^ ayfw"f + /Sy^o?''-'' + yy''a'"' + &c., 

the sum of the dimensions of of and y in each term of M 
and JV being equal to r. 
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Let y ^ XX where x denotes a new function of x, then 

^»y = » + xdgZ. 

Hence, making these substitutions in the given equation, and 
dividing by j?', 

/(«) + 0(«) {« + xd^x\ « 0, 

or - + — — « 0, 

^ / W _, „ 

in which the variables are separated. Similarly, the variables 
may be separated by making w^ysfs-*, and the latter substi- 
tution will be more convenient when iV is a more complicated 
expression than M. 

Hence it is easy to e£Pect the separation of the variables 
in equations which are either homogeneous, or can be made 
homogeneous; besides these, the number of equations in 
which that separation is possible, is very limited. 



K. !. 


si^x + a** + 


y»d,y - 0. 


Here 


M 

N 


Sy'a + 


2** X . /x\ 
=.S- + 2 -) 

y \yl 






1 
.-. - + 

X 


'''^ 




S 2 °' 






1 
or - + 


.^''^ ^0, 



^ 3 9. 
-- + -; 



X »* + 3%^ + 2 
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or 



x W + i sr+if 



.-. l<^a? + log («* + 2) - ^log (is* + 1) = log C; 



vVTT 






^*' 



or y* + 2** = C \/af* + »*• 



y« + (a;y + «*)d,y = 0, y = CV» + 2?' 






1 V 

3. (a?-y)d,y = a? + y, tan"^ ^ = log 

5. xdgy - y = V^a?* - y*, sin"'- =log —. 

6. a?dya? + y = v^ + j^, 2Cy + C* = a?*. 

7. \/y + (\/y - V^) d*y = o, 

y — V 2 , /2 \/flp - y/y\ ^ 
^og{y - Vo^y + j») + ^ tan-» ( -^= — j = C 

12. In the following instances, the equations are not ho- 
mogeneous, but are made so by easy substitutions. 

Ex. 1. aj? + fey + c + {mw + wy + p) d,y = 0. 
Let ao? + fey + c = «, mw + wy + p = tJj 
9s and t> denoting functions of x\ then 

a + fed,y = dgXy m + nd,y = d,v ; 
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10 
.-. m + nd^y = d^v . d,» = (a + 6rf,y) d^v ; 

or mv -nx + {bx -^ av) d^v = ; 

this being homogeneous, assume v s xw, w being a function 
of «, 

1 (atr — b)dgW 

in which the variables are separated. 

Ex. 2. a^ (ay + hwd^y) = y" (ay + (i^d^y), 
or (6^-/3y) — «(ay»-aa;^)-. 

Let af^ ^Xi y" = V, ar and v being functions^ of a?, 

m d,« nd^y d,v na? , « 

.-. - = , = — , .-. d,y = - d^v ; 

iv X y V my v 

,\ (bx - Qv) — dgV = av - ax, 
nv 

which is homogeneous. 

Ex. 3. d,y + oy*a?P + fey'a?^ = will become homogene- 

£+1 

neous by making y = «i-», the equation of condition between 
«i, «, p, g, being (p + 1) (l - g) = (w + l) (l - n). 

Ex.4. d,y^ ^ ^ . 
or + afy 



This may be written 



— + . 

y 

and therefore becomes homogeneous when x is written for -. 

y 
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Linear Equations of the First Order. 

13. The next important class of inexact equations of 
the first order which admit of being integrated, are linear 
equations, the general form of which is 

d.y + Py^Q, 

P and Q being functions of w ; they are called linear because 
they involve no pow^r of y above the first. 

Assume y = v«f, v and « being functions of w, 

.\ vdgZ -f xd,v + Pv% = Q. 

Now % being an indeterminate quantity, may be assumed so 
that the equation last written down may resolve itself into 
two others, each of which admits of the separation of its 
variables; to this end let 

vdg,« + P%v = 0, 

d « 
or, dividing by t>, d„z + PiJf = 0, or -^ + P = ; 

So 

.*. logi8f=-J^P, or x = e~f'^. 
The remaining part of the equation gives 

%dgV = Q, or, substituting for », d^v = Qe^"^; 

and y^e"f'''{f,Qef'''+C], 
the complete primitive involving one arbitrary constant. 

Obs. It is unnecessary to add a constant after perform- 
ing the integration indicated in the equation % = e'f"^ ; for 
let X = e-f"^-^"" = C^e'f'^ ; then 

which is the same result as before, since CC^ is equivalent 
only to a single constant. 
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12 
14. If we differentiate the result 

we get ef'^ (d^y + Py) = ef'^Qj 

which shews that if we multiply the proposed equation by 
ef'^9 each member is separately integrable; and this is the 
most convenient practical mode of integrating it. When it 
is once known that the factor which makes the equation in- 
tegrable is a function of w only, its value may be immediately 
found; for let it be denoted by JC; then 

Xd,y + (Py - Q) JT = is exact, 
,-. d,X^d,,,(Py-Q)X^PX, or ~ ^ P. 

.'. iogjr = /,p, or jr=^-^ 

16. It must be observed that if the second member of 
the equation d^y + Py = Q be multiplied by any power of y, 
it is still reducible to the standard form of a linear equation 
of the first order. For suppose the equation to be 

d,y + Py^ Qy% 

then dividing both sides by y", and multiplying by - (w - 1), 
we get 

Hence the factor which makes both sides integrable is 
g-{»-i)/.p^ and the result is 

_L_ = ^»-"/.''{- (« - l)/,Qe-<"->V.'"+ C}. 
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.'. ff = aat + C vl + a»*. 
Ex. 2. <i,y + y = afy*, 



7"' 



■i = a7 + i + Ce^. 



5. d,y = ^y' - a?y, - = a?^ + i + Ce*". 

6 sin <r + cos w r 

6. d,y = a sm c^ + fcy, y = - « . 7—70 + Ce • . 

1+0 



Riccati^s Equation. 
16. There are certain cases of the equation 

(called Riccati's Equation, after the Mathematician who first 
considered it) in which the variables are separable. 
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First, let «i = 0, then d^y ^a-h^, or — ^^f— = 1, where 
the variables are separated. 

Secondly, let wi = - 2, and assume y = — , u being a 
function of w\ 



•*. ^dgU — —- + — •*- — , 

w a? cf 3^ 

where the variables are separated. 

1 u 
Thirdly, let m = - 4, and assume «=-—+--, 

6a? ^^ 

, 1 2u dgU 1 2t^ bu^ a 

where the variables are separated. 

17* Besides the above, the variables are likewise separable 

-4i 

in the cases when m = — : , i being any integer from to 

2 e sfe 1 cf ^ o 

infinity; all which values of w evidently lie between and - 4. 

- 4fi 
First, let m = 



2t-r 

Assume y = — + -r— ; 
oa? aru 

1 2 d,.«^ 



12 6 

6y^ - aai^ = — i + -^— + -r-r - fltX''"; 
^ bco^ aru oo^u^ 
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. adding 


these 


together, 










b 


-aa?" - 





or a^dgU + au^w'*'^* -6 = 0. 

1 m + 2 

Now, let a?«««+3, then d,t*«d,t*d,»s=(m + 3)«^+3 d,w; 

a b "»+^ 

.*. d,t* + u^ = ;8f"'*+3, or d,«^ + fei«^* = aj»'">. 

m + 3 m+ 3 



-4i ^ m+4 4i-4 4(i-l) 



Hence by these substitutions the equation is transformed 
into another of exactly the same form, with i - I instead of i 
in the index of the variable in the second member. 

1 1 -J- 

Similarly, by substituting — - + — -r for w, and Xi^^"^^ 

for j^, we shall transform the equation into another of the same 

form where m^ = —r- ;; — ^ ; and consequently, after i sub- 

^ 2(^-2)-l ^ ^ 

stitutions the index of the variable in the second member will 
become zero, and the variables will be separated. 

-4i 
Secondly, let m 



Assume y 

u 



2i + l 
1 



, d^u b 
then — + -^ = aoT, 



Digitized by LjOOQ IC 



16 

1 m 

let a? = «^+ » , then d^u = d,M d^x = (m + 1) «»«+ * d^u ; 



.-. -(w +!)»« + 


1 d,«^ + 6 = o»«+* tt% 


or a,«^ + 

m + 1 


m + 1 


) 






-21 + 1 
... ^ + 1 = -^ — . 

2» + 1 


2i + r 


m 


-4i -4» 



m + 1 -2i+l 2*-l 

Hence by these substitutions this case is reduced to the 

— 4i , . 

former; and therefore when m = -: , the variables in the 

2i±l 

equation d^y-^hy^^aoT can be separated. It may be observed 
that the more general equation, d,y + b't^ofi"^ = aai^^ is reducible 
to this form by putting a?* = %, 

18. We shall now give some other instances of equations 
in which the variables are separable by particular substitu* 
tions. 



Ex. 1. a (pffd,y - y) = (ar + yd^y) y/w^ + y« - a^ 
When an equation contains the expressions 



yd^y + a?, xd^y - y, \/^ + ^, 

the introduction of polar co-ordinates will sometimes effect the 
separation of the variables ; that is, to assume 

w== p cos Of y - p sin 69 

p being supposed a function of 0, for then 

+ ofdeof = pdep^ xd^y - yd^w = p^. 
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Hence the proposed equation, which considering w and y 
as functions of 0, may be written 

a {wd^y - yd^ai) = {wd^w + yd^y) y/a? + y* - a*, 

becomes ap* = pd^p w f? - a*, 



or o a — ^ \//o^--a* ; 

.-. a0 = \//o?^ -- o* - « sec"* - + C, 
"^ a 

or a tan-> ?^ = y/lFT^^^ - a sec^* "v^ + tt ^ ^ 

/ w 

Ex. 2. y - ird,y + d,y \/ (ax + by) - = 0. 



/ ^ 

or ydyW - a? + \/ (ao? + hy) - « ; 

let - = «r, j8f being a function of y, 

then ydyX - or = y^rfy^f ; 

.-. f^dyX + \/(a«f -f 6) jzry s 0, 
where the variables are separated. 

_, ^ o* + a;* - y* 
Ex. S. d,y :«= ^— ii- . 

Since this is satisfied by y = «, assume y » j? + z'^j 
or rf,« a -- , which is a linear equation ; 
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.-. 




1 


I.e'^^ 


c. 












or 




1 


/.« 




c. 








y — w 




Ex. 


4. 


(1- 


-wy) 


d.y-k-y' 


+ aj? 


-0. 










Assume 


» = 


X - aa? 


9 


so 


that 


» = 


y_ 

1 


+ Oi»* 




1 +tr» 


-«y 








then 


d.« 


+ 




07 
a 


= 0, 







where the variables are separated. 

n(l + y«)l 



Ex. 5. (y-^)dxy = 



assume y 






1 +^» 



(l+««)(i?? + »\/lT^ 1+^ 



0. 



Euler*s Equation. 
19. To integrate the equation 

or, considering o? and y as functions of a new variable ty 

dty \/5^ + d^wy/y^ 0. 

Let the function of i which expresses w be determined 
■" by the equation d|07 = \/5^ and therefore that which ex- 
presses y by the equation d^y = - \/F» ^^o let w + y^^Pf 
07 - y as 9, p and q being functions of /. 
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Then since ((^ar)*B^; .: ZdtW d^x^dtXy at dfmm^d.Xx 
similarly dfy^^d^Y. 

••• d,*p^\{d,X^d,Y) 

= i {26 + 2c (» + y) + Se(a»» + y^ + 4/(«» + y»)^ 

= 6 + c (^ + y) + — {(a? + y)« + (a? - y)H 

^-^fip + y)' — ^^—r^ — — 

4 

Ss 
- 6 + cj» + — (P* + 9*) + i/P (P* + 3?»). 

and dtp . dj^r-JT- F=.fc(»-y)+c(4f«-y») +e(«»-y»)+/(**-y») 

H ^ 

.'. dtp^qx/C + ep-^-fji^, 

or ^/X-'\/Y^(af^y)^/c + e{a; + y)+f{w + yy, 

the integral required. The discovery of this integral, which 
IS due to Euler, was of great importance, as being the first step 
towards the foundation of the Theory of Elliptic Functions. 

20. To integrate the equation 

\/l — c* sin* yjf + y/l - c* sin* (p d^y^ = 0, 
or considering.^ and \^ as functions of. another variable t^ 
VI - c* sin* >^ dt(lf + \/l -c*sin*0 d| \j> * 0. . ' 



or 
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Let d^<f>^\/\^^s\v?ip, and .-. d^x//* - \/l- c^sin*x|r; 

r. d^(f> + d?y\r = - i^ (*^" ^^ + sin 2 >;^), 
rft*^ - d?y\f ^ -• ^c* (sin 2^ - sin2x|^). 

Let pa^ + >/f, q^<l>-'^9 

.'. d^p = — c^ sin p cos q^ 
d^q = - c?*cosp sin g, 

and dtp . d^q = (dj^)* - (d^x/ry = - c^ (sin'»^ - sin^ x|^) 
=s — (cos20 - cos2>^) « — c* sin j> sin q. 

dt^p cos qdtq , /j x i /^^ • x 

dtP sin q 

.-. dtp = C sin gr ; similarly d^gr = C sin p ; 

.-, -y/i - c* sin* =f \/l - e*sin*\|/^ = C sin (^ ^ >^) 

18 the integral of the proposed equation; which is only 
Euler^s equation under a different form. 

The equations \/l +^ + \/l + ^* d,y = 0, 

are immediately reducible to the above form, vi^. 

\/l - ^ sin* >^ + \/\ "• ^sin*0 d^>/^ = ; 

the former by making ^»tan^0, S^-tan^x^; the latter 
by making y/oi « cos 0, y/y = cos \|r. 
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On the Factors which render integrable a Differential Equation of the 

First Order. 

21 . The most natural way of obtaining the complete in- 
tegral of a differential equation of the first order, is to prepare 
it so that its first member may become an exact differential 
coefficient ; for then we shall have only to integrate and add 
a constant. This preparation is always possible by means of 
a factor, when the equation is reduced to the form d,y + jr=0. 
For let an equation /(a?, y, C) = be resolved with respect 
to C, so that 

C«0(a7, y); 

P 

.*. by differentiation, p = P + Qd^y, or d^y + -- = o. 

Now the equation M + Nd^y = may be put under the 
form d,y + iT = 0, which agrees with the preceding, and may 
consequently be supposed to have arisen from the elimination 
of a constant between the primitive /(a?, y, C) = 0, and its 
immediately derived equation. On this supposition, therefore, 

P 

d^y + JT = is identical with d^y + — = ; 

.-. d^y + K = Q~^' 

or d,0 (.Vj y) = Q (d^y + JT), identically. 

The second member therefore is an exact differential 
coefficient, which proves that there always exists a factor 
proper to render the expression d^y + K integrable. 

22. But although the existence of the factor in every 
case is thus established, the investigation of it is usually 
attended with greater difficulties than the solution of the 
original, equation. 

For let P + Qd^y = be an exact differential equation ; 
and let %, a function of a? and y^ be a common factor of 
P and Q, so that P = Mxy Q = Nx, by the removal of 
which, the equation is reduced to the inexact state, 
M + Nd^y = ; 
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then because P + Qd^y » is exact, 

or xd^^M •{■ Md^^x^zd^g^N^Nd^g^^x^ 
or iVd!^,^»-ifd^y^» = «(d^y^if-d^,^JNr); 

an equation between <r, jf, j^r, and the partial differential co- 
efficients of Xy for determining the factor x. The consideration 
of this equation in its general state must be reservied till we 
come to treat of partial differential equations of the first order ; 
but the following particular cases may be noticed. 

23. First, suppose that the factor is a function of only 
one of the variables a?, then d^y^x = 0, and the equation be- 
comes 

d X 1 

which, being integrated, gives x ; for the hypothesis requires 
that the second member should be independent of y. 

Similarly, if the factor be a function of y only, it will 
result from the integration of 

of which the second member is independent of w. 
Hence, if in any equation M + Nd^y = we find 

a function of w only, or 

a function of y only; the factors which may make it integrable 
are respectively c-^'^, e^*^. 



Digitized by LjOOQ IC 



23 

Ex. 1. d^y + {Py -0)^0^ the linear equation of the 
first order. 

This compared with M + Nd^y « 0, gives 
M^Py-Q, J\r=l; 

.-. d^^^M - d^,^N = P^ and 1 (d^y^ JIf - d^^^N) « P, 

a function of x only; therefore the factor is e^'^> 

Ex. 2. »* + (!- ^y) rf,y = 0, 

M=iy\ JV «= 1 - ^y, 

therefore the factor is ,. 

24. In the case of homogeneous equations, a factor 
proper to render them integrable, is readily discovered by 
^ means of the property that if t^ be a homogeneous function of 
n dimensions of the independent quantities t and x, then 

nu ^ tdtu + xdgU. 

Por suppose F, a homogeneous function of w and y of m 
dimensions, to .be a factor which makes M + Nd^y an exact 
differential coefficient, M and N being homogeneous functions 
of w and y of r dimensions ; then if U denote the primitive, 
it will be homogeneous and of 9i» + r + 1 dimensions, and we 
fihall have 

VM^ VNd.y « d,Uy d^,^U^ VM, d^^U ^ YN, 

^Fif + y FJV » (w + r + 1) (7; 

M+.Nd,y^ 1 d^U 

Mw + Ny w + r + 1 ' (7 * 
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md as the second member is an exact differential coefficient, 
it follows that the first is so likewise, and consequently, that 
M + Nd^y is made exact by means of the multiplier 

1 
Mof + Ny ' 

26. The property of homogeneous functions assumed 
above is easily proved. Let uhe a. homogeneous function of 
the independent quantities t and z of n dimensions ; then if we 
change t into t{l + h) and % into i^r (l + A), u will become 

w (1 + A)" = w + nuh + &c. 

But by Taylor's theorem, u. will also become 

u + d^u . ht + dgU . hx + &c. 

therefore, equating the coefficients of A, 

nu = tdiU + xdgU. 

And, generally, if t^ be a homogeneous function of 9i di. 
mensions of any number of independent quantities t, », w?, &c., 
and we change them into ^ (1 + A), x (l + A), w {l + A), &c., 
the new value of u will be equally expressed by w (l + A)* or 
by e(** + '^ + '")^w; and equating the coefficients of A' in 
these two identical expressions, we get, separating as above 
the symbols of operation from those of quantity, 

w(n-l)...(7i-r + 1) w = {tdt + xdg+ tvd^ + ...)r^- 
Ex. /py + y* + {ofy - a^) d,y = o. 

The factor is 

1 1 



(wy + y^)w + (afy -ai^)y 2y^af' 

2y*a? 2y*^ 

is an exact differential coefficient, and gives the primitive by 
Art. 7. 

^ + |log(j?y) + C«0. 
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26. Whenever the variables can be separated in an equa- 
tion, a factor which makes it integrable can also be found. 

For suppose that M + Nd^y = 0, by the introduction of 
two other variables u and «, is transformed into R + Sd^u » 0, 
so that 

M + Nd^y - JJ + Sd^u ; 

and suppose F to be a function of u and x, such that if 
we divide R + Sd^u by it, the variables are separated, i. e. 

R S 

— contains % only, and - contains u only; 

is an exact differential coefficient; and consequently — , which, 

upon restoring the values of u and %y becomes a function of 
w and yj is a factor which makes M + Nd^y » integrable. 

Ex. 1. a + b^f^ + o^d^y = 0. 
Assuming ^ » — > we find 

OB 

a + ba^i^ + cfd^y = a + few* + wd^u - «, 
and dividing by j? (a - 1* + btf) we get 

a + ba^i^ + (x^d^y 1 d^u 



d? (a - w + 6w*) ao? - a^y + 6^y* 
is a factor which makes the proposed equation integrable. 

Ex. 2. y* + oo? + (1 - xy) dj,y « 0. 
4 
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By assuming y = , it may be shewn that a factor 

which makes the proposed integrable is 

1 

f^ + Say Of - a + a^a^ 

Ex. 3, M + Nd^y = 0, a homogeneous equation. 
In this case we know, that making y^wssj we have 

and M -k- N d,y = arf{%) + ofip (x) {% + aid»%\ ; 

consequently, dividing by af'^'^ {/ W + »0 (*)} = ^^ + ^V^ 

we get 

JIf + Nd^y _ 1 (g) d^% 

Mce + JVy a? f{%) + « ^ («) ' 

1 

JIf 07 + Ny 

is a factor which makes the proposed equation integrable. 

Obs. That -rrz r^ is an exact differential coefiSdent^ 

Mx'\-Ny 

provided M and N be homogeneous functions of m and y of 
the same dimensions, admits of an easy proof as follows. 

We must shew that d(«(j^^f^) - \-'\jt^^' 
Now putting - = ;», we have 

M 1111 



Mw + Ny w Ny ^ * 1 + F(») 

iV^ 1 / _ J/g? \ _\ _\ 1 

+ iNry"y\ " ilf4? + Ny) " y" yX-^r F{%) ' 
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d.F(«).- 



( M \ 1 ^ ' a> 

''' ^'^\Mai + Ny] "~x' {l + F (»)}*' 

which expressions are evidently equal to one another. 

27. There always exists an infinite number of factors \/ 
which render an equation of the first order and degree inte- 
grable. 

FcMT let « be the factor by means of which the equation 

is made integrable, and u^O its complete primitive, so that 
« ( Jf + Nd^y) - dgU ; 

multiply both sides by -F(w), where F{u) denotes any func- 
tion of Uj and we have 

« F{u) (M + Nd,y) * F {u) d,u ; 

and as the second member is an exact differential coefficient, 
it follows that the first is so likewise; therefore zF{u) is 
a factor which makes the proposed equation integrable, what- 
ever form be assigned to F{uy 

28. The following geometrical problems are added to 
illustrate this part of the subject. 

I. To determine the trajectory of a given family of 
curves. 

Let AA\ BB!^ (fig. 1.) be two of a family of curves 
resulting from the equation /(JT, F, c) =0, by given par- 
ticular values to the constant c ; and let AB be a curve which 
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cuts AA\ BB^y and all the curves resulting from the equation 
bj giving all possible values to c, at the same angle ; then 
AB is called the trajectory of this family of curves. Let 
Wj y, be the co-ordinates of the point A in AB^ between which 
we are required to find a relation ; AT'f AT^ tangents to the 
curve and trajectory at A^ tan TAT' ^ a; and let a value 
>^(jr, V) of d^Yhe obtained from the equation /(-X", F, c) « 0, 
not involving c; then at the point A, tan-^T^JV* >^(a?, y), 
and tsai ATN 'B d^y, 

d^y-yj/ (a^, y) 
'^ l^d^y.ylf{w, yY 

the differential equation to the curve AB ; and as it does not 
involve c, AB will cut every curve in the series at an angle 
whose tangent >= a. The equation when integrated will in- 
volve an arbitrary constant, and consequently will represent 
a system o^ curves, every one of which cuts the former system 
at the same angle; the constant may be determined, if a 
point through which the trajectory is to pass, be given. If 
the angle TAT' be a right angle, or (a) infinite, the differ- 
ential equation to the trajectory, which is then called ortho- 
gonal, is 

1 ^d^y.y\f{wy y) = 0, 

Ex. 1. To find the orthogonal trajectory to all curves 
resulting from the equation y* (c - j?) = aP^ by giving all 
possible values to c. 

Here c- 47= — , .•.-! = -— j d^y; 

y" y" f 

.-. d,y = xf.(^, y)= ^^ , 

therefore, substituting for y\f (a?, y) its value, the differential 
equation to the trajectory is 

2a?^ + (y^ + Sa;^y)d,y = 0; 

a homogeneous equation whose integral is 
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Similarly, let /(p, 0, c) = be the equation to the curve 
A A referred to polar co-ordinates, and let it give for pd^Q 
the value ^ (p, 9) independent of c ; then considering p and 
as co-ordinates of the point A in the curve AB, 

tan SAT = yff (p, 0), tan SAT '^^ pd^d^ 

which is the differential equation to the trajectory; or if it 
be orthogonal, 

l+y\f{p,e) pd^e = 0. 

Ex. 2. Let the curves be a system of circles touching a 
straight line in the same point, then taking that point as the 
origin and measuring d from the line, their equation is 



1 COS0 J /I ^^^ I / m 

p smd *^ '^ ^ COS0 ^ ^ 

1 + ^pdpO^O, or cosd + &in0pdp6^Of 



or do I ^ 1 = ; .-. p = C cos 0, 
'^ \cos 0J '^ 

the equation to the orthogonal trajectory, which represents 
a system of circles passing through the given point and having 
the given line for their diameter. 

We may generalize this problem, by finding the orthogonal 
trajectory of all circles described through two given points. 

II. To determine a curve such, that the locus of the 
extremity of its polar subtangent shall be a straight line. 

The polar subtangent is a line drawn from the origin 
perpendicular to the radius vector to meet the tangent. 

Let pf 0y be the polar co-ordinates of any point P in 
the curve sought (fig. 2) ; then those of the extremity T of 
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its polar subtangent wfll be p*dp9 and - - , which must 
satisfy the equation to a straight line, via. 

^ '^ \ 2/ sin(0-a) 

c c 

.-. — = sin {Q - a) d^0, and - = cos (0 - a) + C; 

P P 

c 
''' ''"C + cosCe-a)' 

the equation to curves of the second order, having the pole 
for one of their foci. 

III. To find a curve in which SG varies as iSP, PG 
being a normal at P, and SG a fixed line through S. (fig. 2). 

Taking SG for the axis of ^, the equation to the normal 
at P is 

therefore making 

F«0, X^SG^of + yd^y, 

and iS^P ■» \/a7^ + y*, .*. a? + yd,y = c \/ of ^ y^, 



or VV + j^=ca? + C, 
the equation to a curve of the second order. 

IV. To find a curve which is always cut by its radius 
vector at an angle proportional to the corresponding angle of 
revolution; that is, I SPT oc z ASP^ (fig. 2). 

Let p, d, be the co-ordinates of any point in the curve, 
then the angle at which the radius vector cuts the curve, has 
for its tangent pd^O; 

^ /I * ^ ^ep cosnd (pY . ^ 

.'. pci^0 = tann0, or —J^^- -, .-. \^\ =sinw0. 

^ '^ p Siund \CJ 
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SECTION II. 

DIFFERENTIAL EQUATIONS OF THE FIRST ORDER, BUT NOT 
OF THE FIRST DEGREE. 



29. When a differential equation of the first order is of 
a higher degree than the first, we know that it is not obtained 
by the direct differentiation of its primitive, but results from 
eliminating a constant, (which enters into the primitive in a 
dimension above the first,) between the primitive and its de- 
rived equation ; the degree of the differential equation and the 
dimension of the constant eliminated above the lowest dimen- 
sion in which it appears, being always the same. The general 
form of such equations free from radicals, is 

(d^yY + Pi (^,y)""' + P2 (rf*y)"'* + •.. + Pn-idsV + p. = 0, 
the coefficients being functions of x and y. 

If this can be resolved with respect to d^y into its {n) 
simple factors, it will assume the form 

{d^y + ^i) {d^y + 92) — {d^y + ?«) = ; 

then each of these factors put equal to zero, will be an equation 
of the first order and degree, whose integral may be found by 
the methods of the preceding section ; and any one of these 
integrals, as well as the continued product of any number of 
them, will evidently satisfy the proposed equation. If, there- 
fore, we integrate the n equations, 

d»y + gi = 0, d,y + gg * 0, ... d^y + g» « 0, 

and complete them all with the same constant C, as the 
proposed equation is of the first order, we shall obtain the 
required primitive invcJving C in the «*'* power, by equating 
their continued product to zero. 
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Ex.1. (d.j^)* + — d,y-l = 0; 

vy^ + o?^ 

.-. + \/^^r^= cT + C, and -V^^r^=^ + C; 

.-. (\/fT^ - C - a?) (V^T^ + C + 0?) = 0, 
or y2 + a?^-(C + a?)* = 0, 

or j/* = 2Ca? + C«. 

Ex. 2. a?^ (d,y)' - 2a?y d^y^f- a^tf - ^ = 0. 
.-. a?d,y - y == ± ^ \/^TyS 



d. 

or =1- 



:{^^^ 



... iog{^ + Vi+ - h^+c, 



.•. ^/a^ '\-f^ + y = ca?e*, changing the constant, 



and VV+y - y == c^e* ; 
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30. When the resolution of the proppded equation into 
its simple factors is impossible, there are still various forms 
for which the complete primitive can be determined, or its 
determination made to depend on elimination; this is done 
by means of substitution, or differentiation, or other analytical 
artifices, of which we shall now give some instances. 

Obs. The arbitrary constant in what follows is often 
reserved under the sign of integration. 

31. When the equation contains only one of the variables, 
<r suppose, and can be solved with respect to that variable, 
so that w -fid^y) ; let d^y be denoted by p, then 3B-f{p) ; 
and integrating the equation d^y = p by parts, we get 

y^xp^ f^ad^p = pf(p) - fpf(p) ; 

between which and the equation a? =/(j?), eliminating p, we 
shall obtain the required integral. 

Similarly, if we have y«=/(p), since 
dpOff^d^w.d^y^^d^fip), 

we shall have to eliminate p between y^fip)^ 

/I 
- dpf(p). 
pP 

Ex. 1. 0? + j?(d,j)*= 1, 

1 

.'. « - \/m (1 - «) - tan-' V -^ + C. 

w 

Ex.2, y « a -x/i + p\ a? + C - a log (\/y* - a* + y) . 

This is the solution of the problem in which it is required 
to find a curve such that the perpendicular on the tangent 
^^from the foot of the ordinate shall be constant. 
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f / 



Ex.3, yV'l+p*=ap, 



(1+P')i' 

or x + C = \/a*- y* + a log ( . ] . 

This is the solution of'the problem to find a curve such 
that the tangent terminated at the axis of w shall be of a 
constant length. 

32. An equation not coming immediately under this 
case, may sometimes be reduced to it by putting p *= wzj or 
p^yx. 

£x. {d^yy + aa d^y -^01^ = 0. Let p = a?«, 



then a? («* + 1) + a» = 0, or j? = - 
.-. d^y = d^y . rf,d? = - 



ax 



1 + 



(z" + I)' ' 
and X must be eliminated between the integral of this, and the 

ax 
equation a? = . 

33. When the equation contains both the variables <r 
and yy provided it be homogeneous with respect to them, we 

may assume - » 2? ; then the equation will take the form, 

(which is not solvable with respect to p by supposition,) 

f{x,p)^0. 

Suppose this capable of being solved with respect to x^ and 

let it give x = (p (p) ; now y = wx gives p^x -^ wd^x^ 

\ d^ 
or - = 



CO p - X 

substitute the above value of Xj and integrate this equation ; 
then p must be eliminated between the result, which will be 
of the form log ^ = jP (p)^ and y =saf(p (/>). 
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Ex. y - wd^y ^nxy/l^ (d^yy ; 
y /- 

.-. - = p + 7lVl +j[>*; 



.-. log.v = - - {log (;> + \/l+pl) + wlog\/l+p«} + logC; 

V 1 +J!>* ^ 

between which equations p must be eliminated.' 

This is the solution of the problem, to find a curve such \/ 
that the perpendicular upon the tangent from the origin 
shall vary as the abscissa to the point of contact. 

Letw=i. -=l+p2+pv^l+p% l^p + y/T7p; 

OB X 

or '^ ^9.CoG -sf. 

34. Another integrable form is y = ood^y +f{d,y), which 
is called Clairaut's form, after the Mathematician who first 
considered it. Substituting p for d,y, and differentiating, 
we get 

y = a}p+f{p), 
drV-^p^p-^- wd^p + d^f{p) . d^p; 
••• {^'^dpf(p)}d,p^O, 
which resolves itself into the two 

•^ + dpf(p) = 0, d^p = 0. 
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The first of these gives p^^ <f>{ai) supppse ; this value 
substituted for p in the proposed equation, furnishes a relation 
between x and y which satisfies the proposed equation, but 
which involves no arbitrary constant, and cannot therefore 
be the complete primitive. The other equation must there- 
fore lead to the complete primitive ; but this gives p = C, 
and by substituting this value of p in the proposed we find 

y^Cw+fiC). 

Hence Clairaut^s form has the property, that the complete 
primitive is obtained by substituting the arbitrary constant C 
for p, in that form. If we integrate p=»C, we find y=^Cx+Cf; 
but the condition of the proposed equation being satisfied gives 
Cf = fiC)y the same result as before. 

We shall afterwards return to the consideration of the 
other solution, which is called the singular solution, and is 
not derivable from the complete integral. 

Ex. 1. y = d?p + — ^ ^-^. 

P 

Differentiating, we get p = p + wd^p + o ( 1 — A d^py 
or fa? + a ;]rf,p = 0; 

.*. d^p « gives p^ C, and y^Ca^- — ~ , 

the complete integral, and 

a? + a--T=0 gives p = ^^ ?L- 

P^ w + a 

which, substituted in py = o + p* (a? + a), gives 



the singular solution 



=*= y V « 2o, or ^ « 4.« (a? + o), 

^ + a 
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Ex. 2. y ^ xp-^- y/V + a^p*, 

y = Ca?+\/6* + C®aS the complete integral. 

a*y* ^b^w* a a^b* the singular solution. 

This is the solution of the problem to find a curve, such 
that the product of the perpendiculars dropped from two 
given points upon the tangent may be invariable ; for taking 
the line joining the two given points (whose distance 
suppose a 2o) for the axis of /r, and their middle point for 
origin, and /r, y the co-ordinates of any point in the curve, 
the equation to the tangent at that point will be 

Y-y^d^yiX-w), or F = rf,yjr + (y-a?d,y), 

and the lengths of the perpendiculars dropped upon this 
line from the points {0^0)9 {-c^o) will be 

^ pc — (y — wp) pc — (y — wp) 

and the product of these is 

- — ^ o — ^ = 6 J suppose ; 
1 + p 

.*. y- xp ^ y/}? + ©^p*, putting a* = 6* 4- c'. 
35. A still more general case is the equation 

which by differentiation is reduced to a linear equation of 
the ^rst order in a\ for we get 

P "^/(P) + ^dpf(p) ^»P + dp<^ {P) ^xP ; 
••• {P -fiP)] dpOf = ^dpf(p) + dpip(p) ; 

' f(p)'P f(p)'p' 

which gives w^F(p); then p must be eliminated between 
this and the proposed equation. 
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Ex, 1. y =■ wp^ + 2p, 

or rf^a? + 



p — 1 p^ — p^ 

which is made integrable by the factor (p-l)'; 

/•2o — 2 
.-. a.(;>-l)*=-. /--^ = -2;> + Iogp» + C. 

But p =s: ± \/ - + -- ; therefore substituting this 

^ iV or 

in the preceding, we obtain the complete primitive between 
Of and y. 

Ex. 2. y = /pmp 4- n v^l +p^ ; 

2m— r 
»» .^ *v »»— 1 

jiTTT^ n r p^ 

Ex. 5. y'\-p((i-ai)^nj^ \/\ +p^\ 

^ " 2(^ (n + 1) "^ 2 (n - 1) (o - a?)"-» "^ 

This is the solution of the problem of finding the path of 
a point P which moves uniformly towards another point Q, 
also moving uniforiply in a straight line. 

For taking A (fig. 3.) for the origin, and AB^ which is 
perpendicular to By the line in which Q moves, for the axis 
of a?, we have, supposing P and Q to start together from A 
and B, BQ = nAP, or if AN = a?, NP = y, ^J? = a, 

y ^ (a - .v) p == n l,y/l -^ p^. 
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36. In the following examples the method of substi- 
tution succeeds. 

Ex. 1. (1 - p*) a?y = jp (a?" - y* - c*), 

which expresses that the normal bisects the angle between the 
focal distances; 2c being the distance of the foci, the origin at 
the middle point between them, and the line joining them the 
axis of X. 

xz % 

Let jp = — , .*. y* = a^x - c^- 



y 1 4-« 

therefore, differentiating, 



this resolves itself into 
c 



]to 

, which gives y* + (^ - c)* = 0, 



1 4- «' 

the singular solution; 

and dgX = 0, or jjr = C, which gives y* = C f a;* — j , 

the complete integral. 

If we integrate p = — , we get y' = Ca^ 4- C'y where C 

must be determined by the condition of .the proposed equation 
being satisfied ; and by this condition, in general whenever the 
method of solution raises the order of the equation, must the 
number of constants be reduced. 

By the same substitution may be solved the more general 
form 

awyp^ + p(a^ ^ «y* - 6) — a?y = 0. 



Ex. 2. xd,y - y = Jr \/{d,yy - — d,y + i , 
or d,y"??«- vfd,y-?^V+l-?[. 

XX \ x) c2?* 
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y 

Let y^w«9 then d,y = iif + a?d,« = - H- «pd,»; 



If Jr= 1, we have sin-'« = sec"'* + C, 



or sin 



■» ? -= see-'* + C. 



Ex. S. and 4. 



=/(vvT?), 



y - a;p / ^ y 

Introducing polar co-ordinates, we get for the first, 
-^=4^==/(p), or d,e=- //?^, M. 
The second gives 



The former expresses that the perpendicular on the tangent 
from the origin is a given function of the radius vector ; the 
latter, that the sine of the angle at which the radius vector 
cuts the curve is a given function of the cosine of the angle 
at which it is inclined to the axis of w. 
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SECTION III. 

ON THE SIN6I7LAR SOLUTIONS OF DIFFERENTIAL 
EQUATIONS. 



37. From the complete integral of a differential equation 
we can deduce as many particular integrals as we please, by 
giving to the arbitrary constant particular values. But some 
differential equations are satisfied by a relation between x and 
y not containing an arbitrary constant, and not deducible from 
the complete integral ; such a relation is called, as has been 
said, a singular solution of the differential equation. The 
existence of such solutions depends upon the fact, that when 
a solution of a differential equation has been obtained in any 
manner, it will still be a solution after a quantity of any kind 
has been introduced in any way, provided the same derived 
equation result. This is merely an extension of the principle 
on which the arbitrary constant is added. 

38. Before entering upon the general theory, it may be 
useful to consider the following particular instance* 

Let the equation 

y = sed^y + a + a (d,y)* (l) 

be proposed, which, since it falls under Clairaut^s form, has 
for its complete integral 

y^^Cx^a^-aC^^ (2) 

C being the arbitrary constant. If we now regard C, not as 
a constant, but as a function of ^, and differentiate, we get 

d^y^C^'{ai'\'%aC)d^C\ 

and if we eliminate C between this and y ms Caf + a + aC^^ we 
shall obtain a differential equation, but not the proposed one, 
6 
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for that arises by eliminating C by means of the equation 
d*y - C. But if C be so determined as to make the co- 

Of 

efficient of d,C vanish, that is, if C = , then the derived 

equation will be d^y - C, and the result of the elimination 
of C will be the proposed equation. 

Substituting for C its value, we get 

«5= +a + — = + a, 

Sa 4a 4a 

a relation which manifestly satisfies the proposed equation; 

ffar it gives d^y « — — , and these values of y and d^y, being 
2a 

substituted in the proposed equation, make it identical. But 
this solution contains no arbitrary constant, and yet being the 
equation to a parabola it cannot, either by making C = 0^ 
or any other constant quantity, arise from the complete in- 
tegral which is the equation to a straight line ; it is conse- 
quently a singular solution, and arises from the complete 
integral by changing C into a fcinction of /r so determined as 
to make the term involving d^C disappear from the value 
of d^y. 

Thus we see how the singular solution arises from the 
complete integral ; next, let us consider its geometrical signi- 
fication. The proposed differential equation expresses that 
the curves to which it belongs have the property, that the 
tangent at any point is intersected by a perpendicular upon 
it from a given point in a given straight line. 

Take the given point S (fig. 4.) for the origin, and TS^ 
AS respectively parallel and perpendicular to the given line AC, 
for the axes of w and y. Let TC be the tangent at a point 
who^e co-ordinates are w and y, then its equation is 

F-y = d,y(X-j?), 

and the equation to the perpendicufar upon it from S usk 

Fd,y=-^, 
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and for their point of intersection 

but this point is always in AC for which F ■= a ; 

the same as the proposed equation. 

The complete integral jf = C^ + a (l + C'), which repre- 
sents a series of straight lines, evidently satisfies the problem 
for all values of C\ for let Tt be any one of these lines, 

then y as - _- is the equation to a line through S perpen- 

dicular to it; and combining the equations to get the co- 
ordinates of their point of intersection, we have 

y = -Cy + a(l + C*), or y = o; 

the intersection consequently falls in AC- And as a straight 
line is itS own tangent at every point, the equation 

y s Cos + a (l + C*), for all values of C 

represents a line such that the intersection of the tangent 
at any point, and a perpendicular upon it from S^ falls in 
the given line AC. Now the curve which is generated by 
the perpetual intersections of these lines will also satisfy the 
problem; for each of the lines will be a tangent to it, and 
therefore perpendiculars from S upon its tangents will inter- 
sect them in the various points of AC. To get the equation 
to this curve we must, according to the usual method, dif^ 
>/ ferentiate with respect to the parameter C, which gives 
sz cp + 2Cay and eliminate C between this, and the equation 
y = Cof + a (l + C^), which gives 

y=-^+a, or 4a (a - y) = or^, 

the equation to a parabola, vertex A^ focus Sj of which 
curve it is a well known property, that the perpc^ndicular 
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from the focus intersects the tangent »t any point, in the 
tangent at the vertex. 

This result being obtained by exactly the same process 
as the singular solution was obtained, of course coincides with 
it; hence it appears that the singular solution belongs to 
the curve which touches the family of curves resulting from 
the complete integral by making the arbitrary constant assume 
all possible values. The conclusions arrived at in this par- 
ticular instance, we shall now shew to hold generally. 

39. Having given the complete integral of a differential 
equation, to find its ?ingular solution. 

Let f{cOy y, d^y) s be a proposed differential equation, 
and suppose it to result from the elimination of the arbitrary 
constant c, between the equation F (^, y, c) « 0, and its 
immediately derived equation M + Nd^y ^ 0. Now change 
c into c\ any function of a^ and y, then our equation becomes 
jP(a?, y, c) =«0; and its immediately derived equation 

Sr ^N'd^y-\-Cd,c^O (1), 

c entering into M* and JV' just as c did into Jf and Ny 
and dgC of course denoting d^^^c + d^^^c' . d^y ; now C is the 
differential coefficient of F(a?, y, c') with respect to c', regard- 
ing w and y as constant, and will therefore usually involve 
0?, y, and c ; and if put equal to zero, will give such a 
value for c as makes the last term of equation (l) disappear ; 
and then the elimination of c' must evidently produce the 
proposed equation /(a?, y, d,y) = 0. Let this value of c' be 
substituted in F (a?, y, <?') = ; then this equation is changed 
into (f> (cr, y) « 0, and furnishes a relation between w and y 
which satisfies the equation f{wy y, d^y) »= 0, but contains no 
arbitrary constant, and is not deducible from the complete 
integral by giving a particular value to the constant; since 
it results from the complete integral by substituting for c 
a variable value deduced from the equation d^F(w, y, c) = 0. 
Consequently, the relation (^, y) » is the singular solution 
required. 
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40. To explain the geometrical signification of the sin- 
gular solution of a differential equation. 

Let jF(a?, y, c) = (l) be the complete integral of a dif- 
ferential equation between two variables; if we differentiate 
it with regard to c, we have dcF(w, y> c) « (2); and if 
between these equations we eliminate c, we get 0(a?, y)=0 (3), a 
where c does not appear, and which is a singular solution of the 
differential equation, of which (1) is the complete primitive. 
Suppose equation (l) to be the equation to a system of curves, 
in which the position and dimensions of any particular curve 
is defined by a particular value of the parameter c ; also let 
equation (3) be the equation to a curve referred to the same 
co-ordinate axes. Now equations (l) and (2), when e receives 
a certain value, are satisfied by the same values of w and y« 
Hence from the manner of its formation, equation (3) is satis- 
fied by the same values ; or the curves which are represented 
by (3) and (l), with a particular value of c, have a common 
point. But equations (3) and (l), being each a solution of 
the same differential equation, furnish the same value of d^y 
for the same values of a? and y; consequently the curves 
touch one another at their common point. The same thing 
happens for every one of the system of curves which equation 
(1) represents. Therefore the curve represented by the sin- 
gular solution touches in a point every curve represented by 
the complete primitive. This is the geometrical interpret- 
ation of the singular solution of a differential equation of the 
first order. 

41. Having given a solution of a differential equation, to 
find whether it is included in the complete integral or not. 

Let d^y ==/(a?, y) be the proposed differential equation, 
and y ^ F ^Wy c) its complete integral, c being the arbitrary 
constant ; and when c s c\ let this become y^u^u contain-* 
ing DO arbitrary constant ; then y^u is a particular integral 
of the proposed. 

Hence, since jF (a?, c) ^u becomes zero when c = c\ we 
have 

F (a?, c) - t« = (c - c')'". z = a%i suppose, 
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X being a function of x and c, which is neither infinite nor 
zero, when c^c\ or when « = 0; and m expressing the highest 
power of c - c' which enters into every term of F (d?, c) - w. 
Consequently the complete integral becomes 

which being substituted in the proposed equation, 

gives d^u + adgX -/{ofy u + ax) (l). 

Now since x is neither infinite nor zero when asO, we 
may expand it in a series of ascending powers of d, in the 
form 

« = JT + ^a« + BaP + &c., 

a, jS, &c. being increasing and positive, and iT, A, B^ &c. 
functions of /r; 

.-. d,tt + adgX = d,w + dgKa + d^Aa^"^^ + 8^c. 

Again, the development of /(a?, u-k-ax) will be of the 
form 

/(d7, w + a») «/(a?, w) + Jf (a»)*" + N(axy + &c., 

m, n, &c. being increasing and positive. Hence, by substi- 
tution in equation (l), observing that d^u =^f{pSi t«), we get 

d,jr.o+d,-4.a«+* + &c. 

= MaT" (JT + ^a« + Sec)"* + iVo" {K + Aaf' + &c.)" + &c. 

Now unless this equation is identical, y = u cannot result 
from y ^fiue^ c) by changing c into c ; and .•. y = w cannot 
be a particular integral of the proposed; and if it satisfies 
the proposed, it must be a singular solution. Now the indices 
w, w, &c., are known, for they result from writing u^ax 
for y in /(a?, y), and expanding according to power of a»; 
and we must endeavour to determine a, /3, &c. so that the 
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two mei&bers of the equation may be identical. If fi» be > 1, 
this can easily be effected ; for we must have c2« JT « 0, or 
X ^ a, constant ; a + I ^ m^ and d^A ^ MK^ ; and so on 
for the other terms. Consequently it will be possible to make 
the two members identical, and y^u will be a particular in^ 
tegral. In the same way the identity may be established if 
m = 1. But if m< 1, there is no term on the first side cor- 
responding to JI/JT^a"; and since JT cannot be equal to zero, 
it is impossible to satisfy the identity ; and therefore y ^uii 9l 
singular solution. Hence to discover whether a given solution, 
y ^u^ of a differential equation d^y =/(^9 y)^ is a singular 
solution or not ; we must write u -^h for y in the value of 
dj^, and if the expansion in ascending powers of h involve 
a power of hy whose index is < 1, the solution in question is 
a singular solution. 

42. To deduce the singular solutions from the differential 
equation, without knowing its complete primitive. 

Let y ^u he a, singular solution of the equation 

then by the preceding article, substituting u + A for y, we 
get 

f{pDy u + h) ^f{aiy u) + Mh^ + ATA" + &c. 

where fit, n, &c. are proper fractions ; 

••• d^f(wj u + h) = dj^f(w, u + A) = wJlf A*~^ + nNK"-^ + &c.; 

consequently, when A=:0, d«/(^, w)=oo. 

But d^f(iVj u) is what dyf{aiy y) becomes when y^u; 
and therefore, conversely, every value, u of y, which satisfies 
d^y =/(a7, y), and makes dj^f(w, 9) » Co, is a singular solu- 
tion of the equation d^y =/(a7, y). 

4f3. It is not essential to give the equation the explicit 
form d^y-f(w,y). For let d^y = p^ and let F=0 be the 
given relation between <r, y and p ; then we may regard p 
as a function of <r and y determined by the equation FsO; 
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Hence the condition dyf{x^y)^^ is equivalent to 
rf^^jFsO, provided d^^^y remains finite; and therefore the 
singular solutions of the equation V - F {w^ y^ p) = 0, are 
determined by eliminating p between r«0 and d^^^r=0, 
provided always that these solutions do not make d^jF^O. 
It is evident that if we consider y as the independent vari- 
able, and put the equation under the form 

r= F {w, y, d^ai) ^ 0, 

the same reasonings would shew that singular solutions may 
be obtained by eliminating p = d^Xy between F = a and 
dyr=0. 

Ex. 1. To find the singular solution of 
y — ivd^y + w — , — = a. 



becomes 



Here rf^)F'= - jj + — = o ; .-. P^-^y and the proposed 



(af + y-a)p = y + ^pp*, or (a? + y -a)jo = 2y, 
or (ct? + y - ay = 4«ry, 
which may be reduced to the form v ^ + y/y » \/a. 
Ex. 2. To find the singular solution of 

(y - wd^y) Lv - — j « aK 4>wy = a^ ' - 
Ex* 3« To find the singular solution of 
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These three examples determine respectively the curves 
which have the properties that OT -^OT" is constant, that 
area of triangle TOT" is constant, and that TT" is constant, 
TT' bdug the tangent at any point meeting the axes Oofj Oy 
in T and T. (fig. 5.) 

Ex. 4. y + l^p - tan-* (d,y)} rf,y - 1 « 0. 



Here d.^V « a? - tan-^p - 7-^-3 « ; 



1+y 



a? = tan~*p + 






ar - COS-' V^ + \/y - J/*, 



which represents a cycloid whose base coincides with the axis 
of ay the origin being in the centre of the base. This is 
the solution of the problem to find a curve always touched 
by the same diameter of a circle rolling along a straight 
line. 

Ex. 5. rf,y 



.*. the relation ^+y'-a*sO makes d^p infinite, and satisfies 
the proposed equation; it is <x)nsequently a singular solution 
of the proposed. 

Ex. 6. To find whether or not o^+y^-a'aO is a sin* 

gular solution of rf^y- — ^ . The solution 

V^* + y* - a» - y 

gives y = \/d* - ^, therefore substituting \/a*-a;* + A for 
y in the value of d^y, we get 

^ ^cs . \/iA 

when developed according to powers of A ; and as the index 
of A is a proper fraction, ^H-y^-^K'^Oisa singular solution, 

7 
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44. Every factor proper to make a proposed differential 
equation integrable, is made infinite by the singular solution. 

Let d^y +/(a7, y)^Ohe the proposed equation, F(47, y)^c 
its complete integral, and z the factor which makes it in- 
tegrable, so that 

« {^mV + /(^j y)} * d,F{^j y) ; 

also, let y ssu he the singular solution ; then this is not 
deducible from the complete integral, and therefore if u be 
written for y in F(wj y)^ the result will not be constant; if 
therefore we substitute u for y, in the preceding equation, 
since the second member will have a finite value, and the 
factor d^y + /(a?, y) of the first member will be zero, the 
value of z corresponding to this substitution must be in- 
finite. 

This property will sometimes lead to the discovery of the 
factor which makes an equation integrable ; as in the example 



(a* - a/^) d,y + afy = a v^ + p-^j 

a singular solution of which is a?* + y* - a^ = ; if we try a 
factor of the form 

we arrive atms— 1, ne-j^; and the factor which maken 
the proposed integrable is 
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SECTION IV. 

DIFFERENTIAL EQUATIONS OF THE SECOND OKDER, AND OF 
HIGHER ORDERS. 



45. Every differential equation of the n^ order admits 
of a primitive with n arbitrary constants. 

Let f(wy y, c,, Cg, ... c«) = be an equation between the 
variables <r and y^ containing n constants c^, e^, •.. c.. Let 
the first n derived equations be 

/i(^5 y^ d,y, c,,,.. c„) = 0, 

/2(a?, y, d^y, dly, c,, ... c„) = 0, 



/.(a?, y, d,y, dly, ... d>, Ci, ... c„) = 0. 

Between these n equations and the original, the n con- 
stants may be eliminated, and the result will be 

^(^, », dsy, dly, ... d^y) = 0, (l) 

a differential equation in which none of the constants enter. 

Conversely, a differential equation of the «* order being 
proposed, it must admit of a primitive containing n arbitrary 
constants, because this number of constants, and no more, can 
be eliminated in its formation. Hence every differential equa- 
tion of the n^^ order admits of a primitive containing n arbitrary 
constants. 

46. Again, between the original equation and its first 
n - 1 derived equations, n - 1 of the constants may be eli- 
minated, and a differential equation of the (n - 1)^ order with 
one constant will result. 
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Every such differential equation, having the same pri- 
mitive with equation (l), is a first integral of that equation ; 
hence a differential equation of the n^ order has n first in- 
tegrals, each a differential equation of the (n - 1)^ order, 
and containing one constant. 

Also, between the original equation and the first (r) of its 
derived equations, (r) of the constants may be eliminated, and 
a differential equation of the r^^ order containing n — r con- 
stants, will result, which is an integral of equation (l). Now 
(r) constants can be eliminated in a number of ways equal to 
the number of combinations of n things taken r together, or 

-A n(n-l) ,„(n-r + l) 

This then is the number of integrals of equation (l) of 
the (w - r)^ order, each a differential equation of the r^ 
order, and containing n -r constants. ' 

47* Of the general equation of the second order 

we shall first of all consider the following particular cases, in 
which dtly is involved with only one, or two, of the other 
quantities ^, 9, d,y ; and which admit of integration, or rather 
of reduction to forms of the first order. ' 

I. F(dlyf d?) = 0. Let this by resolution, give 

dly-fi^h ••• d,y^fsno^) + C; 

and integrating again, and adding another constant, we obtain 
the complete integral. The same process applies to d!^y^f(^)» 
Also, if we have F(d]J^*Sf, c^y) = 0, and put d^^y^Uj we 
get F (d^Ui, f^) B ; and if this can be integrated, and givea 
u^ f{ai)y it is reduced to the case just noticed. 

.'. y « a log - + Cj? + C 
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II. -F(djy, d,y)taO. Let this by resolution, give 
djy " /(^#y)> or ^mP - /(p)j putting d,y = p ; 



1 



{J /p B — 1 — jff ws I : 

p must be eliminated between these two equations. 

Ex. ad;y « (i,y ; .•• d,p - - ; .". logp « - + C, 

a a 

and d|,y « d,y.dpar«p.-« a; .-. y«o]!> + C'; 

P 

.-. -+C = log^ . 

a a 

III. F(<^y, y) = 0. Let this give djy «/(y); 
.-. 2d,ydiy-2/(y)d,», and (d,»V=C + 2jt/(y); 

*••• d*y « V^cTiX/W. and ^«/ ^— — ■ 

Also F(dJ+*y, d]Jy) =0, putting dT^y ^u, becomes 

and this, treated as above, gives x ^ tp {u) ; and if ^ a (ti) 
can be solved with respect to u, it is brought under Case I. 

Ex. 1. aM*y + y = 0; .-. 2o*d,yd*y + 2»d,y « 0; 

. , y ^/ 

.-. a sm"^ —7= Bs w -^ C . 
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Ex. 2. v^djy=l- 

IV. F(djy, d,y,^) = 0. 

This becomes of the first order in p and dr, by putting 
p for (i,y and d,p for dj y ; let its integral be (a?, p, C) = 0. 

If this by resolution give p or d,y =/(^)> then y = ^/(a?); 
among other cases, this will happen when the proposed is 
of the form djy + Pd,y - Q ; for it can be solved as a linear 
equation of the first order. 

If it gives (B =/(p), then 

y^ Lp^^p- f»^^sP = ^p - fpf(p) 5 

and p must be eliminated between these equations. 

Ex.1. Cy + -d,y = 0; .'. ^djy + d,y = d,(a?d,y) = 0; 

..-. ofd^y = C ; .-. y = Clog a? + C. 
Ex. 2. (1 + a^) cCy + 1 + (d,y)' = ; 

.-. — *:?- + = 0, tan •"*!) + tan"* «r = tan"' c; 

1 + p8 ^ 1 + ^p2 * ^ 

c~a? 1 /I +c^ 






y = -^logO +ca?)-- + C. 



E,.3. il±§=l«' ./(-). .-. *^ 



^ •'^••" ■■ 0T?jl /(..)• 
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/— — ; = (77-^"^* suppose; 



I' ">«* y-l:j^ 



This is the solution of the inverse problem of the radius 
of curvature, in which it is required to find the curve whose 
radius of curvature is any given function of the abscissa. 

Ex. 4. «?y + (6* - 1) rf,y ^e^y y - e* + Ce'"^ + C 

V. F{dlly,d,y,y)^0. 

Putting d^y ^py we get djy ■= dyP . d,y = pd^p ; and the 
substitution of pd^p for d^y, and of p for d^y^ will make 
the proposed of the first order in p and y ; let its integral be 

(py y> = 0. 

If this by resolution give p ^ d^y-f(y)y then 

r ^ 
V(y) 

If it gives y9Sif(p)y then 

J.p p Kf '^ p K p' 

and p must be eliminated between these equations.. 
Ex.1. ydJy-^id^yY^o, 
or d,(yd,y) = 0; .-. yd;,y^C\ r, y^^^Cw^tC. 
Ex. 2. yd^y + « (d,y)* + n = 0. 

•'• ypdyP + np^ Jfn^O\ 
.-. ^log(l+pO+»logy==^logC; 
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.-. (l+p»)s^ = C; 

.•. pa V^Cy"*"-l, and or* /--===. 

This is the solution of the problem in which it is required 
to find a curve whose radius of curvature shall vary as its 
normal; for this ' condition gives 

y\/lT7-^^^|^, or yd,p + (Ti»)(H-li*)«0, 

- or + according as the curve is convex or concave to the 
axis of <v. If n « 1 the curve is a circle, if n « 2 a cycloid, 
if ns - 1 a common catenary. 

Ex. 3. l + (d,yy'-2ydljfmSay{l'¥(d^yy}i. 

X . ayi+C 

this gives dyW* 



y/yiayi^cy 



Ex. 4. !/*d;y+\/l + (d,y)*=0, gives dyW^ ^—L === . 

VI. F(rfjy, y, a?) = 0. 

As there is no substitution by which this can be generally 
reduced to an equation of the first order between two variables, 
the artifice to be* employed in any case will depend upon the 
nature of the example proposed. Among other substitutions 
for J^yj the two following may be noticed, 

afdly-djfla^dg. -U and ^d^y »dS(^y) - 2d,y. 

Ex. 1. af^dly^2y; 

.-. dH^y)-2d,y = !*', 

Of 

or cP (wy) = - d, (wy) ; 
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Ex.2. (a;» + ff»)»<ey + a*y"0, 



«.-»-i» 



or 






let ^««; 



/ ^ 

1 +«* 



or 






Ex. 5. d^y 8 007+ 6y; 

or . dl {aof + 6y) = 6 {am + Ay), 

which becomes dP^z^hx^ putting aw -\-by ^ z^ and so falls 
under Case III. 

48. When, the equation 

F{wj yy d^yy d*y) = 0, or F{xy y, p, d^p) = 0, 

is homogeneous, reckoning the dimensions of p and d^p to 
be and - 1 respectively ; it may be reduced to an equation 
of the first order by putting 



y = w«y and d^p = - . 



8 
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For each term, if r denote its dimensions, will consist of 
some function of p multiplied by a factor of the form 



g)"^-(d,p)% 



m and n being any numbers from to oo ; therefore upon 
making the substitutions stated above, every term will be 
divisible by or'', and the equation will assume the form 

F{jpy 9, «) = 0; 

and if this can be solved relative to 9, we shall have 

d p 
g « (sf, p) ; but q « wd^^p -= (jt> - «) ^ = (p - «?) d^p ; 

• • («> ?) + («- P) d^P = 0; 
let this give 

p « ^ («), then ^ (») « sf -f wd,Xy 

which will give the required integral 



Ex. n|l + (d,y)«}»=:d;»\^^nr7; 

putting d^y « p, y ^wzy d^p = - , we find 

w 

n(l+p»)»«g\/i +«•; 
which is the same as Ex. 5. Art. 18 
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Linear Equations of the Second and Higher Orders. 
49. The linear equation of the w*** order is, 

all the coefficients being functions of Wy and each term of the 
first member involving either y, or one of its di£Ferential co- 
efficients, in the first power. The first step towards the 
integration of this equation is the establishment of the following 
theorem. 

If there be n particular values u^ ^> ^ ••• ^n) functions 
of Wf which, when substituted for y, satisfy the equation 

a^y + pidl'^y + Pidl'^'y + ... + P«» - 0, 
its complete integral is 

y = ©!«, + 02«2 + (hfh + ... + «•«!•» 
ai, Os, ••« a, bdng arbitrary constants. 

For let this value of y be substituted in the expression 

and it becomes 

d* (aiMi + Ostia -f ... + a.w.) 

or, collecting the terms multiplied by the factors a^ a,, ... a«, 

+ 02 1^«*« + Pid^'^Ui + ... + p«w»} + ... 

+ a, {(?«, 4- p,<C-*w« + ... p,tt„} ; 
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Now, since Uj, t«„ ... t«„ satisfy the equation, each of the 
quantities within brackets is equal to zero, and therefore the 
whole is identically zero ; and therefore the assumed value of 
y satisfies the equation ; and it contains n arbitrary constants, ^ 
consequently it is the complete integral of the equation. 

Thus the equation c^y - n^y = is satisfied by y = e"*, 
and since it is not altered by changing the sign of n, it is 
fllso satisfied by y = e""'; .*. y = OiC"* + agC"** is the com^ 
plete integral 

Linear Equations with constant Coefficients. 
60. To integrate the equation 

d;y + Pi<C"*yH-;>2d!J"*y+ — +p«y = «^ 

all the coefficients and. u being constant. 

First write y + — instead of y, and the equation becomes 

«Cy + Pi<ir^» + -. + l>«y = 0- 

Let y = c"*, .*. €"'(m* + p,i»*~* -f ... + p«) is Jthe value 
of the first member; now this will vanish if m be any root 
of the equation 

f»" + pi w*"* + pgW""* H- ..• + p, = 0, or ilf = /(m) = 0, 

which is called the auxiliary equation of the proposed linear 
equation. 

Hence the n real or imaginary roots of this equation, 
«9»i9 f?i2) ^ •«• ^119 provided they be all unequal, will give n 
different particular values of y, «"*»', e"**, ... c**', which satisfy 
the proposed equation ; and therefore its complete integral is 

y « OjC*!* H- OgC*** + ... 4- One"**. 

51. But if any of the roots are equal to one another, 
as, for instance, wii = m^, the value of y becomes 
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which contains only n - 1 arbitrary constants (because Ui -f a^ 
can be reckoned only as a single constant), and therefore 
cannot be the complete integral of the proposed equation. 
In this case, in order to discover the complete integral, first 
suppose the two roots mi, fTts, to be only very nearly equal 
to one another, so that Wg = mj + A, where A is a very small 
known quantity ; then the part of the value of y correspond- 
ing to these roots is 

• 1 r f^^ h sf 
OiC"'»'-fa8(e*«'«= c^J'Cai+age*') = c'"»'(oi +a2+«2 — + (h + «c.) 

a^ ^ 

= c*"i' (cj + c^x -f Cg A — + &c.), 

replacing the constants Oi + ag, and a^hy by Cj, Cj, respectively. 
Now let A = 0, then this becomes e"»' (cj -f c^^) ; and the 
complete integral consequently is 

y =. (c, -f C8^)e"»' + CgC*** + ... + c^^**. 

62. Generally, if we suppose r roots of the auxiliary 
equation to be nearly equal to one another, and therefore 
to be represented by 

nij + Aj, i9ii + A29 ... n»j + A,., where A^As, ...&r 

are very small quantities, the complete integral takes the 
form 

y = e"i' |aie*»' + <h^*' + ... + Or^*''} + ^r+i^"*^'' + &c., 
or, expanding e*»', e***, &c., 

y««"iM2(a) + S(aA).^p + ^^^^^+... 

1 • 2 



r-r [r 



+ ar+i^'^»'+.... 
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or, replacing the constants 2(a), 2 (a A), ... — , by 



r- 1 



y s c"i'(Cj + 02^?+ ... + CrX""^ + terms involving A^, A,, &c.) 

Now let Aj = = Ag = ... = A„ in which case the auxiliary 
equation has r roots each » ^i ; then the solution becomes 

y » e"i' (ci + Cga? + c^a^ + ... c^-a?'-^) + c,+ie"-+i' + ... + c.c"*', 

which contains n arbitrary constants, and is consequently the 
complete integral. 

53. Of the correctness of the above modification for the 
case of equal roots, we may assure ourselves by the following 
reverse process. Let y «= ^'w, then since 

cC {uv) = cC<> . « + ni^l'^v . d,w + ^ ^"" ^ d;"*v • dlu + &c. 

1.2 

cC(e^tt) = 6*^ {w\ w + wm»-^ d, «* + ^^^^^.^ w"-»d;w + &^ 

= e"' (to + d,)" «, 
separating the symbol of operation from that of quantity. 

Hence the first member of the equation becomes 
c""(mH-d,)*w+Pie"*(TO+d,)""''w+p2e'"'(TO+d,)*-*w+&c.+p»e"'w, 
= e""/(m + d,) u 

Now let TO = Wi, then since /(to) « has r roots equal to 
TOj, each of the quantities f(fn)if{m)y .../^*'"^^(to) becomes 
equal to zero, and the first r 4;erms vanish; and if 

w = Oi + a^w + OgO;* + .4. + Ufixf"^^ 
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then d;^«, cC*"***, •-. dlu likewise become zero, and all the 
remaining terms vanish. Hence for each group of r equal 
roots in the auxiliary equation, there will be in the value of 
y a term of the form 

e*i*(Oi + Ogd? + a^w^ + ... + a^tr*""^). 

54. Let h ^k v — 1 be a pair of imaginary roots in 
the auxiliary equation, then the corresponding terms in the 
value of y will be 

or e'*" \C (cos kw + \/^l. sin kw) + C'(cosA;a? - yCTi sinA;^)} 
= e^{ciCOskw+c^8inkw^^ (changing the arbitrary constants); 

or = v^cfTc|e'^|— ?==L===co8 kw + . ^ sin kaf\ 
Ivof + cJ vq + c| j 

= /3e*'cos(Ap<r + a), 
again changing the constants by putting 

\/c? + c| = fl, and - tan a = — . 

And if there should be r pairs of imaginary roots in the 
auxiliary equation, equal to h ^ k y/- 1, the corresponding 
part of the value of y will be 

c*'(ai + ©2^ + ... Or^'""*) (cos Ap/t + \/^ sin fed?) 
+ e*' (6i + fegO? + ... 6,0?'"^) (cos/p« - \/- 1 sinkw)^ 
or, changing the arbitrary constants, 

e**(ai + a2«» + ••• ara?*""^) cosAca? 
+ e*'(/3, + /32d? + ... (ir^'-') sinkw. 
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55. We shall now give some examples of integrating 
linear equations with constant coefficients. 

Ex. 1. djy + 4y-2y = 0. 

Let y = c"% .-. m*«"' + meT* - 2e"*' « 0, 

or m* + m-2 = 0; .-. m = 1, or -2; 

Ex. 2. a*e^*y = dly^ y = CjC* + C^e ^ + Cga? + C4. 

Ex. 3. . djy + n^y •= 0. 

Let y = e"*' ; .*. wi^ + n^ = 0, or m « 4s n\/- 1; 

= Ci (cos «ti? + \/- 1 sin wa?) 

+ Cg (cos na? - \/- 1 sin »a?) 
= (Ci + Cjj) cos wo? + (Ci - Cg) \/— 1 sinwo? 
-aicosfiw+n^ sinnofj (changing the constants) 
= )3 (cosno? cosa - sinner sin a) ^ (icos {nx + a). 

This equation, djy + n^y = 0, pf which the solution ia^ 
y = a, cosnw + a^ sin »«!i?, or y = )3 cos (w«r + a), 
is of very frequent occurrence. Hence also the solution of 

Ex. 4. dly + fi?y + azr? + 6 = 0, or 

/ aw b\ ( aw b\ 

y + -T + -l = pcos {nw + o). 
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Ex. 5. djy - Qmd^y + (m^ + n*) y = ; 

.-. y = e**' (C cos no? + C' sin n/r). 
Ex. 6. c^y + 2me^«y + n^y = 0; 

y « e""' )3 cos (^ y/n^ - m* + a), (« > w). 

Ex. 7. <y— 2d^y + 2d^y-2d,y + y = 0; 
the equation for determining m is 

(w^ + 1) (m - 1)* = ; 
y = e* (ttj + ajiX) + /3 cos (^ + a). 
Ex. 8. d*y + ^n^d^^y + w*y = ; 

- y =s (a + OjO?) coswa? + (6 + fej^?) sin »a?. 

Ex. 9. d*y + 4mdjy + 2 (»* + 3m*) d^y 

+ 4m (w* + n*) d;,y + (m* + »*)* y « ; 
the equation for determining *r, when y = ^% is 

.'. y a e-«* (a + a^of) cos w^ + €""*' (6 + fci^?) sin »a?. 

Ex. 10. d*y - 4djy + Gdlj^ - 4d,y + y = 0; 
y = e* (^ + B^ + Co?* + Da;^). 

66. To integrate the linear equation of the n^ order 
with variable coefficients, 

^iy + Pidr'y + Pa^'^y — +p«y = J^- 

9 ' 
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Let Ui^ U29 ••• »n be the n particular values of u which 
satisfy the equation 

d;w+Pid]J"*w+P2d;"**w+ ... +p„« = 05 (1) 

which coincides with the proposed if we suppose its second 
member to become zero; then, as already proved, 

If we now divide both sides by u^ and differentiate, we 
shall eliminate c^ ; next dividing both sides by the coefficient 
of Cg, which suppose v^, and differentiating, we shall eliminate 
Cg; again dividing by w^, the coefficient of Cg, and differen- 
tiating, we shall eliminate C3 ; and proceeding in this manner 
till all the constants are eliminated, our final result will be 
of the form (where each d, affects the whole of the expression 
after it) 

,1,1 J 1 J «* 

a„ dg . . . o, — dg — =s 0, 

in which, the coefficient of djw is evidently ; 

therefore dividing by this, so as to make the coefficient of 
d^u unity, we get 

which must be equivalent to the first member of equation (I). 
Therefore the same expression, only with y instead of u, 
must be equivalent to the first member of the proposed 
equation, and therefore equal to JT. Hence equating these 
equals, and reversing the operations, we find the integral of 
the proposed equation, 

y = «*iMi;»2 — f-^n-if—— — — -; 

(each symbol of integration affecting the whole of the ex- 
pression whidi follows it), which is a general formula for 
deducing the solution of a complete linear equation of th^ 
n^^^ order, from the solution of the same equation when the 
term independent of y h supposed to become zero. 



Digitized by VjOOQ IC 



67 

67. We shall now give some applications of this method. 
Ex. 1. dj jf + 4d,y + 4y = aa?. 

If we suppose u to denote the value of y in this equation 
when a = 0, since the roots of the auxiliary equation are 
-2, -2, 

in which expression the coefficient of d\u is e'*', therefore 
dividing by this coefficient, and changing u into y, we get 

.-. y = e"** j[^* {asd^^') ^ ay— + f ) + ^"*' (^i + ^a^)- 

Ex. 2. djy - 2wd,y + (m* + n*) y = X. 

If we suppose t« to be the value of y in this equation 
when Jr=sO, we have 

u = C^d^ cohnw + Cgc'^sin Tia?; 

.'. dividing by e^'cosn^, and differentiating, we find 

d^e"^* sec nx,u— C^n sec* nw ; 

. •. rf, cos* w a? d^c""** sec Tia? . w =. 0, 

in which expression the coefficient of dlu is e""' cos n^, there- 
fore dividing by this coefficient and changing u into y, we 

c"* sec »a?d,cos* nwd^e'"^' sec nw .y ^ X\ 

,\ y = c"'cosn<!Pjtsec*wa?jie~"'* cosw«j»-y, 

or, if we suppose the constants to be added after each inte- 
gration, 

y=e'«*(C' cos wa? + Csin ncT?) +e"" cos Ti^i' j, sec^w^j^c'""cos w 07^. 
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Hence the solution of djy + «*y = JT is 

y =s CcoswcTH- C AnfKv + cos nw fgSec^nxfgCOsnxX. 

Ex, 3. Let JT = ^ cos (mof + a) + B cos («« + fi), 
then fgCOsnwX 

= — T -_ tm sin (wo? h- a) cos na? - n cos (mo? + a) sin Wtp| 

nr -rr ^ 

+ — I — sin(2w/r + )3) + a?cos)3>; 

J co8(fii^ + a) 

••. Lsec*wa? LCOSW^-A = r r 

•^ •' mT'-nr' cos «d? 

B fa? sin (nw + fl) sin fl ^ I 

+ — < ^: ^-^ ^tan»d7>; 

2 [ ncosTi/r 2nr j 

.\ y =2 C cos wa? + C sin »a? + -^ cos {mw + a) 

+ — w sin («^ + i3). 
2n '^^ 

Ex. 4- d^y - lOdJy + 62d^y - 210d,y + 26ly = e*. 

The roots of the auxiliary equation are 

2 ± 5 \/^, 3, 8 ; 

.*. t< s (a + aiO?) e^ + e** (6 cos 5^ + c sin 5/i?) ; 

.-. d*(e"*'w) = -c"'{(246+ 10c)cos5i? + (24c + 106) sin 5^}; 

.\ d, cos*5^d,e'sec5a?d'(6"^w) =0; 

.'. c^' sec 5w d, cos* Swd^ef^ sec 5/i?dJ («"^y) = e* ; 

- cos 5^7+ 5 sin 50? 



.-. cos* 5«d,( ) = ^-*- 



26 
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26 \ cos 507 COST 5a?/ 



1 C-" 



.'. e* sec 5wdi (c^y) = —■ ; 



*v ^^^ 26 / ^ 104 
y = — - + (a + ajO?) e*' + e** (6 cos 5^7 + c sin 5w). 



Ex. 5. djy + p^djy + p^cgy + p^d^y + ^4^ « ^; 
where the coefficients are such that the auxiliary equation 

^ + Pi^ + P8** + j»s* + P4 = K* + ^y + ^'}* = ^» 

has two pairs of imaginary roots. Then the solution of 
d^u + Pidjw + p^dlu + ps^^,!^ + P4«^ = 0, is 
€^*u = (a + a} Of) cos »^ + (6 + b^a) sin «d7 ; 

.'.. d^(secnw€^'u) = a^ + 6* tanwo? + «(6 + 6**17) sec'n^^ 

a* 6* . 

.•• cos^nwdgi ) = — (14- co8 2n^) + — sin29»j7 
2 ^2 

+ n(6 + 6*a?); 
••• d, ( ) = — na^ sin 2wa? + nb^ (l + cos 

.-. sec?nafd:,( ) = - 2«a* tan«47 + 2w6^; 

/. dg9e(fnwd,( ) a= - 2n*o^ sec*wa?; 

.•. dg iXM? nw dg sec? nxdg( ) = 0; 
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and as the coefficient of d^u is e"' cosn^, dividing by this 
coefficient and changing u into y^ we get 

c""'8ecna?rf, cos'nozd, sec* wa?rf, cos* «a?d, sec »d76*'y = Jf ; 

or, if we suppose a constant to be added after each integration^ 
we may add to this value of c"'y, the terms 

(a + a^a) cos nw + (6 + h^a) sinn^. 

Ex. 6. Suppose the linear equation to have constant 
coefficients, and let — ai, — a^, ... -a« be the roots of its 
auxiliary equation, so that 

m* + PiWi""^ + Pt«»""* + ... + p, = (m + Oj) (to + Oji) .•. (to + a,), 
then u = Citf-'»»* + Ct6-««*+C3r'«»* + ... + c^e-*''; 
d, (c«t*tt) - o, (a, -- o,) c(«» -^t)* 

+ Cs («i - «3) e^""* " '^^' + .-. + c^ (flTi - a.) ^(«i -«•)*; 
d,tf(«a-«i)»d,(e«i*ti) «Cb(o, - oj) (a^ « a,) «(«»-«»)» + ... 

+ c. (a, - a,) (a, - a^ e<^"^>S 
&c. = &c. 

d,e(««-«-i)*d,e(^» -«-»>* ... d,c(««-«i>*d,(e«i'«) «0; 

therefore, dividing by e^*, which is the coefficient of d^u, 
and replacing u by y, we get 

e-«-'d,ef«-- «-»)*... d,e(«»-«i)»d,(e«i*y) » ^; 
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Method of Parameters. 



58. There is also another mode of integrating linear 
equations, which deserves to be mentioned, called the method 
of Parameters, which we shall now explain. It may be stated 
thus. The complete integral of the linear equation of the 
n^ order will be of the form 

y « v^u^ + v^u^ + ... + <>,ti„, 

where Uy^^ u^^ ... u^ are the n particular integrals of the 
^nation when the term independent of y becomes zero, and 
<^i, t>g, .., r, are functions of a?, determined by equations of 
the form 

Let the proposed equation be 

and suppose y = v^u^ + v^u^ + ... <>«w« = 2 (vw), 

and as we have made only one assumption respecting the n 
independent quantities ©,, u^ ... tJ«, we may make n-1 more; 

now d,y = 2 {vd^u) + 2 (wd,«) ; 

.-. d,y = 2 (vd,«), 

putting for the first of our additional assumptions 

2(«d,i>) = 0, (1); 

similarly d^y « 2(tJrfJ«^), putting 2 (d,w d,tj) = 0, (2) 
&c. « &c. 
and Cy = 2 («(C«*) + ^» putting 2 {dl'%d^v) = X^ (w). 
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Substitute these values for y, d^y^ dly^ ,^. diy^ in the 
proposed equation, and the first member becomes 

-I- &c. +«,(d;w«+Pidi;"*w«4-... + p,Wn) + Jf> 

which manifestly reduces itself to ^, since each of the quan- 
tities within brackets is equal to zero ; hence if the parameters 
^19 ^f > ••• ^11 be determined subject to the above n equations 
of condition, the assumed value of y will satisfy the equation ; 
and since there are n equations in which the n quantities 
dgV^f dgV29 ... d„v^ enter to the first degree^ each of the 
parameters v^, v^^ &c. will be determined by an equation of 
the form d,t?j »/j (ai), and so n arbitrary constahts will be 
introduced. 

Ex.1. €^y ^n^y ^ ^ecnw. • 

The solution of djy + n^y = is y^a^ cos nw + a, sin nw. 
Let ,*. y = v^cosnaf -^ v^sinnw; 

•*• d#y = "" ^1^ sin nw + v^n cos nw^ 

making d^v^ cos no? + d^^Vf sin n^ s 0. 

dj J/ = - v^n^ cosn/v - v^n^ sin nw + sec nw^ 

making - d,«j n sin nw + d^v^n cosnw== secnw; 

I sinw^ 1 , _, 

.-. d,«i = , or «! « — log cos na + C,. 

n cos nw n 

1 w 

d;pt>2 = - 5 or tjg = - + Cg r 

n n 



y B |— log cos wa? + Cj j cos wa? + | - + d] 



smno?. 
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Ex. 2. d> + Pidr'y + P9«"*y+ ... +p„y = Jr, 

where the coefficients are constant, and the auxiliary equation, 
JIf =0, has none of its roots — Hx, -^ (hi ••* -^lu equal to 
one another. 

Here tti«:c"*^', t«2 = e""»', ,.. w» = e"*^; 

and in order to determine the parameters v^9 v^9 &c. assume 



then, by substitution, the n equations of condition become 

2(c) «0, 2(ac) = 0, 2(a*c) = 0, ... 2(a"-*c) = 1; 

and from these equations, we have to determine each of the 
constants e^, c,, ...c^, by eliminating the rest; for this pur- 
pose multiply them in the reverse order by the indeterminate 
coefficients 1, fj, q^i ••* 9ii-i» and take the sum; then 

Ci «-* + ?,«""* + ?2«r' + ••. + ?.-i) 

+ Cft (aj"' + 9i«r' + 92»r' + ••• + 9-1) 

+ .•. + c,(a:"^+9i«r* + ... +?i.-i) = i- 

Now assume q^^ q^j ... 9«_i so that the coefficients of 
<^, Cs, ••• c,, may each equal zero; then 029 ^9 ••* o* a^ 
the roots of the equation 

therefore, replacing the coefficient of c^ by the product of its 
simple factors, 

e, (a, - ttg) («! - ©3) ... (Oj - a,) = 1, or c, «- j^; 

Xe^'J'; similarly, «,=— i_/,e«.«Jr, &c.; 



10 
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• '-'^{iS^H- 



This method eyidently fails when M^O has equal i^oots^ 
since one or more of the quantities d^a^J^j d«=a,-M', &c. 
in that case becomes zero. 



Ex. 3. 
Here 



diy-y^e'. 



JIf = m*-l =(iw + l)(m- l)(fii + \/^)(m-\/^) = 0; 
therefore the four values of fg^X are 



and the four values of 



1+V-l 



are 






l-V-l 



grV-l 



4 "-4 -4,y^«i 4V-1 






= ■- + c^e"' — -- + Cjje* + — + Cs cos/v + C4 sm ^ 

8 44 

= — + Cje"* + CjC* + Cs cos 07 + C4 sm ^. 

8 4 

Ex. 4. di^-^scos^; 

th6 four values of f^e^X are 

--Xcosa?4- sin^) + Ci, — (-cosd7 + sin«) + c,, 
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^ (a? + sin .re*^^) + Cj, ^ (^ + sin we"''^) + C4, 

and the values of -; are the same as in the preceding 

example ; 

1 c 1 c 

•'• y = -(cos^r + sino?) + — e"*+ -(cos a? - sin a?) — -^ 

8 4 8 4 

(jffer^^^ + sin w) y==r er'^^ 



8\/-l 4\/-l 

+ p== (/re* >^ + sin 0?) + — ^ e*^^i 

8\/-l 4V-1 

cos J? 07 sin ^ , 
= ~^ + — p- + i(cie-'- e^O 

+ — 7=(cosa? + v/^ sin a?) ^=(cos<r-\/-l sin a?) 

4v -1 4v -1 

07 sin or 
= : -|- «!«• + €L^e~* + Oscoso? + a* smo?. 



59. When the auxiliary equation contains equal rpots, 
we may, by the following theorem, depress the proposed 
equation to one whose auxiliary equation shall contain only 
unequal roots, and then complete the solution by variation 
of parameters. 

If djy + Pi<C"V + P2^ *y + ••• +P«y = -y« be a linear 
equation of the n^ order, with constant coefficients; and 
— a^, — «2> •••—«« the roots of its auxiliary equation, so 
that 

then if we multiply both sides of the proposed by e"*, - a 
being any root, and integrate, the result will be a linear 
equation of the (n - 1)*** order with the same auxiliary equa- 
tion, except that it wants the factor m + a. 
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For let w»-^ + ^j^"-* + ... + 9,_i -=0 be this auxiliary 
equation. 

Now d^e^ (dr'y + 9i<C"*y + .- + ?»-i») 

« <r (djy + Pid:;"'» + P2<"*y + ... +p«») « «^^«; 

the auxiliary equation being m*"* + q^m*~^ + ... + q'.-i = 0, 
cleared of the factor m + a ; and if this factor should occur 
r times, then repeating this process r times, we should get 

drry + t,drr''y + ... + t.^^y = e'^'i:e^x^, 

whose auxiliary equation differs from JIf = o, only in wanting 
the factor (m + a)*". 

60. If we know a particular integral of a linear equation 
of any order that has no term independent of y, we may 
reduce it to another equation of the same kind, of the order 
immediately inferior. 

Let y B t^ be a particular integral of /(d,) y = 0, 

or <y +pid;;"'y + P2rfr'y + — +Pi.y = o, (0 

and assume y^uf^x^ % being a new variable, a function 
of w ; then, separating the symbols of operation from those 
of quantity, we have 

d, being understood to affect u only, and dl to affect f,z only* 

Hence, substituting for the differential coefficients in the 
proposed, by this formula, and again separating the symbols of 
operation from those of quantity, we get successively 

(d^ + diyuf^x + pi (d^ + dlY'^uf^ss + &c. +pnU^« 
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= f,«.f(d,) » + - «/'(d,) « + — - d,zf"{d,) u + ... 

1 1 • ss 

+ — — dr*«/""' (d,) « + tfj-'sr . tt. 

W — 1 



But since t« is a particular value of y, /(d,) w = ; hence, 
reversing the order of the terms, and observing that /^'^(d,) 
denotes the same function of d,, that (F„f(v) does of t?, we 
get for the depressed equation, 

udl''^% + (nd,u +PiU)di'^« 

Similarly, if we know another particular solution, u^, of 

equation (l), then d, f — j will be a value of j^f in equation (2), 

and we may depress this equation to another of the same 
form of the {n - 2)*^ order ; and if we know r particular 
solutions of equation (l), we may in this way depress it to 
an equation of the same form of the (n - r)^ order. 

61. I{ y =^u be a particular integral of the equation 

dJy + Pd,y + Qy = 0, 

the other particular integral will he y ^u j — r— . 

For if we denote the two values of y by t^ and t?, we 
have 

d^w + Pd^u + Qt* = 0, dlv + PdgV + Qi> = ; 
►•. «d*« -t>d^«* + P(wd,t? - wd,tt) = 0; 
.-. ud^v - vdgU = Ce~'f'^; 



.-. V ^ Cu — — 
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Change of Independent Variable. 

62. Besides linear equations with constant coefficients, 
very few equations of the second and superior orders admit 
of being integrated ; and those only by particular methods. 
Sometimes the integration of an equation may be facilitated 
by changing the independent variable. 

Ex. dly-^a {d^yf •{• bw (d^yf = 0. 

Make y the independent variable, and consequently for 
^*y, d^y, write 



1 




d?^w 










dyW' 




(d,a>y 


9 










- 


id,«y 


a 

1 


bw 


= 0, 






id,a 


')* idyWf 








or 


dla- 


adyO) - bm 


= 0, 




tegrab 


le. 






t 


-- f 


^^4.1^ 



63. In the above example there is no difficulty in fixing 
upon the new independent variable. In other cases we must 
consider a and y as functions of a third variable t^ and sub- 
stitute the values of d,y, dly^ 8ec. corresponding to that 
supposition. When the equation is thus generalized, we must 
assume for co or y some known function of ^, according to the 
circumstances of the case, so that there may arise for deter- 
mining the function of t which expresses y or a?, a differential 
equation simpler than the proposed one ; between the integral 
of which and the assumed function, if we eliminate ^, we 
obtain the required relation between a? and y. 

Ex. 1. (1 -^)d^y-a?d*y + w^y = 0. 
The generalized equation is 

„, dtO} d?y -- diy dt w dty , 
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Let m = cos t^ then d|a? = - sin ty dtW^- cos t ; 
therefore, by substitution, d^y + n?y-(^j which gives 
y^ A co^nt + B sinn^; 
••. y-A cos {n co^'^w) + B sin (n cos"^ ai). 

^ A B 

Ex. 2. d!y + r- d*y + 7 r-Ti y = 0. 

The generalized equation is 

dtfvd^y-djyd^x A dty ^ B ^ ^ 

(d^ofy a + bwdtOf {a + bopy 

Let diW^a-^-bXy or 6^ = log (a + 6a?), or c*' = o + 6^; 
then d^x^bdt(X!\ 
therefore, by substitution, 

the solution of which is 

, mm* 

y ^ ci€r* + CgC*"'* = c, (a + 6a?)* + Cg (a + 65?)"^, 

j» and m' being the roots, supposed possible, of 

w* + (i< - 6) m + i? = 0. 

If the roots be impossible, and of the form m ±n\/-l, 
the solution is 

m n 

y = /3e"* cos (n^ + o) = /3 (a + bw)^ cos {log(o + 6a?)^ + o}. 
Ex. S, (a + 6j?)'^dJy+A(a + 6a?)*d;y + i?(a + 6j?)d,y + Cy 
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A« in the preceding Example, 
let a + bx ^ e^ i .\ dtOf = c*', 

.-. by substitution we get 

a linear equation with constant coefficients. 

The same substitution of course succeeds for the equation, 

(a+6a?)"d;y+J(a+6^)"""'dJ"'*y+&c.+Jr(a+6^)d,y+Ly«/(a?), 

which, by putting a + 6a? = jj?, may first of all be reduced 
to the form, 

«"d;y + Oj«"-»c5-^y + &c. + a^^^zd^y + a^y =/i («) ; 

and then to a linear equation with constant coefficients, by 
the formula 

s^dly = d?y - p^O^^^y + &c. =t j>«-iy, 
where z = e\ and pj, p^^ &c. are such that the roots of 

ft" - Pi Af -* + Pgft""* - &c. ... ± p^.iAp = 0, 
are 0, 1, 2, 3 ... (n - l). 



Ex-. 4. a?^djy + a?d,y - y : 



^-l' 



y = £^^C,.H-% ?(!-.) log (^) 
The substitution ^ = c* gives 



d?y-y = a + 



- 1 
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••• y = C,e-' + Ce* + ^ + ^ (e"' - ««) log (^) , 
which agrees with the above, upon writing a for e*. 
Ex. 5. ^djy + o^d^y - y =» a^> 



y = — ^ + C^oB + 



a? ' n*-l' 



Solutions expressed by Definite Integrals. ^ 

64. Sometimes the complete integral of a differential 
equation may be expressed by a definite integral, as in the 
following instances. 

Ex. 1. dl'^y tswy-^m. 

Let t>i = Oi / e^i*^.e «, a^ being a constant quantity ; 

then di-^v^ « a? Y"^"'^**''^"* =^*(^ + ^^i)> 

change a^ into a^, os, &c., and let v,, «8, &c. be the cor- 
responding values of Vi, 



Now suppose ttj, ttg, ...a„ to be the n*^ roots of unity; 
then multiplying each of these equations by an arbitrary 
constant, and taking the sum, 

^"^ (Cl»l + CgVg + ... + CnVn) = Ci + Cg + ... 

+ <?« + ^ (Cit?i + Cjj«, + ... + C,W,), 

let c, + Cg+ ... +c, == »i. 



11 
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then y^Ci^i + c^v^ + ... + e^^v^ is the complete integral, con- 
taining only ^ - 1 arbitrary constants by reason of the equa^ 
tion of condition ; or it may be written 

ya J «-•• {c6** + c,ac«** + c,aV'*' + — + <'-«-ia""'e«""**'}, 
a being a primitive root of Ap^ - 1 »& 0. 

Let m » 0, n = 3, then d^y^wy, the solution of which is 



I cos- 



-iU + 5Vi)m^}. 



Let m » 0, n B 2, then d«y « ^y, and 
y«C / e «(e^+«"'*)=B-y= / c* cos a? «, putting ^=»v-l, 

« -7= . ^-:- = C v^ire* , (Integ. Cal. Art. 112). 

which agrees with the result obtained by direct integration. 
If the more general form were proposed, 

assume a + j3a? « /3 • ^ ; .•. dj,y = dty dj = /3*dty, 



and cO = i3«d?-'y«/3«^y; 
.-. d^-^y^ty; 

and to this form may the still more general case 
di'^y = a J?y + 6a? + cy + c 

be reduced. (The above solution is taken from Crelle^s 
Journal, Vol. x.) ' 
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Ex.2. ci^y+ [n«-^] y = o« 



Let y^ ~ where m (m + 1) s a, and suppose a to 
be such that m is a whole number; then 

w (dj« + n*w) == ^mdgU. (l). 

Differentiate this S, 4, &c. times successively, and add 
each result to the preceding multiplied by n*j in doing which 
the symbols of operation and quantity may be conveniently 
separated, then 

ofdl (dj + n*) « + 2d, ((^ + «») f^ « ^mdluj 

and t»«* (dj + n*) t* =» 2fiin^d«t£ ; 

.-. a?(dj + n«)»w = 2fiid,(c5 + n*)tt-2d,(c^ + n«)a 

= 2(m-l)d,(c5 + n«)tt; 

similarly, tf? (c^ + n*)^ w = 2 (wi - 2) d, (dj + n^y w, 

and repeating the process m times, the second member becomes 
zero, and 

(dJ + nT'''« = 0. 

Hence we have, for the determination pf u, a linear 
equation with constant coeiBcients, whose auxiliary equation is 

.'. u «■ (a + UiW + ©2^ + ••• + «m^") COS no? 

+ (6 + 6i47 + ^^ -► ... 6a,aj^) sin«d7. 



* This equation, when a = 6, presents itself in investigations relative to the 
figure of the Earth ; the solution in the forms here noticed was given by Mr Gaskin 
(Cambridge Problems, 1839). 
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this gives the form of the solution; the arbitrary constants 
must be reduced to the proper number, in any case, by 
actual substitution ; thus let m = 2, so that a » 6 ; then 

t^ s (a + diOf + Og^ cos nof + (b + biW + b^sf) sin nor, 

and substituting this value in (1) in order to reduce the con- 
stants, we find 

a^y s t^ ss a {cos nmil J + nw sin no? J 

+ 6 {sinn^ |l J -n^vcosn^}, 

for the complete integral of 

But when in the equation m (m + 1) ^ o, m is a frac- 
tion, 

make y = z»*+^ w, then w (d/u + n*«) + 2 (i» + l) d,« « ; 
the complete integral of which is 

u^(i "TC^ - »T cos (wt + a), 
a and )3 being the arbitrary constants ; for this gives 

d^u = /3 "*r (e - T^T (- ^)-sin {mt + a) 
•'I 

<e« =/j "y " (<» - «»)■ (- /») cos (j>< + o) 

=/8 ' Y'- (<* - n*)" (<•-«' + »*) COS (/»< + a) ; 
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Similarly, when m is an integer, it will appear by sub- 
stitution that 

is the complete integral of w {dlu -^ n^u) ^ ^md^u, the re- 
duced equation when m is an integer. 



Simultaneous Equations. 

65. In these equations, which sometimes occur in the 
higher parts of Dynamics, instead of one equation between 
J7, Pf and the differential coefficients of y with respect to a, 
being given to determine the relation between w and y; we 
have two equations containing w, y, t (of which w and y are 
considered as functions) and the differential coefficients of or 
and y relative to tj to find that relation. 

66. To integrate the simultaneous equations of the first 
order, • 

aw-^ by + diW = T, a'w + b'y ^d^y^Ty 

T and T' denoting functions of t. Multiplying the latter 
by an indeterminate quantity m, and adding it to the former, 
we get 

aw •{' by -{-m {a'w + b'y) + d^ (^ + my) = T + mT^ 

,, vy /*• b -^ fnb ^ __ __, 

or dj (w + my) + (a + ma ) (cV + -, y) = T + mT. 

a + ma 

Let > = m, which will give two values of w, m, 

a + ma' 

and m, ; then the equation, under this condition, becomes a 
linear equation of the first order; and we obtain by inte- 
gration 

(a? + my)«^^«'*')* « jf (r + mr')e(*+'»«')*; 
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and by substituting successively the two values of i», we 
obtain two primitive equations, which will furnish valuer of 
X and y in terms of t^ and the relation between w and y, 
if ^ be eliminated. If the two values of m are equal, we 
shall obtain only one equation between w^ y^ and t\ but if 
this can be solved with respect to w at y^ and we substitute 
the value so found in one of the given equations, we shall 
obtain a second relation either between w and t^ or between 
y and i\ and then t may be eliminated as before. 

Ex.1. 5a - 2y + dtx := ef, 6y - w + dty ^eF^i 
(5 - m) 07 + (- 2 + 6f») y^dtiw-^ my) = e* + me**, 

let =m, ori»* + m-2 = 0, or m » 1, or -2; 

5 -m 

.-. ^ + y = 7-f -T +Ce-"; 

5 O 

^ 2 

similarly j?-2y = c** + C\e"'*, 

8 9 

which determine w and y in terms of t 

Ex.2. dtaf + 5af + y^e\ dty-^Sy-^w^^; 
.-. dj (a? + wy) + (5 - m) a? + (1 + 3 w) y = e* + we** ; 
1 + Sw 



5 — wi 



= wi, or 1 - 2m + m*= (l -m)^ = 0. 



Hence the values of m are each » l, and integrating, 
we fiDd 

^ + y = Cifj-** + -€* + -^. 

5 
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By means of this, eliminate y from the first equation^ and 
we get 

4 ^ 

5 o 



25 36 ^ 



and y«le' + lc^+C,(l+0^"'-C',c-^'- 
2o 3o 



The more general form 

aor + 6y + JdtW + i?d|y = T, a^ + b'y + -rf'dj.i? + ffdty « 7^, 

may evidently be reduced to the above by successively elimi- 
nating diyy d^w. 

67* To integrate the simultaneous equation of the second 
order, 

ao? + 6y + c + <?^ « 0, a a? + fe'y + c/ + c^y * 0* 

Multiplying the latter by an Indeterminate quantity m, 
and adding it to the former, we get 

{a + ma') ar + (6 + mb') y 4* <? + mc' + <^ (a? + my) = 0, 

or c? (^ + my + c^) + (a + ma ) {w + my + cj = p (l), 

_ 6 + m6' ,^ c-mic' 

if ,^m (2), > = Ci; 

a-^ma a + ma 

therefore, integrating equation (1), and substituting succes- 
sively the two values of m given by equation (2), we obtain 
the two required primitives ; or if the values of m be equal, 
we must proceed as in the former case. 
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Ex. 1. dfw - (Sx + 4y - 8) = O, dt*y + (w + y+s)mO. 

.-. <?(a? + my) + (m-S)(a? + -— -y+ — — -1 -0. 

- m '-' 4t . ^ vo 

Let sin, or m*- 4m + 4 = (i» — 2r« 0; 

.-. d? (a? + 2y) -r (ar + 2y - 13) = 0; 
.-. ^ + 2y - IS = ce* + c c~*, 

and eliminating a from the latter of the given equations, we 
find 

.-. y ss 18 - - (^ - ^) e* + - ^«~* + a«* + a c""*, 

and ^ = - 2S + c (^ + ^) e* - c'(^ - 1) c"* - 2ae* - 2oV*. 
Ex.2. dlaf+adtW^Oj d!^ y -k- adtff + h :» 0. 

68. If we have three variables w^ y, z^ which are func- 
tions of if and if 

d^of + aw ^ by + c X = Tj 
dty + aifl? + h^y + Ci» = Tj, 
dj« + 0,0? + 6,y + c,» = Tg, 

then if we multiply the second and third by indeterminate 
constants 97», nh\ and add them to the first, and assume 

4 + b^m + b^ml = «« (a + aitn + Ogm'), 

c + Cim + Cjw' = m' (a + OiW* + flj^O? 
we have 

di{w -^my '\- m'z) + (a + moi + m'a^ {x-^my^ m'z) 
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To determine m, m\ we obtain •cubic equations; hence 
if nil, wig, Wg, «i/, m/, wig' be their roots, we have, by 
solving the linear equation, 

from which equations a?, |^, % may be found in terms of t. 

Approximate Solutions of Differential Equations. r\-,.i \ 

K 
69. When all the, known methods of integrating a 
proposed differential equation fail, we must endeavour to 
resolve it approximately, that is, to obtain from it the value 
of y in terms of a?, in the form of a series. The first mode 
which presents itself of effecting this, is to assume for y a 
series arranged according to powers of zr?, with both its co- 
efficients and exponents undetermined; for in most cases it 
happens that the exponents do not follow the progression of 
the natural numbers, and that particular artifices are requisite 
for discovering their law. When the form of this series is 
known, we determine its coefiicients by substituting it and its 
differential coefficients, for y^ d^y^ &c. in the proposed equa- 
tion. The following application of the method to Riccati^s 
Equation will give an idea of the mode of obtaining both the 
exponents and the coefficients. 

Ex. d^y + 6y^ - aa?" = 0, 

let 6y = -^; .-. d^ss - ab% af" ^ Oy or putting c for - a6, 
dlx + cxdf = 0. 
Assume ^ = ^« (-4 + Bafi + Cw^^ + &c.) ; 
.-. a(a-l)Ja?«-^ 
+ (a+)3)(a+/3-l)Sa?«+^-*+(o+2j3)(a+2/3-l)C^«+^^-«+&c.] 



+ Ci<i?«+" + cB 

= 0. 
12 
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Hence /3 must equal w + 2 ; and then to determine 
a, -4, By &c. we have 

o (a - 1) -rf = 0, 

(a + » + 2) (a + « + 1) i? + c J = 0, 
(a + 2w + 4) (a + 2n -H 3) C + ci? «= 0, 
(a + Sw + 6) (a -H 3n + 5) D + cC = 0, 



Hence a may be either zero or unity, and A remains 
undetermined; calling therefore A and A' the two values of 
A corresponding to a = 0, a = 1, we get 

^ (n + l)(n + 2) (w+l)(n+2)(2n+3)(2n+4) ^'^ 

■*■ ^^ '"(n+2)(n-|-3) "*" (n-|-2)(n+3)(2n+4)(2n+5)"" ^'^' 

d z 
and substituting this value of « in the expression y » -^, 

we shall obtain the value of y, involving only one arbitrary 

A 
constant -^7. 
A 

As the terms of the above series have divisors of the 
forms (w + 2) i =f Ij where i is an integer; if n be such that 
(» + 2) i =f 1 = 0> one or other of the series will be illusory, 
and we shall only obtain a particular value of %; and if 
n + 2=:0, both series become infinite, but in that case the 
equation may be exactly integrated by Art. 16. 

70. In the preceding instance we arrive immediately at 
a complete result ; but it often happens that the solution we 
obtain by the method of indeterminate coefficients involves 
no arbitrary constant. To supply this defect, we must in- 
troduce, instead of the arbitrary constant, a value of y cor- 
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responding to a given value of w ; that is, supposing these 
to be 6 and a, we must substitute in the given equation, 

then determine u in a series all whose terms vanish when 
^ = ; and replace u and t by their values y — 6 and a? - a ; 
in this way it is evident that the arbitrary constant will be 
involved implicitly; for, from the complete integral 

f{x, y, C) = 0, 

C may be expressed in terms of a and 6. 

Ex. d^y-^y^gaf^, 

this becomes d^u + b + u = g{a '\- 1)"^; 
assume u^t^(A^ Bt^ + Ctr + 8ec.) ; 

+ b + J^ + Bt^-^P +&C. 

-gar - mgdr-H ^ ^ga"-^^-8EC. 

.-. a = 1, )3 = 1, 7 = 2, &c. 

A = gdr - 6, 2i? = gmoT-^ - g-a"* + 6, 

6C = ^m (w - 1) a"»-^ - gma"^-^ + ^a*" - 6, &c. 

71. The approximate solution of a differential equation 
may sometimes be obtained in the form of a continued fraction 
by assuming 

^Aw^ Bafi Cwy 

^~TT TT 1 + &c. ' 

First, suppose. 0? to be very small, and for y substitute 
Aw^ in the given equation; then, retaining only the terms of 
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lowest dimensions in w^ A and a become known by equating 

coefficients and exponents. Next, write y— in the 

proposed equation, and in the result put x = Ba^i and de- 
termine B, j3, as before, by supposing ^ to be very small ; 

Bof^ 

then in the transformed equation in », put as = ; and 

so on for the rest. 

Ex. my + (1 +a?)d,y = 0. 

Let y^Aof^; .-. (m + a)-4/i?« + -4a^"^ «0, 
or Aaaf^'^=^0; .-. a = 0, 

and A remains undetermined. Next, put 

J 
y = , and x - Bafi. 

and we get successively 

w (1 + «?) = (l + ^) dgZ ; 
»» + J5a;^(wi-/3) = i3Ba?/*-S or m^Bfiaf^-'; 

.'. /3=1, B^m; 
similarly, putting 

z as , and t = C^y, 

we determine C and 7; and thus we obtain 

m — la? m ■{- 1 Of m^2 w m+2 w m-^S w 



A mw 1 '2 S *2 3 '2 5 *2 



y=T-T7- 



1 + 1 - 1 + 1 + 1 - 1 - 1 + &c. 
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Since the proposed equation when integrated gives 
y^A^l + wy^, 
the above continued fraction is the development of A (l +^)"". 

72. We may approximate to the integral of a di£Per- 
ential equation by successive substitutions, in a manner 
similar to that invented by Newton for the solution of 
algebraical equations, as in the following instance. 

Ex. dly + n^y + aj^ + a = 0, where a is a very small 
quantity, 

we may assume 

y = « + o«i + a^Ufi + &c., 
which gives 

d^u + n^u + a + o (d^u^ + n^Ui + u^) + a* (c^t^ + n^u^ + ^uu^) 
+ a^ (d^Us + w*i^ + 2uUi) + &c. = 0. 

Hence, equating the coefficient of each power of a to 
zero, we get 

d^f* + n^u + a = 0, 

d^W2 + w^«^ + 2wWi = 0, &c. (3.) 

a 

The first give u ^ + c cos nw + c^ sin na; ; 

n? 

and this value substituted in the second reduces it to the 
form 

the integral of which by Art. (58.) is 

Ui = cos no? (c. LX^ sin »«») + sin nw (c/ + — t JTi cosniv). 

n n 
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Similarly, these values of u and u^ substituted in (3) 
reduce it to the form 

which may be in like manner integrated; and in this way 
the coefficients of the powers of a may be deduced one from 
the other by a uniform process. 

73. We have seen (Art. 57.) that the solution of 

dly + n^y = A cos (mof +a) + B cos (nof + )3), 

is y = Ci cosna? + Cg sinwo? 

+ T cos (mo? + a) + — a? sin (na + /3). 

71— m 2n 

Hence, if from the proposed equation we had to deter- 
mine y approximately, we could not neglect the term 

A cos {mw + a) 

even when A is exceedingly small, provided m and n are 
nearly equal to one another ; because in the value of y this 
term is divided by n^ - m^ which is very small. With re- 
spect to the last term in the value of y, we remark, that 
it is not periodical, but may increase indefinitely, as w in- 
creases. The equations which present themselves for solution 
in physical Astronomy are usually of the above form; and 
upon the peculiarities just noticed depend some of the most 
interesting results in that subject. 
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SECTION V. 

ON DIFFERENTIAL EQUATIONS INVOLVING TWO OB MORE 
INDEPENDENT VARIABLES. 



74. In this part of the subject, it is convenient to 
employ a system of notation which we have not hitherto found 
it necessary to notice, and which may be explained as follows. 

If t^ be a function of any number of independent variables 
^> Vi %y ty &c., and if w + 5t^ be the value of u when these 
variables simultaneously become 

X + ^Xy y -i- Sy9 ^ + S%9 &c., 

then Su = d^u.Sw + dyU.Sy + dgU.Sz + ... 

+ terms of two dimensions in ^<v, Sy, &c. 

If we BOW agree to neglect all terms of higher dimensions 
than one in 5a?, 5y, &c , and denote the value of Su cor- 
responding to that supposition by dw, we have 

du = d^fU.Sx + dyU.Sy -^ d^u.^x + ... ; 

according to which definition, it appears that du denotes that 
part of the increment of w, as- given by Taylor's Theorem, 
which involves only the first powers of the arbitrary incre- 
ments of the variables on which it depends ; hence, in con- 
formity with this definition, 5<2?, hy^ hss^ &c., must be repre- 
sented by dwy dy^ dz^ &c. ; for if f{x) = a?, then 

fijx}-^ Sw) = .T? 4- o»r; 
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consequently, that part of the increment of f{ai) which in- 
volves the first power oi hx is, in this case, the whole of it, 
^<r, therefore, dw = Sx ; and so on for the others ; 

.'. du = dgU . dw + dyU . dy + dgU . d% + &c., 

where dx, dy, dx^ &c., are the arbitrary increments of the 
independent variables <r, y, z, &c., and are entirely inde- 
pendent of those variables; and du is that part of the 
corresponding increment of the dependent variable w, which 
involves only simple dimensions of dx, dy, d%, &c. ; and 
which approaches nearer and nearer to the value of the whole 
increment of u, the smaller dw, dy, d%, &c. are taken. 

The quantities d^.u.dw, dyU.dy, &c., are called the 
partial diflFerentials of u with respect to w, y, &c. respectively ; 
and du h called the total or complete differential of u, or the 
diflFerential of u, merely. 

75. According to the above definition, the differential 
oi du, or the second differential of u, will be, supposing it 
to involve only two independent variables, 

(Pw = d, {du) dw + dy {du) dy 

= dg{dgU . dx + dyU . dy) dx + dy {d^u , dx + dyU . dy) dy 

= d^u . dai^ + 9,d^dyU . dxdy + d^z^ . dy^, 

because in this process dx and dy are entirely independent 
of X and y ; being, in fact, the arbitrary increments of those 
variables, which might have been denoted by h and k, did 
not a due regard to the precision and symmetry of the 
notation require it otherwise. And in general, if ^^ be a 
function of any number of independent variables x, y, z, t, 
&c., the n^ differential of u will be given by the formula, 
(where the symbols of operation are separated from those of 
quantity,) 

d"«^ = {dx .d:,-\- dy ,dy-\- d% .d^ + &c.)"w. 
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76. Hence, we see that the difiPerential of a function 
has no reference to one variable rather than to another; 
but in its formation, all the variables of which it is a func- 
tion are supposed to undergo simultaneous unconnected 
alterations; and the value of the differential depends upon 
the values of all those alterations. Whereas in forming a 
differential coefficient, one particular independent variable 
only is changed, the rest remaining unaffected ; and a quan- 
tity is produced, wholly independent of the value of the 
alteration which that variable may have received. 

77- It is evident that the rules for deriving differential 
coefficients, suffice for finding the differentials of functions. 

Having, for instance, formed all the partial differential 
coefficients of the first order, if we multiply each by the 
arbitrary increment (or differential) of •the corresponding in- 
dependent variable, and take the sum, we shall obtain the 
first differential of the function; and similarly, the differ- 
entials of the second, and higher orders, may be formed by 
means of the formula at the end of Article 75. 

78. Total integration of a proposed differential, is the 
finding a function whose differential is the quantity pro- 
posed. This operation is denoted by the symbol / without 
any affixed variable; this process, like the former of differ- 
tiation,^ having reference to all, and not to one particular 
variable. 

Ex. If t^ = cr* + y^ 

du s= 2 («rcf a? + ydy) is the differential of u ; 

also f(wdaf-i-ydy)-'--9 since the latter quantity differentiated 

2 

produces the quantity under the symbol of total integration. 

When a quantity is presented for total integration, it 
must be reducible to the form f(v)dv; and then its integral 
'^ fvfi'^)' After this observation, the rules of ordinary in- 
13 
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tegration will enable us to perform the process of total 
integration. 

Ex. To integrate lxL^Z£^ ^ 



y {aidy - ydoi) I at wdy - 



r ^ . y 

/ J dvy putting - = t> ; 



= f Z- = Vl -Hj' 

= \/o^ + y* + C. 



+ C 



79. DiflFerential Equations involving more than two vari- 
ables admit of division into two classes, Total and Partial. 

A total difiPerential equation is one which expresses the 
differential of the dependent variable in terms of the other 
variables and their differentials, and sometimes also of the 
dependent variable itself; it is consequently equivalent to a 
system of equations in which each differential coefficient of 
the dependent variable is given explicitly. Thus, u being 
a function of the independent variables Xy y, %^ the total 
differential equation 

du = Pdx + Q^dy + Rd% 

amounts to the' same as the system of equations, 

d^u^Py dyU^Qiy dgU=^R. 
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Also, z being a function of the independent variables 
w and y, the total differential equation 

Pdw + Qdy + Rdx = 

amounts to the same as the two equations, 

P + Rdj^x « 0, Q + RdyX = 0. 

A partial differential equation, on the contrary, is only 
a relation between all or certain of the partial differential 
coefficients of the dependent variable, and the variables; as 
in the instances, 

(«r - a) dg% + (y - 6) dyX ^z - c^ 

ar^d!j,x + ^wy d^dyX + J^dy« = 0, 

X in both being a function of the independent variables x 
and y. 

Total Differential Equations. 

80. To integrate du = Pdat + Qdy, u being a function 
of the independent variables x and y, and P and Q being 
functions of x and y. 

Since Pdx + Qdy must be identical with dgudx + d^u dy^ 
we have 

dgU^Pj dyU ^ Q9 with the condition d^P-dgQ; 

.-. u^f^P-^ f(y), f(y) involving y only ; 

.•; Q^dyU^dy(f,P) + dyf(y); 

consequently w is known. 

Ex. du = 7 r^ dx - ;^ dy, 

y f '* 

d,u =-7 To 5 d^t^ = - 



Or + y)«' " (^ + yr* 



Digitized by VjOOQ IC 



100 



.-. u ^+/(y) 

ai + y 

w 1 V 

<fl? + y 

81. To integrate dw^ Pdw -\' Qdy ■¥ Rdz^ u being a 
function of the independent variables a?, y, «> 

Since Pdo^ + Qdy •\' Rdz is identical with 

d,w . dw + dy w . dy + d;ft^ . d», 

we have d^u = P, dyU = Q, d^t^ = R ; 

together with the equations of condition Y^\^\v^ t ^^ i^^ 

dyP^d.Q, d^Q^dyR, d,R^d,P, 

which are found by supposing each of the quantities z, w, y 
to be constant in succession (as we are evidently at liberty to 
do, since those variables are independent of one another), and 
then taking the corresponding equation of condition for du 
being the exact differential of a function of two independent 
variables. Hence 

u = fgP + w, w being a function of y and z ; 

also Q^dy{f^P)-\-dyW, 

R^d,if.P)+d,tv, 

which two equations giving the values of the partial differ- 
ential coefficients ottv, its value may be found by the preceding 
Article; and so the value of u completely determined. 



Digitized by LjOOQ IC 



101 

Ex. 1. du = y%dw + wzdy + xydx ; 

dgU = yx^ dyU =s w«^ d^u = wy^ 
.*. u = wyz + w; 
tvz = a?« + dyW, or dyW = 0, 
wy^wy-^- dgW^ or d,w = 0, 
.-. w=C; 
.'. w = wyz + C. 



Ex. 2. df^ = —^ — dw + dy + - — ^ dx ; 

a - » a — «r {a - j^)'* 



dgU = — ^ — 5 d«« =s , dgU =5 



<2?y 



a - « ^ a - «f (a - 5f)* 

Ex, 3. du = a(apdw + ydy) + bdx'^ 
u - - {a^ + ^) + bx + C. 

82. We next come to the consideration of the equation 

Pdw + Qdy + Rdx = 0, 

one of the three variables ^, y, », being a function of the 
other two, which are independent. 

Since the proposed equation may arise from combining the 
results of difiPerentiating two separate equations, we have first 
to examine whether it can be satisfied by a single primitive 
relation between w^ y and x. If w^ y, x be co-ordinates of a 
point, the cases will be distinguished according as the proposed 
equation is the differential equation to a series of surfaces, or 
a series of curves. 
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For example, the equation 

(x - c) ofddff + (z -c) ydy * {^ (<«? - a) + y (y - 6)} dx^ 
arises from combining the results of difiPerentiating the equations 

X -c ^ \x -cl 

for these equations give respectively 

ydgX - wdyX = 0, (a? - a) d^x + (y — b) dyX ^x — c\ 

from which if d^x and d^^r be determined, and substituted in 

dx = d,« . dw 4- dyX . dy, 

we get the proposed equation; which, consequently, cannot 
in general be satisfied by a single relation between ^, y^ x. 
It is the analytical expression of the conditions of the problem, 
to find a surface belonging at the same time to conical surfaces^ 
and surfaces of revolution about the axis of x* 

83. To find the equation of condition for 

Pdw + Qdy + Rdx = 0, 

admitting a solution of the form /(a?, y, x) = C. 

If the proposed equation can be satisfied by a single 
relation between a?, y, and ar, 

P Q 

is thie difiPerential of a function of two independent variables ; 

P Q 

with the condition dy |~j = ^* f-^) 5 

or, since P, Q, R may contain x^ which is a function of jp 
and y, 

^iy) Q + ''o (I) dy^ - dw (I) + di^y (I) *«-* ; 
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therefore, substituting for d,z, dyX, their values, 

M,,(|)-«'^.,©-J"'«(|)-'"'..,(|) = 

and performing the difiPerentiations, and reducing, we get 

the equation of condition for 

Pdw + Qdy + Rdx = ; 
admitting a solution of the form /(a?, y, z) = C 

84. If the above equation of condition is not satisfied, 
then the equation 

Pdw + Qdy + Rd% = 

cannot, by being multiplied by any factor, become susceptible 
of a solution of the form 

/(a?, y, z, C) = 0. 

For suppose F to be a factor which renders 

Pdw + Qdy + Rd% 

the immediate difiPerential of some function u of w^ y, z, con- 
sidered as independent; 

d(,,(FJ2) = d^.,(FP); 
or rd^,,P + Pd(„F=rd^,,Q + Qd(,/ 

Frf(,,Q + Qd^,,r = Vd^^yR + Rd^^^V . (1). 

rd(.,i? + fid(„r = rd> + Pd(./ 

Hence, multiplying the first of these equations by if, the 
second by P, and the third by Q, and taking their sum, the 
factor V disappears, and we find 

P(d(,,i?-d^,,Q) + Q(d^„P-d(„i?) + fi(d(„Q-d(„i*) = 0, (2) 

the same equation of condition as in the preceding Article. 
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85. If this equation be satisfied, which will be the case 
only when the proposed admits a primitive of the form 

f{w, y^ Xj C)^ 0, 

equations (1) afford a means of determining F. Thien 

du = VPdw + VQdy + VRdz, 

whence u can be found by the method of Art. 81, and 
t^ + C B is the required relation between Wj y, and %. 

Or, without determining F, we may integrate considering 
one of the variables constant, and add an arbitrary function of 
that variable ; then differentiate the result with respect to that 
variable and compare it with the proposed equation, and so 
the correction will become known. 

Obs. In the majority of cases which present themselves, 
the factor F is capable of being determined by inspection. 

Ex. (ay - hss) d,v + (ex - aw) dy + {bx - cy) dz « 0. 

Divide by (ay - b%) (bw - cy), and the resulting equation 

dx (cz — aai)dy dz 

bw - cy {ay - bz) (bx — cy) ay --bz 

satisfies the conditions, considering cT, y, z^ as independent, 

dyP^d.Q, d^Q^dyRy d,R^d,P; 

and therefore the first member may be regarded as the dif- 
ferential of some function, u, of <v, y, z considered as in- 
dependent ; 



d^u • 



bw — cy* 
\ 



w = - log {bw - cy) + w ; 
o 
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1 



ay — b% 
cz — aa c 1 



I 



(ay - bz) {hw — cy) b ba — cy 
or a«w = — T 



^ b ay -^bx 

1 adv — 6dir 

•. dt<^ - — . — ^ ; 

6 ay^bx 

.-. w; = - - log (ay - bx) + C; 

1 - fbof - ci/\ * ^ bw — cy , 



1 , fbw -cyS ' 



Ay - feW oy - bx 

Obs. When in an equation of this sort, the di£Perentials 
enter above the first degree, it is not integrable unless it can be 
resolved into rational factors of the form Pdw + Qdy + Rdx ; 
for whatever be the integral, it must upon di£Perentiation pro-' 
duce a result of that form. 

86. If the equation 

Pdw + Qdy + Rdx = 
/ 
be susceptible of a primitive of the form /(d?, y, x, C) « 0, 
and be homogeneous and of n dimensions with respect to 
^9 9) ^9 then, putting nf^vx^ y^wx^ and dividing by «", 
it becomes 

S(vdx + xdv) + T{wdx + xdfv) + Udx = 0, 
#9, r, U being functions of v and w ; 

dx Sdv + Tdw 



or - 



X Sv^Tfv-^U" 



hence the second member is an exact differential since the 
first is so, and it may be generally integrated by inspection, 
or by the method of Art. 80. 
14 
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Ex. 1. (y + z)d3ff+{iv + z) dy + (a? + y) d» = ; 
put w^vz, y^wzy 
dz (w? + 1) dv + (v + 1) dw 

dz . d (t> + w + vw) 
or — + i — ^^ ^ = ; 

z * V + ti? + vw 

.\ log «* (© + w + vw) = log C, or a?» + yi^f + a?y = C 

Ex. 2. (y* + y« + «*) da? + (af^ + a?» + «*) dy 

+ (^ + a?y + y^) d» = 0. 
Make w a v;^, y as fi^jir, 

d« («?• + «r + 1) dt? + («>* + 1? + 1) dw 

" «f V (w^ + w + l) + w{v^'^v + l)+v^ + vw + w^ 

but of the latter fraction, the denominator 

^vw(w + v -^1) +v(w -{-1 -^-v) +w(v -^-1 +w) 
= (vw + V -^ w) (v + w + l)y 

and the numerator 

= (l + u + w) d (« + w + vw) - (v -^ w + vw) d (v -^ w) ; 
dz d(v + w + vw) d(l +v + w) 



— + 

z 



v-^w-^vw 



= 0. 



1 +v + tl? 

^\ 1+v+w j ^ 
.-. «ar+i8?y + a?y= C(«r + a? + y). 
Ex. 3. (a?* - y* + «*) da? - «*dy 
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putting of^vzy y^wxy this is reduced to 
dw « (1 + tj* - w^) dv ; 

: jtg-tH» ^ c^ (Ex. 3. Art. 18.) 



e-^ 



w ~v 



or = Le f» + C«, 

Jir being constant under the sign of integration. 



Total Diffisrential Equations that do not admit of a Single PrimitiTe. 

87. We have seen that .the equation 
Pdof + Qdy t R^^ = 0, 

when the equation of condition (Art. 83) is not satisfied, does 
not admit of being derived from a single primitive equation 
involving two independent variables. The integral in this 
case will be exhibited by a system of two equations; and 
the proposed equation cannot be regarded as the di£Ferential 
equation to a surface, but to a system of curves in space, all 
endowed with some common property. 

Ex. 1. d» s aydof + bdy. 

Since the equation of condition in this case is not satisfied, 
w and y cannot be independent, and we may assume y-f{(xi)% 

••. d%^af(af)daf'{-bf'(w)daf; 

which, with y-f{(xi)y constitutes the integral of the proposed 
equation. 

In general, if F be a factor which makes Pdw + Qdy an 
exact difierential, considering z as constant, and we find 

JiVPdcB ^VQdy) ^w -\- <t>{z)x 
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it is evident that 

w + («) « 0, together with d^w + (j/ (x) - Fi? = 0, 

satisfy the proposed equation, where <f> (%) denotes any func- 
tion of z 

Ex. 2. zdw + ^dy + ydz » 0, 

y + « log« + {%) = 0, log^ + d,d>(%) « -. 

Ex. S. {^ (^ - o) + y (y - ft)} d%^(x - c) (wdw + ydy)^ 
0^ + y* + 20(«) « 0, /p (j? - a) + y (y - 6) + (« - c) 0'(«) « 0. 

FWial Differential Equations. 

88. In partial differential equations of two independent 
variables, the differential coefficients of the first order d^z^ 
dyZ, of the dependent variable z^ are usually denoted by the 
symbols p and q\ and d^z^ d^d^z, d^z, the differential co- 
efficients of the second order, by r, s, t^ respectively. A 
partial differential equation is said to be of the n^^ order, when 
it involves one or more of the partial differential coefficients of 
the dependent variable of the n^ order ; but none of a superior 
order. To be the general equation of the n?'^ order, it ought 
to contain the independent variables, and the dependent va- 
riable together with all its partial differential coefficients from 
the first order to the n^^ order inclusive. To integrate a 
partial differential equation, is to find for the dependent 
variable, an expression between the differential coefficients of 
which, that relation exists which is indicated by the proposed 
equation ; and under the most general form possible. 

Equations of the First Order. 

89. The complete integral of /(zp, y, «, j», q) ta o,. the 
general equation of the first order, will involve one arbitrary 
or general function. 
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For let « = /'{a?, y, «f, ^C^?)} «0, be an equation by 
virtue of which j^r is a function of the independent variables 
w and y, v being a known function of w^ y, and x. Then 
p and 9 are given by the equations d,u » o, d^u s 0, each 
of which will involve d«0(f) or ((/(v), and may involve 
(f>(v); consequently, between the three equations u^O, 
dgU 8s 0, dyU B 0, it will be possible to eliminate (v) and 
0'(«), and there will result a relation /(.i?, y, «, p, g) = 0, 
wholly independent of the form of the function (v) ; and it 
IS evident that in general more than one function cannot thus 
be eliminated. Conversely, an equation of the form 

f(^y y> «> Pj q)^o 

being proposed, its complete integral to have all the gene- 
rality possible, must be of the form F {w^ y, at, 0(v)} "^O, 
where the form of <l>(v) is perfectly arbitrary. 

For example, let 

w = » + mof + «y - (t>) = 0, 

where © = (^ - a)* + (y - 6)* + (z — c)* ; 

then dgU « p + 1» - 0'(tj) {2 (d7 - o) + 2 (» - c)p] = 0, 

dyU ^q + n - 0'(») {2 (y - 6) + 2 (« - c) q} = 0. 

Hence, transposing, and dividing one result by the other 
to eliminate (p'ip)^ we find 

p + m (/r -- q) + (y - c) j? 
. ? + » "" (y - 6) -»- (« - c) ?* 

or p {y - 6 - « (» - c)} - g' {a7 - a - w (ar - c)} 

e n (^ - a) - w (y - 6), 

the partial differential equation of which the complete inte- 
gral is 

% + ma: + ny = \{*v -€i)* + (y - by + (« - c)*}. 
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90. To integrate an equation in which only one of the 
differential coefficients of the first order enters with w^ y^ 
and z. 

Let the equation be 

Integrate it, considering y as a constant, and in place of 
the arbitrary constant C, add a function of y of arbitrary 
form. The resulting solution, containing one arbitrary func- 
tion, will have all the generality that can be attained. * 

Ex. d.«(^' + y*)«i^ + »»; 

d^z 1 



»« + y» I^Jrit 



= 0; 



.'. tan"* — tan-* - = tan"* (y), 

y y 

or z^af^(j^ + az)(f>(y)y 
(f>(y) being arbitrary in form. 

The equation / (a?, y, «, 9) = is similarly integrated, 
the correction in this case being an arbitrary function of «r. 

Ex. DoydyZ + nar « 0, 

91. To integrate the linear equation of the first order, 

P, Q, and R being functions of w, y, and z. 

Let the primitive he F(af, y, ar) »= 0; therefore, denoting 
d^^^F(wy y, z) by F\ai)y and so on for the other co- 
efficients, we get 

F'{ai)^F'{z).p^Oj 

F'(y) + F'(z).q«0; 
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.-. PF{ai) + QF{y) + RP{«) * 0. 
But dF{afy y, «) = /^(zr)da? + JP'(y)dy + F'(af)d«-0; 
.-. PFiai) dof + PFiy) dy + P/^(«) d» = 0, 
and PF(ai) dw + QF'(y) do? + RF(x) dw^O; 

which is satisfied by 

Pdy - Qdw A 0, 



-0,| 

r (2). 



Suppose that by integrating these equations, either se- 
parately or conjointly, we obtain M ^ a, JV «= 6, two relations 
between the three variables and the arbitrary constants a and 
bf which satisfy them. By means of these two equations, 
any two of the variables as y and « can be expressed in terms 
of a and 6 and the third variable of. The complete primitive 
becomes 

F(a^ y, «) « (^j «> 6) = 0, 

^ and the differential of (^, a, 6) must by virtue of equations 
(2) be identically equal to zero, therefore (p (^, a, 6) cannot 
contain J7, and 

O«0(a, 6) = 0(Jlf, JNT), 

or M^f(N). 

The assumptions (2) satisfy equation (1) independently 
of the forms of F{fj)^ F(%)j that is, independently of the 
form of F\ therefore the form of 0, and consequently of/, 
is arbitrary. Hence M ^f{N) being a solution of the pro- 
posed equation and containing one arbitrary function, is the 
general primitive. ^ 
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92. For the success of the method, it is in general 
necessary that of the three equations 

Pdy - Qdw = 0, Pd« -• Rdm = 0, Qd» - Rdy = 0, 

one at least should contain those two variables only whose 
increments or differentials appear in it. They are called the 
reducing equations of Pp -¥ Qq^ R\ and may be easily re- 
membered under the form 

dw dy d« 

93. By a similar process, partial differential equations of 
three or more independent variables can sometimes be inte- 
grated. If jv be a function of a^ y^ t^ and 

r» + Pp + Q? - Ry 

where n » di», we have 

R-Pp- Qq 

ns ^ i; 

T 

which, substituted in 

dz m ndt + pdof + qdy^ 

I gives Tdz - Rdt = p {Tdx - Pdt) + q (Tdy - Qdt). 

dx dy dz dt 
And if the equations — = — = —=—, or 

Tdz - Rdt « 0, Tdw - Pdt « 0, Tdy - Qdt = 0, 
give L^a^ M ^h^ N ^Cj 
then the solution is 

(i, M, N) = 0, or L «/(ilf, N). 

94. To explain the geometrical meaning of the solution 
of a partial differential equation of the first order. 
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If we regard a, y^ x as the co-ordinates of a point, and 
the proposed differential equation as the equation to a system 
of surfaces, M^a, Nmb, are the equations to two surfaces, 
which, being satisfied by values of w, y, x which satisfy the 
differential equation, conjointly represent a line on a surface 
represented by the primitive. By giving any alterations to 
a and b we obtain other lines in space; but in order that 
these lines may always lie on the surface in question, these 
alterations must not be independent but connected, according 
to some law expressed by the relation 

a-/(6), or M^f(N). 

For example, the equation 

pw + qy ^x 

is the differential equation to a conical surface, and 
- = a, - = &j 

X X 

are the equations to any straight line through the vertex. 
If a and b undergo all possible variations consistently with 
the restriction imposed by the nature of the directrix, the 
assemblage of straight lines thus formed is the conical sur- 
face. 

We shall now apply the preceding method to examples. 

Ex. 1. pof + qx-k-p^Of 
dx » pdw + qdy^ 
wdx s pwdw + q^dy 
- yda B padw 4- qxdco ; 

.-. wdx + ydaf = q (wdy - xdw) ; 

.% wdx + pdw « 0, 

wdy - xdw s* 0, 

xdx + ydy =» 0. 
16 
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From the last, «■ + y* « «*; 



d« dw 

+ — = 0; 



.-. sin"'^- + log^= logft; 
a 

Ex. 2. p (a? - o) + gr (y - 6) B « - c. 
This compared with Pp -h Qg = iJ, gives in the place of 
Pdx - JBd^ = 0, Qdx-- Rdy « 0, 
the reducing equations 
(a? - a) d« - (« - c) do? «= 0, (y - 6) dar - (« - c) dy « ; 



" )»-€ "^ \«-C/ 



Ex. S. fuji + nfa 1; 

.". mdiir - dd? e 0, nd«-dy«0; 

.\ y-nx^f(af-mz). 

Ex. 4. (a?-w»)p + (y-n«)g«0; 

.-. (w - f»«) d« Bs ; .% « = a, 

(a - m«) dy-^iy- nx) dw «= 0, 
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or (af-wa)dy-(y-na)da?«0; 



w — ma 

1^::^= 6; but 0-/(6), 



-, /a7 - mz\ 
••• ««/ ). 

7. p^+gy».,^, l-i+/fi_I). 
« a '' \y xl 

10. (aor + 6y) p + (a'^ + Vy) q « c«. Assume Xme\ 

95. The following examples require artifices in the com- 
bination of the reducing equations. 

Ex. 1. {y-6-n(»-c)}ji-{dr-a-m(»-c)} J 
«n(^-a)-i»(y-6), ^ 
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the equation to surfaces of revolution. Here the reducing 
equations are 

{« (07 - o) - m (y - 6) I doz - {y - 6 - n (» - c)} d» = 0, (l). 

{y — b - n (z ~ c)] dy -^ \w - a - m(« ~ c)} dw ^ 0, (2). 

I^p - a - m (« - c)} d» + |« (4? — o) - f»(y - b)]dy = 0, (3). 

Multiply (1) by 07-0, and (3) by y-6, and add^ 

••. |« (^ - a) - 1» (y - 6)} X 

{(w " a) dw -k- (y - b) dy-^iz-^ c) d«] « ; 

••. (w - a) dd7 + (y -^ 6) dy -i- (x ^ c) dz = 0; 

.-. (w - ay + (y - 6)* + (« - c)* = a. 

Multiply (1) by m, and (3) by n, and add, 

.\ {n (07 - a) - w (y - 6)} (md^ + «dy + dx) = ; 

.-. wd^ + ndy + d« = 0, 

w^ + «y + ar = fi^f(a) ; 

.-. « + m^ + ny «/f («? - «)* + (y - 6)' ^ («f - o)»}. 

Ex. 2. (^ + y + «) df« + (^ + y + «) d,» 

+ (^ + 07 + ») dyX ^t-^-w + y. 
The reducing equations are 

(07 + y + «) d« = (/ + 07 + y) d/, 
(o7 + y + «) do7 « (/ + y + «) d/, 
(07 + y + Jir) dy = (* -f « + «) d^; 
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••• (* + » + ») (dx - dt) » (^ • «f) dti 

(w + y ■{• «) {dw - dt) = (t - w) dtf 

(a? + y + «) (dy - dt) ^ (t -y) dt. 
Also, 
(or + y + «) (do? + dy + d» + d^) = 3 (^ + «? + y + «) dt; 

d% - dt dt dw -^ dy -k- dx + dt 

^-« « + y + « 3(a? + y + » + 0' 

» 

.«. (a? + y + « + 0'(« - 0' * «> 
Similarly, (/p + y + » + t).(w - ^)'* = /3, 
(j? + y + « + 0-(» - 0' •* 7' 
••. (ar + y + ar + Q. (»-/)» 
^F{(w + y-\'X + t).(af -t)\ (a? + y + « + t).(y - ^)«}. 

96. To integrate F(w, y, )ir, p, q) « 0, when it contains 
terms of more than one dimension in p and 9, 

In order that 

dx^pdw -{■ qdy 

may be a perfect differential, we must have 

^ip)P + 9d^z)P = ^i»)9 + P^c,)?- (2). 

If from the proposed equation we determine 

equation (2) becomes an equation for the determination of p 
of the form 

d^^^ + Af d^y)P + iV'd^.jp = Z, 
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the integration of wBich depends on the integration of one 
of the equations 

dp-Ldw^Oy dy-'Mdw^Oy d«-iVd^ = 0. 

Let p=fi^^ y^ «i a) be found from these equations, a 
being an arbitrary constant. This value of p and the corres- 
ponding value of q found from the proposed equation, being 
substituted in d«=pdw + qdy render it an exact differential; 
and thus a value of x will be obtained involving two arbi* 
trary constants a and 6; and this will consequently be the 
complete primitive. The general primitive may be obtained 
by putting 6 = >|r (o), differentiating the equation with regard 
to ttf and eliminating a. The result containing one arbitrary 
function is as general as any solution which the equation 
admits. 

Ex, 1. p» + 9*= 1, 
q^\/l-p^ d^q = " !^ d^p, d,q^ ^ —f==:rd^p\ 

/ pdmP P^dgp 





-/" 


- da? = 0, 


p 


./' 


dm 


v^ 


Vi-p' 


. p = 


a, 9: 


= y/l^a\ 


dz^ 


adw-^- 


'\/\-a*dy. 



0; 



» = aa? + \/l-«*y + 6, 
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or z^ aof + y/\ -a*y + («)> 

If a be eKminated between the two latter equations, the 
general primitive results. 

This is a simple case of the problem of finding surfaces 
of equivalent area to a given surface, that is, so that any 
cylindrical surface parallel to axis of x may always intercept 
equal areas in the required and given surfaces. 

If P, Q, be the values of d,ar, dyX in the given surface, 
thie general condition is 

which, if the given surface be a plane, leads to the equation 
of this example. . 

Ex. 2. ar = p?. 



% 



y + o 

9S 



'^ + 0(a). 



r»-d.0(o). 



Equations of the Seeond and Higher Orders* 

97* Here, besides the coefficients p and q of the first 
order which may be involved together with Wy y, and Zy the 
equation must contain one or more of the coefficients of the 
second order r, «, t\ so that in its most general form it 
will be 

FifOy y, «r, py qy r, «, f) - 0. 

In partial differential equations of the second order, we 
cannot be certain of the form of the solution, nor pronounce 
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beforehand how many arbitrary functions it ought to contain. 
For let 

be an equation containing two arbitrary functions of v and 
w two known functions of w, y^ and « ; then p and q will be 
given by the equations 

dgU = 0, dyU = ; 

and r, «, and ^ by the equations 

c^ii «B 0, dt,dyU B 0, d^tf s 0. 

These together with te » make six equations, into which the 
six quantities 

(p{v\ >/.(ir), f(i,), >/.'(«;), 0"(f^), >^"(tD), 

may enter. Consequently, it will not generally be possible to 
eliminate these six quantities and obtain a relation between 
^9 Vj ^9 P9 99 ^9 ^9 ^9 iudcpeudent of the forms of the func- 
tions <f} and yj/ ; although particular cases do occur in which 
this can be effected. 

In general, it u = contain n independent variables, and 
m functions of the form (v), where v is a determinate func- 
tion of the variables; these functions can certainly be eli- 
minated between the equations obtained after (r) differentia- 
tions, if 

^ _ n(n + l) ^ n(n-H)... (n + r-1) ^ _^^ ^ ^^ 
— 1T2 — ••• \r >«»(r+l), 

the former expressing the number of equations, and the latter 
the number of functions. 

9& As an example of formings partial differential equa- 
tion of the second order, by the elimination of two arbitirary 
functions, we may take 

w = y - ^0 («) - ^ (ijf) = ; 
dyU m 1 - w<f/{z) q - yf/{x) 9-0; 
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.^ r + «0 (x) + pq^'iss) = 0, 

.\ ^r - 1)« + (g« - pt) («) e 0, 
or (fr - 2p^« + ^t = 0. 

99. The equations 

d^«f + Pdy« + Q = 0, 

where P and Q are functions of w, y^ and Xy must be in- 
tegrated as equations between two variables, y being regarded 
as constant in the former, and w in the latter ; and arbitrary 
functions of those variables respectively, being introduced 
instead of constants. 

100. The equations 

d,dyX + Pdt,x = Q, 
d^dyX + PdyX = Q, 

where P and Q do not contain «, are reducible to the case 
of Art. 90, by considering dgX or dyX respectively as a single 
quantity v. 



Ex. dgdyX + j^ dgX = a^. 






This being a linear equation in d^x which is made in- 
tegrable by the factor —r , we have 



:;^ - °X;^-«S''^^+*<-'^»^'^'W' 
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ax 



101. Similarly, equations of the forms 

f(!»9 Vf «, d.», djj8f, &c., <«) = 0, 

/(a?, y, 5f, dyj8f, dJsB, &c., d;») = 0, 

may be treated as if they contained only two variables; 
arbitrary functions of y instead of constants being introduced 
into the solution of the first, and arbitrary functions of <r 
into the solution of the second. To this case may also be 
reduced the equation 

/(^, y. ^», d^d^x, dld^x, &c., d;d;x)^0; 

for by putting <^jvc=v, it becomes 

/(w, y, V, d,t>, d>, &c., dJ^v) « 0, 

which will give a value of v containing m arbitrary functions 
of y ; and then d^% = v will give x involving n arbitrary 
functions of <r. 

102. To integrate the linear equation of the second order, 

Rr + Ss + Tt=r, 
where JR, S^ T, V are functions of cT, y, x^ p, q. 

By means of the relations 

dp =! rd^ -f «dy, 
dg «= ^dcT + tdy ; 

eliminating two of the three coefficients, r and /, from the 
proposed, we get 

Rdpdy + Tdgda? - Vdwdy ^8{R (dyf ^T{daif - Sdwdy}, 
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which is satisfied by 

Rdpdy + Tdqdai - Vdwdy = 
R (dyy + T {dwy - Sdiedy = 



(0- 



Let M ^ a^ N - b, be two relations between w, y, x, p, q^ 
and the arbitrary constants a^ b^ which satisfy these equations ; 
then M =(f) (JV) satisfies the proposed equation. This will 
be shewn by proving that it can reproduce the proposed 
equation. 

Let dy = mdw ; .-. Rm^ - Sm + T « 0. 

For each roof of this equation, we have 

dy - mdw = 0, Rmdp + Tdq — Vmdx «= ; 

.'. dy — tndWf \ 

n Rm 

(2). 

d% = pdoB 4- qdy. 

Hence M — a gives on differentiation 

= dgM . dai + dyM . mdw + d^M {pdw + qindx) 

+ d^M.dp-\'dqMA-—dai- -7=rdp\ , 

wherein all the known relations (2) having been introduced, 
dw and dp must be independent, 



Vm ^ Rm ^ 
«^="^ do?- — dpy 



Vm 
T 



= d^M + mdyM + d^M . (p + mq) + -;^ (^^Jf, 



0^dpM--~dgM; 
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Vtm. 

.-. d,M =» - {tndjfM + (p + mq) d,M + -=;■ d,ilf}, 

d,M^^d,M. 

Vtn 
So d,N = - \md^N + (p + mq) d.N + -=- d^N}, 

d,N=^d,N. 

By differentiating the assumed equation M =^ <(> (N) we 
have 

Now 

Vm 
dM « - {mdyM + (p + mq)d^M + -=- d^M} dw + dyMdy 



+ d^M{pda + g'cfy) + -— dqM.dp + d^-Sf • dg 

« {dyM-vqd^M) (dy-mdw) + -^ {Rmdp -l- Tdg- Fmrf^), 

and a similar value exists for dN', 

.•. (dyM-k-qd^M) (dy- mdw) +-~{Rfndp + Tdq-Vmdw) 

« 0'(JV) {(dyJVT + gd^JVT) (dy - mda^) 

■ +^(Rmdp+Tdq^ 

which may be put under the form 

Rmdp + ydg - Vrndw = w (dy - mdw) 
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or Rm {rdx + ady) + T (sdw + tdy) - Vmdx=ia{dy --mdai), 
where dw and dy are independent ; 

Rms + Tt^w; 

m 

Hence, M = (f> {N) satisfies the proposed equation. 

According as the roots of 

Rm^-^Sm^ T = 

are unequal or equal, we are thus supplied with a total or 
partial differential equation for the determination of %. 

Obs. As the reducing equations (l) may contain 

^, Vi ^9 P, ?5 

and as these together with dz « pdw + qdy will generally 
lead to an equation containing three variables, which will 
not always admit of a single primitive (Art. 83.), it may 
happen that the first integral of the proposed equation 
cannot be determined; but we must not thence conclude 
that the proposed equation does not admit of being solved. 

103. Hence, to integrate the linear equation of the 
second order 

Rr-\-S8+Tt=V, 

the process is to obtain a value of m from the equation 

Rm^-Sm+ r = 0, 

to substitute it in the system 

dy-mda?=0, Rmdp ■{- Tdq^^ Vmdw^ (l.) 
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to satisfy these, conjointly or separately, by two relations 
between a?, y, «, p^ g, 

then to put M =: (f> (JV), and to integrate this equation of 
the first order. 

104. If JK, Sj T'9 be constant, and V a function of /v 
and y only, then the values of m will be numerical, m and 
n suppose; and the integrals of equations (l) will be 

y - maff = a^ Rmp + Tq^ mf^V+by 

where, previous to integration, mw + a is substituted for 
y in F, and after integration the value of o, viz. y — ma^y is 
restored: consequently calling this value Fi, since Rmn = T, 
we have 

Rp + Rnq = Fi + (p'(y - mw). 

Next to integrate this equation of the first order, we have 
the reducing equations 

dy -ndx^tOy 

Rdz - {Fi + (j)(y-'mw)} dti7 = 0; 

.•. y -nw ^ ay 

Rx^f.Vi'^f»(t>'(y-nia;)+fiy 

na '\- a being substituted for y before the integration is per- 
formed, and afterwards the value of a, y - ncoy restored ; 
this gives ^ Fi = Fg, suppose, 

and /,0'(y-mj?)=/,0'{(/i-w)^ + a} =— — .0(y-m^); 

hence, including the constant multiplier under the sign of the 
function, 

R% = J^2 + (y - fnsB) + >^ (y - nx). 

Ex.1. ^^% = €?^ySiSy or r-a^/eO, 
or rdxdy - a^tdydat = 0, 
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or dpdy - s(dyf - a^ \dqdx - «(da?)^} = ; 
.-. idyY^a'(dw)\ 
dpdy = a^dqdw^ 
dy = :k ados, 

\^. dy =: adx^ y - ax =^ a^ 

dp^adq, p - aq^ fi^f{y ^ ax), (i). 

2****. dy=^- adx^ y + ax ^a^ 

dp = - adq, p + aq^ji^ (^'{y + ax). (2). 
By adding and subtracting (l) and (2), 
^P =/'(y - «<») + 0'(y + ax), 
Qaq^ip^iy-^-ax) -f(y-ax); 

.•. 2ad« = 2apdx + 2aqdy = 0'(y + ax) (dy + ad<r) 

""/'(y - ax) (dy - ad^) ; 
.•• % = <p(y + ax) +f(y - ax). 

Ex. 2. r--^«+^ ^ = 0, 

q (f 

2© ®* 

rdwdy adxdy + ^ tdydx = 0, 

2© ©* 

dpdy - « (dy)^ adxdy + ^ {dgd<r ^ s (dxf] = o ; 

.-. (dy)' + 2?dyda?+^(de»)8 = 0, 
dj^dj^ + — dqdx = ; 
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P 
.: dy '^ — dw, 

q 

p 
9 



qdp-pdq p 

••• P-9/W = 0. 
To integrate this equation of the first order, we have 
dijr s= 0, or z — a^ 
dy +/W dw^Oj or y + /p/(o) = /3 ; 

this is the equation to the surface generated by a straight line, 
subjected to pass through three given fixed curves. 

Ex. 3. r - a^t «= a?y, 

*» = ^ flT^y + /(y + a^) + (» - a^). 

7. r+ (a + 6)« + «6^ = ^y> 

« - g^ ar*y + — (a + 6)^* = (y - aa?) + >/f (y - 6a?). 

8. (1 +pg + g*)r + « (q^-p^) -(i+pq+p^)t^ o. 
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A TREATISE 



CALCULUS OF FINITE DIFFERENCES. 



SECTION I. 

DIRECT METHOD OF DIFFEB£N€£8. 

Definitions and Principles. 

Aet. 1. In the Differential Calculus it is the first term 
only of the series, arranged according to ascending powers 
^'of A, expressing u^^^^-u^ (where u^ is any function of w 
and u^^^ the same function of a? + A), or rather the coefficient 
of h in that term, with which we are principally concerned, 
and which we usually write hd^u^. But in the Calculus 
of Finite Differences, it is the whole of that series which 
forms the object of our investigations, and it is usually written 
At^,, so that 

2. It is common, however, to suppose the finite incre- 
ment which the principal variable x receives, to be 1 instead 
of A, both for the sake of simplicity, and because that is the 
value of the increment when w, is regarded as the general 
term of a series; and in that light it is by far the most 
frequently regarded in finite differences; so that 

Aw, »w,+i-w,. 
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There are, however, in this subject, several important 
theorems which it is advantageous to investigate on the hypo- 
thesis of an indeterminate increment A, for the principal vari- 
able, instead of unity, as the process is the same on either 
supposition, and the result one of greater generality. And 
if in other cases it should be desirable to introduce the same 
hypothesis, the expression must be prepared by first writing 
hx instead of <r, then performing the operation on the usual 

supposition of A«r = 1, and in the result writing - for %. 

h 

3. By a Series is meant a regular progression of terms 
increasing or decreasing in magnitude according to a certain 
law; hence, when that law is given, and also the place of 
any term in the series, the magnitude of the term may be 
found, and thus the successive terms of the series may be 
produced in order. The place of any term, in a series is 
assigned by giving the number of terms by which it is re- 
moved from some one which is considered as fixed. This 
number is called the index of the term to which it belongs. 
Thus in the series 

0, 1, 8, 27, 64, ... aP^... 

taking the first term as the point of departure, we have the 
corresponding series of indices 

0, 1, 2, 3, 4, ... a?, ... 

If the series be continued backwards, the indices must 
be considered as negative; thus the backward continuation 
of the above series gives the terms 

• •• — oi/y ... — 27, ""8, "~ 1» 
with the corresponding indices 

— ,r, ... —3, —2, — 1. 
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4. Since the magnitude of every term is determined 
solely by its index and by the law of the series, it follows 
that any term is a certain function of its index, the form 
of which does not alter in passing from one term to another, 
but remains the same throughout the whole series. Thus 
in the above series every term is the cube of its index. 

This function analytically expressed is called the general 
term of the series; (in the above series the general term 
is w^;) and it is evident that all the terms of the series 
will be produced from it in order, by substituting successively 
for the index a?, the progression of natural numbers 

...-2, - 1, 0, 1, 2,... 

The general term of a series is usually denoted by w,, 
where u, is a certain function of w determined by the nature 
of the series. Thus, u„ denoting the general term, the series 
will be 

...«_25 «*-lJ ^09 «*1> ^29 . ..«#-! J ^»9 ^M+19 .•• 

5. The excess of any term w,+i above that which immedi- 
ately precedes it, or the function u„+i - u, is called, as has 
been stated, the Difference of the function i*„ and is denoted 
by £^u^y Xf^ certain cases, which will however be expressly 
mentioned, we shall take Au^ to mean ^s+hT^x*). 

It is obvious that u^^+i - u, is itself in general a certain 
function of s, the nature of which is entirely dependent on 
that of the original function «,, from which it is derived, 
and susceptible of a difference. 

The difference, consequently, of the function Aw, (which 
must be considered as having Au for its characteristic, in 
the same manner as w, has u) is 

A (Aw,) « Aw,^-l - Aw,, 

which is usually written A^w,. 
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In like manner 

A (A"w,) = A^w, = A*M,+i - A*t^,, 
A (A^w,) = A*w, = A^w,+i - A'w,, 



A*i*, = A""*w,+i - A* * w,. 

6. Hence if in any function of a we change a? into a? + I9 
and from the result subtract the proposed function, we obtain 
the first difference of the proposed function ; and the second, 
third, &c. differences are formed, each from the preceding, 
by a similar operation. To determine these differences of 
given functions, and to investigate the relations which hold 
between differences of any orders and the functions from 
which they are derived, is the object of the direct method 
of differences. 

We shall now proceed to give instances of finding the 
differences of various functions, according to the above defi- 
nition. 

Differences of Explicit Functions. 

7. To find the difference oi au^-if c, u^ being any func- 
tion of 07, and a and c quantities independent of x, 

A (aw, + c) = aw^+i + c - (aw, + e) = a (w,^.i - t*,) = a A«,. 
Hence, making a = 0, Ac = 0. 

8. To find the difference of the sum of any number 
of functions of ao, 

A (w, + «, + «?,) = w,^.i + v,+i + ti?,^.i - {u, + t>, + «r,> 

= W,+ , - W, + tj,+i - «, + ti;,^j - w^ 

= Aw, + Au, + Aw,. 
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9. To find the difference of the product of two functions* 

= w,At>, + u,At#, + Aw,. Atj, =» w,Atj, + i?,+iAw,. 

10. To find the diflerence of the quotient of two func- 
tions. 

^A «*+! «** («** + Aw,) «, - (v, + A«,) w. 






^J^+ltJ* 



v,Aw, - w,Av, 

11. To find the diflTerence of the continued product of 
any number of successive values of a function. 

A (w,w,4i ... w,+„) = w,4.iW,+2 ••• «#+ii+i - «**^*+i ••• «*,+« 

Hence in the particular case where w, = a + 6a? 
A (w,w,+i ... w,+„) = w,+iW,+s ... w,+» . (w + 1) 6. 

12. To find the difiference of a fraction whose numerator 
and denominator are the continued products of any number 
of successive values of two functions w, and t?, respectively. 

/W,W,+ i...W,+»\ W,+ i«^2-^«Vh^ r V 

\^*^*+l ••• ''^*+m/ tJ^V^+i ... «,+„ + ! 

Hence A = j 

»*^*+l ••• ^*+« ^*V*+l ... ^'jf+m + l 

and in the particular case where «, « a + 6^, 

1 (w + 1) 6 
A = . 
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13. To find the differences of any rational integral func- 
tion, and to shew that the n^^ difference of a rational integral 
function of the n^ degree, is constant. 

Let Ug^Aw^-^ Baf^'^ + . . . + Kof -»- Z, be a rational integral 
function; then its first difference is 

Aw, = J{(a?+l)»-a?»}+fi{(^+l)"-^-^-^} + ... + /(2d?+l)+jr 

which is a rational integral function, one degree lower than 
the original function. In like manner, for the difference of 
this, or the second difference of u^, we have 

AX = n (71 - I) Jcr»-» + jB^a?—^ + ... + A, 

and so on; and for the n^ difference, we have 

A"w, = « (» - 1) ... 3.2. l.J. 

Hence the n}^ difference is constant, and the differences ' 
of all orders superior to the n^ vanish. 

Also A" Qxf) = 1 . 2 . 3 ... 71. 

14. To find the differences of a*. 

Aa' « a»+^ - o* = o' (a - 1), 
A«a' = (a - 1) Aa' = a' (a - 1)% 
A»a' = a' (a - 1)". 
Also A a''' =» a*'+^«'' - a'^" = a''* (a^*'* - l). 

15. To find the difference of log ©,. 

A (log*'*) = log<?* + i - logi), = log ^^ « log f 1 + — ^y 
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16. To find the differences 6t isin©, and cost?,. 
Asint>,= sin(v, -»• At?,) -sint>,= 2sm-^cos \v, + —j 

= 2 sin — ^sin {«, + J^(7r + Av,)}. . 

. Av^ . I AvA 
A cost?, = cos(v, + A«,) - cost?, = - 2sin— sin It?, + — j 

= 2 sin cos {t?, + i (^ + At?,)}. 

2 



2sin-] sin|a?0 + a+-(7r + 0)}, 
2sin-J cos{a?0 + a + -(7r + 0)}. 



17. To find the differences of tan t?, and tan'^t?,. 

sin t?,+i cost?, - cost?, + i sin t?, 

A tan t?. = tan t?- a. 1 — tan t?, = — — . 

^ wui C7, taxi 1/,+! « ^^g ^^^^ cost?, 

sin(t?,.n-t?,) sin At?, 

" cos t?, + 1 cos t?, tos t?,+ 1 cos t?, 

sin 



Hence Atan«r0 = 



cos (a? + 1)9 cos wO* 

1 t?,+i-t?. 



Also Atan-^t?, = tan-^t?,4.i-tan-»t?, = tan- ^^^ ^^^ 



tan^ 



Hence A tan"^ wO = tan" 



1 +«*+!«, 




1 + (a? + 1) . 
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Relations between the successive Values and the Differences of a 

Function. 

18. To express A"t^« by Ug and its n successive values, 

We here take h instead of unity for the increment of the 
principal variable, as the investigation is precisely the same 
on either supposition. 

Aw, =«^,+A- w. 

Now suppose this law of the coefScients, which as far as 
we have gone is the same as that of an expanded binomial 
whose index is the order of the difference, to hold for the 
ffi^ difference, so that 

then A"+^w,= 



u 



» «*, + («+ 1)*- + Pi) «*, + «*+ (Pl+P2)w, + (n-l)*-.- 
^ 0+ Pi) «*, + /. T«^„ 
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which is the same alteration with regard to the coefficients as 
occurs in passing from (x - l)" to (« - 1)"+^ If therefore, 
for any value of n supposed a positive integer, the coefficients 
of the expansions of A"t«« and (z — 1)" are the same, they will 
always be the same ; but these coefficients are identical when 
n « 1, 2, 3 ; therefore they are always the same ; 

or, supposing A » 1, 
A« n(n-l) 

Obs. For the sake of being easily remembered, this 
result may be written 

where, instead of the powers of u, from n to 0, the » + 1 
descending values f^«+«, t««4n-.i, ... «^«+i, t«« are to be written. 

19. If in the series just investigated, we assign a par- 
ticular value to Ugj we shall readily obtain an expression for 
its n*^ difference. Thus let u^^af^y 

.-. A"(/»") = (a? + n)"-n(<»+n-l)"'+^^^^(a? + n-2)"-&c., 

in which equation, if n > m, since the former member vanishes, 
the second member is zero for every value of ^ ; and if m s n, 
so that A*^ aa 1 . 2 . S ... n, we get 

l.2.S...n«(a7 + n)"-n(a? + n-l)''+ —^^ ^(«» + »-2)"-&c.; 

and making /» » 0, since the equation holds for all values of a?, 

l.2.3...n«n"-n(»-l)* + ^^^— ^(n-2)»-.&c. 
^ ^ 1.2 "^ ^ 

2 
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20. If af«o, and A"0" denote the particular valune of 
A'a?* when w « 0, we have 

A«0" = n" - n (n - 1 )* + ^^^^^^ (n - 2)" - &c. 

Of the numbers comprised in the form A"0", we shall 
make considerable use in future investigations; whenever 
n> m the value is zero, in other cases it may be computed 
by the above formula; thus, 

AG" « 1, 

A*0» = 2, A*0» = 6, A'O* = 14, . . . 

A'O* « 6, A'O* « S6, A'O* = 150, . . . 

21. Reversing the order of the series in Art. 18, we 
find 

(- 1)" A"w, = tt, - «w,+i + ^ ^ w.+t - ... ± w,+,. 

Also putting Ug^oT^ and then supposing iT^^n, <rsl, 
we find successively, 

(-l)"A"a^ «ar-n(a?+l)- + ^^^^^^\a?+2)'"-...ir(a?+ 
(-l)«l.«.3..«=ar-n(ci? + l)" + ^^^^^^(«r + 2)*-...±(a? + «^^^ 
(-l)"l.2.S..n«l"-«.2" + ^ .g**- ... ±(n + l)". 



22. Hence it may be proved that 1 . 2 . 3 ... (p - 1) + i 
is divisible by p, if p be a prime number. 

Let «, = cT* - 1 ; 
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.-. A"Oi?"-l) = l.2.3...n = (j7 + »)*-l-»{(a?+«-l)»-l} 
+ — ^ {(a? + n - 2)« - 1 } - &c. 

1*9! 

Let dj = 1, and n + 1 <> p, then « -t- »>=p, and 
1 . 2 . S ... (p - 1) + 1 = p« - (p _ i){(p -: I)?-* - 1 } 

(p-l)(p-2) 



1 .2 



|(p_2),-i_lj_&c. 



Now by Fermat^s Theorem every term of the second 
member is divisible by p when p is a prime number; con- 
sequently 1 . 2 . 3 ... (p - 1) + 1 is divisible by p. 

23. To express w,^^;^ by w, and its first n differences. 

We take h instead of unity for the increment of the 
principle variable, the investigation being precisely the same 
on either supposition. 

^*+8** ^* + A w, + A (w, + Aw,) = w,+ 2 A«, + A*w,, 
«*x+3*=«**+2Aw,+ A'w,+A(«, + 2A«, + A'tt,) 
= w, + 3Aw, + 3A*w,+ A'w,. 

Now suppose this law of coefficients, which as far as we 
have gone is the same as that of an expanded binomial, 
whose index is the number of increments which the principal 
variable has received, to hold for n increments, so that 

«*,+!.* =w.+PiAtt,+p2A*w,+ ...+PiA-^w,+A"w„ 

then tt,+(^+i);i = «*,+PiAw,+P2A*w,+ ...+Pi A*"*w,+A*w, 

+ A(w,+PiAw,+ ...+p2A*"^w+PiA''"*w,+A"«^J 

= t^, + (1 +Pi)Aw,+ (Pi +P8) A'w,+ ... 

+ (P2+Pi)A'-'w,+ (p, + i)A"w,+ A"+*w„ 
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which is the same alteration with regard to the coefficients 
as occurs in passing from (l + «y to (1 +»)•+*. 

Hence if the coefficients of the expansions of u^^^j^ and 
(l + xy are the same for any value of «, supposed a positive 
integer, they will alwayar be the same; but they are identical 
when n^ly2y 8; therefore they are always the same ; 

1 .2 
or» supposing A « i, 

«*+. = «,+« A«,+ ^'Z ^ AX + ... + nA-^w,+ A*tt,. 
1 • 2 

Obs. This result in order to be more easily remembered 
may be written 

«.+« = (l + A)««„ 

each term of the development of (l + A)* being understood 
to be prefixed to «,. 

24. Let u^^aT, then w.^., = (a? + w)" ; 

.-. (^ + w)"«(l + A)"a?", 

and making j? = o, «* = (l + A)"0*, 

each term of the development of (1 + A)" being prefixed, 
as said above, to oT and 0*, respectively. By the latter for- 
mula any power of a number is expressed by the numbers 
comprised in the form A*0". 

26. To deduce Taylor's Theorem from the formula 

n(n- 1) „ 
^*+n* = «** + w Aw, +— 1__/ AX + &c. 
1 .2 

Let nAs^, and let h be infinitely diminished whilst t 
remams finite; therefore n i» infinitely increased; and since 
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h is indefinitely diminished, we have, regarding the differ* 
ential coefficient as the limit of the ratio of the simultaneous 
increments of the function and the variable, 

Hence, preparing the formula as follows, 






Aw, nh(nh-h) A*«, 
IT"^ Ti A« 

and taking the limit of both sides by supposing A to be 
infinitely diminished and n infinitely increased, their product 
always remaining equal to a finite magnitude t^ we get 

e 

w,+t - «** + td^^, + __ d* w, + &c* 

The Differential Calculus is a particular case of that 
of Finite Differences; and the above investigation is intro- 
duced to shew how, from results in Finite Differences ob- 
tained with an indeterminate increment for the principal 
variable, we may pass to the corresponding results in the 
Differential Calculus. 

26. The theorems of Arts. 18 and 23 have been proved 
by an inductive process; they may also be established by 
the theory of Generating Functions, the principles of which 
we shall now proceed to explain ; as it is a theory which, for 
its generality and power, especially merits our attention. 

Generating Functions. 

27* Let (p{t) be a function of t susceptible of the de- 
velopment 

(p{t) « ... + W.i#~* -k-Uo + Uit -^ ... 

then Ug may evidently represent any function of of whatever, 
if we regard this equation as the definition of <p{t). The 
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function (p(t) consequently by its development generates the 
coefficients Uq, Ui, ...u, annexed to their proper powers of t, 
and is therefore called the Generating Function of Ug^ and 
is denoted by Guji, so that 

Thus since 

log ( 1 - " ' = ^ + i ^ + i ^ + — + - ^ ' + • - 

w 
Similarly^ since 

^(1 -^)-*«# + 2^ + 3/'+... +arr + ... 
^(l-^)-«=Gcr. 

28.^ To determine the generating functions of u,+^ and 
A"t«« from that of u,. 

Let (p(t) = Gugy 

then (p(t) = ... + w,^ + u^+\f"^^ + -. + w,+n^^"+ -. 

.-. r«0(O * G«^,+n, or t-*Gu,^Gu,,,. 

Again, f*0(O " ••• + ^-n + «*-.«+i^ + ... + 1*,-.»^ + ... 

.-. ^0(0 = GUg.^f or f'Gt^, = G«,.,. 
Hence it follows that the generating function of 



Aw« or 



for this function being developed will produce the difference of 
two series whose general terms are respectively w,+i^ and t*,^; 
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.-. G(A«.)=(^-i)g«.. 
Similarly, 

and G(A"w,) « fi - 1 j Gw,. 



Also 



G (A««,_.) - (i - 1)" <?«.-. - (^ - i) Vg«,= (1 - 0'G«.. 



29. To investigate the expression for A"t^« in terms of 
^«+ii9 ^«+ii-i9 &<^*9 ^y Generating Functions. 



G(A»«,)-(i-i)*Gw, 



= r«Gu, - nr»+' Gw, + ^^^ ^^ r»+» Gw. - &c. 

1.2 



Gw,+, - wGw,^.«.i + V Gu^^^^^ - &c. 

1 . z 



^, »(n-l) „ , 

G (tt,+, - »«*,+._! + p^ W.+._2 - &C.) 



« (w - 1) 
*. A»w, = w,+n - w«*+.-i + — T-z — «**+«-.8 - &c. ; 

1 . z 



for, the generating functions of both being the same, the co- 
efficients of /' in their developments must be identical, however 
those developments have been effected. 
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30. To investigate the expression for ^,4.. in terms of Ug 
and its first n difFerences, by Generating Functions. 

/I \ ^ n(n-l) /I \* 
= Gw, + n f- - 1 j Gu^ + ^ g ( - - 1 1 Gu^ + &c. 

= Gu^ + nG(^u,) ^. ^Y g g(A'u^) + &c. 



M Tfi >— 1) 
= G (w, + n A«, + — ^ ' A^u, + &c.) ; 

n(n-l)_. . 
'• «»+•«= »» + n Aw, + — 7—- — A*w, + &c. 



1.2 



31. It is obvious that by transforming the expressions 

[- - iYgu^j and t'*Gu^j 

In different ways, we may obtain various other expressions 
for A"««, and u^^^ besides the above. 

Thus to express A"«, in terms of A"w,.,, A*+'w,_«-.i, 
&c., we have 






(1 - 0"G^«** + «0 - ty+'Gu^ 

+ ^ ^ (1 -0"^ fif«^^ + &c. 



1.2 
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= GiA"«,.,) + nG(A-+'«.-,_,) 
n (» + 1) 











' 1. 


S 


= G 


lA- 


'w»- 


» + W 


n(n 


-ii-i 
+ 1) 










.2 








.-. 


A"«. 


= 


-1- 


n/> 


,n+l- 


/ 


n (n + 1) 



G(A«+««._._,)+&c. 



A"+»«,_,_s + &c} ; 



A"+X-«-« + fec. 

1 . % 



32. Again, to express u,^„ in terms of w„ Aw,-rj 
A*w,_9y, &c., we must transform ^"" into a series of powers 

of r I 1 j , or develope ^"^ in powers of a from the 

equation 



1 



-ii)'' 



which may be done by Lagrange's Theorem ; and we find 

w(n + 2r-l) 
1 . % 

1.2.3 

(HerschePs Examples). 

Obs. This method is obviously not confined to the 
function w,+, - w,; it is equally applicable to any other com- 
bination of the successive values u^, ^«+i» Ug+ji, &C.9 of the 
first degree. If we take Aw^ to mean ow,+8 + 6i*,+i + c«*,, 
then the generating functions of Au^, and A"t^« will evidently 
be 

and the expression for A"t^, in terms of u, and its successive 
values, might be obtained as in the preceding case. 
3 
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Separation of the Symbols of Operation from those of Quantity. 
33. We have seen (Art. 23) in the formula 

an instance of the system of notation, which consists in sepa- 
rating the symbols of operation from those of quantity ; the 
use of which is not confined to simple cases like that just 
noticed, but may be extended with remarkable effect to a 
great variety of investigations connected with this subject- 

If the expression (l + A)% regarded as a function of a 
certain symbol A, be expanded in powers of A, it will pro- 
duce the series, 

w (n - 1) ^ „ 
1+ n A + — — ^— A^ + &c., 

and if the symbol A be looked upon merely as an instrument 
by means of which we are enabled to produce the numerical 
coefficients of the series affected with their proper powers 
of A, then the expression (l + A)" .must be considered as 
having no other meaning than as an abbreviated expression 
for its development; and when prefixed to the function w„ 
each term of this development is understood to be applied 
separately to that function ; so that the following expressions, 

(1 + A)"w,, 

^ w(n -1) . , . 

(1 +wA + — ^ ^ A + ...)u„, 

u, + nAu, + — - A*««j,+ ..., 

1 « <^ 

may be used indifferently for one another, the two former 
being regarded merely as abbreviations of the latter. 

34. In general, if ^(A) be a function of A capable of 
being developed in a series of powers of A, such as 

^ A« + B A^ + &c. ; then for Aj\''u^ + BA^u, + &c. 
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the expression F(A)Uj, is used as an abbreviation; and the 
same notation is applicable to other characteristic letters, such 
as d,, ^, 2, Ss9 &c. 

Hence also, the successive performance of two or more 
series of operations represented by F (A), JF" (A), &c. upon 
the same function u,, is equivalent to the performance of 
that series of operations which is represented by their product. 

For example, we have seen that 

A (w,t?,) = A w, . tj,+ w, AtJ,+i = (A + At?) u^v^ 

supposing A to affect t^« only, and using At? for an operation 
affecting v, only, such that (Atj/t?, = A'tJ^+ri 

then A^ (w,t?,) = (A + At?) A (i*,t?,) = (A + At?)*«,tj,; 

and, generally, 

A"(w«t?,) = (A + A »)"«,<?, 

= A'w, . tJ, + n A*~^ w, At),+i + — A""*w,A- tJ^+8 + &c. 

1 . 2 

which may be also proved inductively by shewing, as in 
Art. J 8, that the coefficients of the developments of A'*(w,t?,) 
and (l -i- xYj which are identical when w = 1, undergo the 
same changes in passing from n to n + \. 

Using an accent not to imply that the operation denoted by 
A is altered at all, but ^merely that A' affects v, only, whilst 
A affects u, only, the above results are sometimes written, 

A (w,t>J = (A + A' + A A') w,t)„ 

A' (w,t?,) =. (A + A' + AA')'«^,t?„ 

and, generally, A" (u^vj) «= (A + A' -h A A')" u^v^; 

or, if there were more functions w^^ x^y &c. and we use 
A", A'", &c. to imply that these symbols only affect «?,, »^, &c. 
respectively, we have 

= {(l + A)(l + A')(l + A")(l + A'")...-ir^-^*"'*** — 
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35. To shew that A"«, «= (e**- - !)"«,» in which the sym- 
bols of operation are separated from those of quantity. 

By Taylor's Theorem, we have 

•»* ^ 

1 . 2 

and separating the symbols of operation from those of quan- 
tity> we get 

^ «V n'di „ , . 
^ 1.2 1.2.3 ^ 

Similarly, «,+„-i= e^""^^**'!*,, 



w,+i = e^-w,; 



1 . 2 

or, again separating the symbols of operation from those of 
quantity, 

A"«, = (e^- - we<»-^>^' + —^^ ^ c^--*)"- - &c.) w, 

- (e^- - i)«w,» 
a celebrated theorem first given by Lagrange. 

36. To find a general expression for A'w, in terms of 
Ug and its differential coefBdents. 

The development of the second member of the equation 
A>, = (e*'- - 1)X 
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will consist of a series of terms of the form 

(^0 + A^^ + ••• + XC + •••) «^*» 

or AqU^ -\- Aid^u^ + ... ^^„<w, + ..o 

and A„ is evidently the coefficient of t* in the expansion 
of (e'-l)". 

Now (e« - i)« = e^ - nc<-^>* + ^^»^ ^(._,), _ ^^ 
^ ^ 1.2 

and the coefficients of /" in the developments of e*^ e^""*>*, &c. 
are respectively 

"" (^-J)" fee- 
• . J« =s ,— In" - n (» - l)» + -5^ ^ (n - 2)- - &c.| 



A"0" 



(Art. 20). 



Now so long as m Kn^ this vanishes ; and when n^m^ 
A*0"«|n; 

A"0"+* A"0"+" 
.-. A"w, = cCw* + rf;+^w, + -dl+*w, + &c. 



n + 1 



w + 2 



Ex. Let u, = ar, then 
A"(*r"') = w(f»-l) ... (m-n ^ l)<r"»""" 

A"0"+' 



WH- 1 



m(w- 1)... (m-n)a?'"-''-' 4- &c. 4- A"0". 



37. If in Lagrange's Theorem for A"w„ » = 1, we have 
Aw,= («'*'- l)«^,9 or A = c"*'- 1, the meaning of which is, 
that the operation denoted by A is equivalent to the series 
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of operations denoted by e"*' - 1. And, generally, the series 
of operations denoted by /(A) is equivalent to that denoted 
by f{e^' - 1). For let /(A) be developed in a series of the 
form 

/(A) = ^ A« + 5A^ + CA^ + &c., 
then /(A)w, = ^A«w, + 5A^w,+ CAyw,+ &e. 

= J(e*''-l)«w,4-B(e'''-l)^w,+C(e'''-l)rw,4-&c., 
or, separating the symbols of operation from those of quantity, 

38. Suppose /(A) = (l + A)", 

then f{e^' - 1) = (l + e^' - 1)* = e^- ; 
.*. (l + A)*w, = ^"^a = «^*+iij as already proved. 
Again, if /(A) = {log(l + A)}-, 

f{e'' - 1) = {log(l + c^'- !)}•= (logc-0"= «; 
.-. JIog(l + A)}»t^, = cCt^„ 

a formula by which the vH^ differential coefficient of a function 
is expressed by its differences. 

39. To find the general term of the expansion of /(c*) 
in a series ascending by powers of t. 

Writing /(e*) in the form /{i + (e«- i)|, and expand- 
ing by Taylor*'s Theorem, we find 

/(eO =/0) +/(i) C*'- 1) + n^/'O) («'- 1)'+ ... 
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Then taking, as in Art. 36, the coefficient of t^ in each 
term of the second member, and observing that in/(l) it may 

/(l).O'" 
be represented by — j , this quantity being 1 when «i = 0, 

and zero in all other cases; and that in 

1 /»Hi) (e'-iy it is /-/"Ho ^"^"^ 



\n \l In 

we have for the coefficient of f in the expansion of /(e*), the 
value 

^« = A {/(00"+/(i) A0-+ ji^/'(l)A'=o»+ ...} 

-[|/(I.A)0-, 

a remarkable theorem first given by Herschel ; for the applica- 
tions of which, see his Collection of Examples. 
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SECTION II. 

INVERSE METHOD OF DIFFERENCES. 



Integration of Explicit Functions. 

40. The Inverse Method of Differences has for its object 
to determine the primitive function from its given difference ; 
or from given relations between it and its differences. We 
shall begin with the simplest case, 

in which it is required to determine a function whose differ- 
ence is given explicitly in terms of the principal variable. 

41. Since ^u^ is the difference of w, + C, as well as of 
Ugf it will be necessary, in passing from the given difference 
^Ug to the primitive function, to annex an arbitrary constant 
C, in order to give the result all the generality of which it 
is capable. Also C may be a function of a? as well as an 
arbitrary constant, provided its value remains unaltered whilst 
w changes to 07 + 1. For if C« denote such a function of w 
that Cg+i = C,, or AC, = 0, we shall have 

A (w,+ C,) = /\u,. 

It is evident that C, ■= ^(2X7r<») has the property in 
question, ^ denoting any trigonometrical function, sine, co- 
sine, &c., and X any integer. 

42. The symbol IS is used to denote the operation by 
which we pass from the difference At^^^. to the primitive 
function; so that 

2 (AWx) = Ug + constant. 
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Also, as the same function admits of successive differ- 
ences, so a function may be integrated any number of times ; 
the second integral of w,, or 2(2w,), is written 2*w,, and 
the n^ integral 2"w,. 

We now proceed to deduce the integrals of Various ex- 
pressions ; chiefly, by reversing the processes given in Section 
I. for finding the differences of functions. 

43. It is evident that 2 (w, + 1>, + w^) « 2 w, + 2t>^ + 2w, ; 
for if we take the difference of both sides, we get the same 
result, viz. t^, + t9, + tr«. And in the same manner it appears 
that 2 (at^,) = a2«,. 

44. To find the integral of any rational integral function. 

Since the difference of a rational integral function is a 
function of the same kind one dimension lower, it follows 
that the integral of a function of that description is a similar 
function one dimension higher; hence, to find the integral of 

we may assume it equal to 

o^+» + haf* + ... 4- Ap^ + /; 

then upon taking the difference of each side, and equating 
the coefficients of like powers of <v, there will arise n + 1 
simple equations to determine the n + 1 quantities a, 6, e,...Ap; 
the last term I will remain indeterminate, being in fact the 
arbitrary constant which must be added to make the integral 
complete. 

Ex. To find 2(^*+l). 

Assume 2 (^ + 1) = a^ + ba^ + car' + da^ + c^ ; 

.% ^ + 1« o(5^*+ 10a?^ + 10a?^ -f 5cr+ 1) + 6(4«?^+6a?* + 4a?+l) 

+ c{Sa^ -{^Sw + 1) + d(2a? + 1) -k- e; 
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ls:a + & + c + d + «- 



1^1 ^ ^ ^ ^9 

5 2 3' so' 



^ a?* 0?^ 29 
^ ^ 5 2 3 30 



46. To find the integral of the product of consecutive 
terms of an arithmetic progression, we must annex one more 
factor at the beginning, and divide by the number of factors 
so increased and by the common difference. 

For let w^= a + 6^, then we have seen, (Art. 11.) that 

Aw,w,+i ••• w,+« « w,+i w-r+2 ... w,+« . (n + 1)6, 

therefore, taking the integrals of both sides, and writing »-l 
for a?, we get 



(n + 1) 
which proves the rule stated above. 



+ c. 






If one or more factors be deficient in a factorial of this 
kind, it may be resolved into others which are complete, as 
in the following instance; 

(2zp + 1) (2fl7 + 5) (2a? + 7) = (2a? + 3 - 2) (2a? + 5) (2a? + 7) 

= (2a? + 3) (2a7 + 5) (2a? + 7) - 2 (2a? + 5) (2a? + 7). 
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. 46. A rational integral function may often be resolved 
into factorials of the above form, and in this way its integral 
more conveniently found, than by the method of Art. 44. 

Ex. 1. a^ + a^^a^{w + 1) ^ (w-1 + l)w(w+l) 

= (^ - 1) J? (a? + 1) + a^ (a? + 1) ; 

4, S 

And in general, any quantity of the form 

aaf" + 6^" * + co^"* + &c. 

may be resolved into factorials, by the method of indetermi- 
nate coefficients; thus, if we assume 

o/p* + 6a? + c = -rf (^ + 1) (a? + 2) + S (a? + 1) + C, 

making zr = - 1, we get a - 6 + c « C; 

.-. a (/»* - 1) + 6 (^ + 1) « -rf (« + 1) (^ + 2) + S (a? + l), 

or a (a? - 1) + 6 = -4 (^ + 2) + S ; 

make ^ = - 2, .-. - Sa + 6 = S ; 

.*. o (^ - 1) + Sa = -4 (^ + ^), .". J ^a. 

In practice, however, it is generally easier to resolve a 
function by inspection, as in Ex. 1, than by this method, 
which is theoretically certain. 

47. To find the integral of a fraction whose denominator 
is the product of consecutive terms of an arithmetic progresr 
sion, and numerator constant, we must efface the last factor, 
divide by the number of factors remaining and by the common 
difference, and prefix a negative sign. 



Digitized by LjOOQ IC 



For let ti^ at a + bar, then we have teen (Art. 12.) that 

c nbc 

A = ; 

r. taking the integrals of both sides, 

which proves the rule just stated. 

48. If the proposed fraction, instead of having its nu- 
merator constant, be 

(the degree of the numerator being at least lower by two units 
than that of the denominator,) we must reduce the numerator 
to a series of terms each of which is the product of consecutive 
factors reckoning from the beginning of the denominator ; that 
is, assume 

-4a^-* + So^"* + ... + Jra? + Zr= A' + ffu, + Cu,u,+i + ... 

+ K UgfUff^i ... t^jr-fa-39 

then, developing the second member, and equating coefficients 
of like powers of ^, we obtain n — 1 equations for determining 
A\ -B*, C ... JT'; and the fraction resolves itself into the fol- 
lowing, each of which is integrable, 

+ &c. + - 



If the degree of the numerator were only lower by one 
unit than that of the denominator, we should arrive at a term 

, of which we are able to find the integral, only ap- 

proximately. 
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Hence also, if any of the factors of the denominator of the 
fraction in Art. 47 be wanting, they may be supplied by in- 
troducing them into the numerator and denominator at the 
same time; and then the resulting fraction may be treated 
as in the present Article. 

Ex. _1 (^-1)^(^ + 1) 

J?* - 4 " (a? - 2) ... (a? + 2) 

1 3 6 



{a + 1) (a? + 2) w{w-\- 1) {a + 2) {x - 1) {at + !)(« + 2) 

6 

^ (^ - 2) (a? - 1) ar (^ + l) (a? + 2) * 

which is got by assuming 

(a?-l)^(^+l) = a(a?-2)(^-l)^ + 6(«r-2)(^-l)+c(a;-2) + rf, 

and making ^s2, 1, 0, successively; taking care to reject 
the factor common to both sides, after each substitution. 

49. To find the integral of a'. 
We have seen (Art. 14.) that Aa'«(a-l)a'; 

... 2o'--^+C. Also S-o'-- ° 



suppressing the part introduced by the constants, which would 
be a rational integral function of the (n - 1)*^ degree. 

50. To find the integral of log v«. 

If W, « log (ViVjVs ... t>,_i), 

then Aw, «log(i?ir8...v,) - log(i?it>8...tJ,_,) «logt>,; 

.•. 2 log Vg^Ug-\- log C = log (Ct?! Vj ... v,.i) = log CPvg^xj 

using Pvg to denote the product of all the successive values of 
the function v«, from some fixed term Vj, (or more generally 
v^i n being independent of w) to Vg inclusive. 
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61. To find the integrals of cos<v0, sin w9. ^ 

9 ffj^ 

Since A cosa?0 = - 2 sin - sin (ar + ^) 0, 

. ^ 
.-. A COS (a? - ^) « - 2 sin - sin w9; 

I 

.'. S sm tpfl = - J + C 

2sin- 
2 

Also, since 
A" sin (w9 + a) « (2 sin - 1 sin { a'0 + a + - (7r+0)} , (Art. 16.), 

n 
integrating both sides n times, and replacing a by a - - iir+6), 

we get 

2" sin (xO + o) = ^ 

(2 sin-) 



Again, A sin a;0 = 2 sin - cos (w + ^) 0, 

Q 

.'• A sin (a? - ^) = 2 sin - cos 0?^ ; 

^ sin (a? - i) ^ 
.-. 2cosa?0 = i 9^^' 

2 sin - 

• 2 



ft 

cos{a?0 + a-- (7r + 6)} 



and 2* cos («»0 + a) < 



(*^'"i)" 



Digitized by LjOOQ IC 



31 

5ii, The preceding expressions may also be integrated by 
substituting for them their exponential values; as in the 
fol 1 ving instance. 

. acosfo?- 1)0- cos/pfl ^ 

= a* r^ — + C. 

or - 2a cos + 1 

Hence, putting a? ^2a cosO + 1 '^c, and denoting 
a* cos w0 by «^, we have cSw^, = «*««-i - «*r> 

and generally, 

1. • ss 

Exactly the same formulae hold for u^^af AnwO. 

In the same manner the integrals of a' (sin wff)^, 
(f(cosa9)^ may be obtained. 

63. To find the integral of ; —r. 

** cos wO cos (Of + 1)9 

^* ^ ^ sind rA . 

Since Atand?d= — -, (Art. 17.) 

cos w0 cos (a? + 1) & 

^ * 1 tan^O •, 

we have 2 ^ 7 TTS " . ^ + C. 

cos ^0 cos (^ + 1) & sm 
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64. To find the integral of tan 



-1 



p-^qw •\-Tai^ 



Since Atan-Ha + 6a^)=tan-^ ; —— -, (Art 17.) 



,.i 1 



we may assume 2 tan" -i «= tan"* (a + 6a;), 

p +^a? + rar 

and take the difference of both sides ; then if the proposed 
function is capable of being integrated, the indeterminate 
coefficients a and h will become known. 



56. The integral of o'w^w^+i--. Wj,+„-i5 where Ug^p(f-\-qy 
may be determined by assuming it equal to the same ex- 
pression (only with another factor at the beginning instead of 
a') multiplied by an indeterminate coefiicient; for 



So'M^... «,4.,,i = ^ ^ . ^j-l«, ». ^^+,, 



a* 

Similarly, to find 2 , for the assumption 

t^,w,+i ... w,+« 

we must efface the last factor in the denominator, and write 
instead of a* an indeterminate coefficient; for 

1 ^^4.n - u^ pa'ia^-i) 



•. 2 



a* 1*1 



«^««j+l ••• «^*+n P (a" - 1) «^*W*+l ... «^*+ii-l 
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56. In like manner, if t«« « a + bw, expressions of the' 
forms 

(p + qw) f (p -^qw -^ ra^) f 

can sometimes be integrated, by assuming their integrals equal 
to expressions of the same form, except that the last factor in 
the denominator is effaced, and the polynomial in the nu- 
merator is replaced by another one dimension lower with 
indeterminate coefficients. It is of course only when a certain 
equation of condition between the quantities a, 6, p, g, t is 
satisfied, that this method succeeds. 

Ex. 1. Let 2 7 -7^ r: — = ^^' ; 

(2cv - 1) (2^ + 1) 3' %m "X 



(2 



a?-l) (2ar+ 1)3* ^^^ \%(b^\ Qw-l) 

^3; (2a?-.i)(2a?+l)' 



.-. -rf=:-f, and 2w,-:C-| 



Ex. 2. Let 2 



(2a? - 1) 3' 
a;^ + 6cV + 12 -4 + Bar 



0? (a? + 1) (a? + 2) 2* a? (a? + 1) 2* 



then -4»-6, S«-2, and 2w, = C- 



a;(a?+l)2— ^' 

67* Since A (u^vj) = w,At>, + t>,4.iAw,, we have 
2 («^, At?,) = w*t>, - 2 (i?,+, AwJ, 

the formula for integration by parts, corresponding to the 
formula 

58. Change v^ into 2t7,, and consequently Af, into t?,, 
then 2 («,«,) = w, 2«, - 2 ( A w, 2 1?,+ 1) ; 
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h^nce, by sucocasive substitutions, we get 

2(A«. 2».^0 = A», 2*»,+, - 2(A*«.2'»,+,), 
2(A*».2'»,4,) •= A'«,2?«,+, - 2 (A'«,2?«.+3), 



2(A»«,S'»,+.> - A*»,S'"'*»*+--2(A'+'«,2"+*»,+.+,); 

.-. 2 («,»,) = «,S«, - A «,2't>,+ , + A*tt,2'»,+» 

- ... ± A»«,2"*'»,+„ =F 2 (A-^'«,2«* '«,+.+,)• 

Obs. In order to be more easily r^uembered, this for- 
mula may be written 

2(»,«,) = (1 + A2«)-'»,2»„ 

where A affects u, only, and 2v affects 2v« only, and implies 
an operation such that 

(2«)'2«, = 2'+««,+,. 
Hence, 

2*(«,«,) « (1 + A2«)-' 2 («,2«.) - (1 + A2«)-*«.2*r,; 
and, generally, 
2" («.»,) - (1 + A 2«)-"».2"», = »,2"«, - « A«,2"+'«,+, 

+ ^^?^A*«,2«+'«.+,-&c., (1) 

1 . Z 

which may be also proved inductively, by shewing that the 
coefficients of the developments of 2*(t*,«,) and (l+»)"% 
which are identical when n » 1, undergo the same changes in 
passing from n to n - 1. This appears by differencing both 
sides of (l). 

59. The above formula always enables us to find the 
integrals of functions made up of two factors, one of which 
leads to zero, as the value of one of its successive differences, 
and the other admits of successive integrations. Suppose, 
for example, that u^ is a rational integral function of the 
nf^ degree, and that «, = a' ; then 

2t>,= , 2^tJ, + i=- -T-, &c.; A"«^,=const, A"+*w,=0; 

a — 1 (a — 1} 
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Again, suppose u, to be a rational integral function of 
the n*^ degree, and «, = cosd?0; then 

cos {ar0--(7r + 0)} 

2 («, cos xd) = I ' 1 }^ 

2 sm- 

2 

the series terminating with A*w,. Similarly, if «, = a' cos wO; 
or if t>, = a' cos*" a?6 sin" w0j since the product cos* x9 sin* a^0 
may be replaced by simple dimensions of sines and cosines of 
multiples of w6; and it will be noticed that the fraction of 
Art. 48 may be brought under this case. 

60. Since the perfoi*mance of the operation 2 upon any 
series of terms -4A*w^ + SA"tt, + ..., reduces it to 

it appears that prefixing 2 to {-^A" + -BA* + ...}w, has the 
same effect as prefixing A"^; in other words, 2 is equivalent 
to A~^ And in like manner, since integrating A"w, n times, 
reduces it to A*""t*,, S" must be equivalent to A"*. 

The same reasoning is applicable to the symbols f*^ d~"; 
whenever therefore, in separating the symbols of operation 
from those of quantity, as in the expression jP(A)fr«, f(d^Ug^ 
terms containing negative powers of A and d^ occur, tbey 
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must be understood to be replaced by the corresponding 
positive powers of S and j^. This being premised, we pro- 
ceed to investigate a general series for 2w*; preparatory to 
which the following propositions must be proved. 

61. To determine the generating functions of 

2"w„ dlu^, fj'u,, 
from that of u,> 

By virtue of the relations 

G(A«,)- Q-i)g»„ 

we have (--A G (2«,) = G (A2«,) = G«, ; 
.-. G (2«,) - (i - l^Gw,, 

and G (S-m.) =(7-1) ^«»- 

62. Again, since 

•■• — + (»,+»-«,) ^ + ...-(^-IjGtt, 

A* 
= {1 -hlogt + (logty-iiC.-l}Gu,; 

I . S 

therefore, dividing both sides by h and then making k « 0, 
which we are at liberty to do since h is indeterminate. 
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... + | ^*^*/"M ^ + ... - - log tGu^ = 1<« 7 . Ou^y 
I ^ J»=o * 



or ... + dfgUg ./*+... 



',-.Gu,; 



.'. G (d,«,) = log J . Ou^ 
and (3l(d;«,) - |logj) (?«,. 
Again, log - G (JC«,) = G (d.j&w.) = €?«, ; 
.-. «(/.«.)-= (log J) «»-> and G (/>.) = (log i] "(?«,. 

63. To investigate a general series for 2u«, involving 
on]y ffUt, »,, and the differential coefficients of «.. 

G (2»,) = (j - l) G«. = (e"*' - ly'Gu, 
logj 



logi 



e '-1 



(logi)'G«.{l-ilogA + A(logi)'-&c. 

.(-ir'^(iogi)%...h 



assuming, as will be proved in the next Art., that 



'-I 



can 



be expanded in a series of the same form as that within 
brackets, and denoting by 



B. 



1.2' 1.3.3. 4 



-^, &c., (-1^1*^-^ 



2n 
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the coefficients of «*, «S ... o^ in that expansion. Hence 
G (2«.) - G (jt«,) - ^ G«, + A G (d.„,) _ &c. 

.B B 

1.2 [4 






Ex. 1. Let w, = d^; 



mMI + I 

wiH- 1 * * 



a^ a^ (J? X 

Hence Sa^* = + + C; 

5 2 3 30 



6 2 12 12 



Ex. 2. Let u. = 



o + 6d?' 



1 


--log(o + 


hx)- 


1 
2(o + 


6») 




o + 6j? 






5,6 


- + - 


fisi* 


^. - &c. 


+ C. 



2(o + 6a?)' 4(a+fta?y 
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Numbers of BeriKmilli. 



64. The numbers J?i, j?39... which are required in the 
general value of ^u^ in the preceding Art., are called the 
numbers of BernouiUi^ and are of great importance in the 
theory of Series. They are defined by the equation 

a l^lt + —^f = <* + ... 

c*-l ® 1.2 1.2.5.4 

and their values may be computed in the following manner. 

t 
The development of -j can involve no odd powers of 

e — 1 i 
' t 

t above the first. For if <p{t) = -^ , we have 

.-. 0(0 - 0(- == - *» an<* ^«*0(O - ^**0(- = 0, • 

or d|*0(O«rft'0(-O. 

Hence, when the development is twice differentiated, it will 
not be altered by changing t into (~ t)^ and therefore con- 
tains only even powers of t\ therefore (p{t) can only contain 
one odd power of t viz. the first. Also since 

t 1 



c*-l t f 

2 1.2.3 

we may assume for (^(ty the form of development given 
above, viz. 

^«i.^^^A^-^v+...+(-ir^^^-i&c. 

e^-l * 1.2 [4 [£!! 
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and if J^^ be the coefficient of ^" in this development of the 
second member, 

^|. «= — * 1 + -ff • i — &c. 

[£^ [2n 

1 (2ny-2n(2n-iy"+... „ 

+ . ^ — r ' T 8ic. 

2n + 1 [2n 

[2n* * ^ 2n + l * 

all the terms after A*"0** vanishing, since A*0** is zero when 
m > 2n. 

Hence, 

*«•.! - (-l)"+M-i ^0««+ ^ A*0*-- &c. ... + — L_ A««0««}. 

By this formula J?i, J?3, J?6, &c. may be readily com- 
puted, supposing the numbers comprised in the form A"0^" 
to be known; we find 

11115 
^'^6' ^'"i5' ^*==i^' ^^'^iS' ^•'66'^'^- 

65, Also, since 

A"(^ = in"" - f» (f» - 1)"" + ^ ^ (m - 2)"" - &c. (Art. 20.) 

we may, if we please, eliminate the numbers A"0** from the 
expression for J?2fi-i« ^^^ ^^ ^^^ 

|2n 
(-l)"+'Bs,.i = -i + J(2»--2)-i(S*'-3.2^+3) + &c.+ -^ 



2n + l 
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Besides serving to express the general value of Stt«» the 
numbers of Bemouilli haive various other uses, of which we 
shall now give one or two of the most remarkable. . 

66. To find the general terms of the expansions of 
cot and tan 9 in powers of ft 

Now the general term of the expansion of 



is 



e2ev=i_i " V V |2^ 



(-i)"^'%^(2^v/-^r. 



or j- 0^; 

\9,n 



.*. the general term of the expansion of cot Q is 

and cot = B^Q !— 0* -- &c. 

Q 1.2 1.2 .3.4 

Also since tan = cot - 2 cot 20, 

the general term of the series for tan Q is 

— i + 2 . i 

6 
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Hence, by differentiating the above expressions for cot 9 
and tand, we may deduce the general terms of the expan- 
sions of cot'0 and tan*0; and by integrating them, the 
general terms of the expansions of log sin d, log cos d. 

67. To find the sum of the series 



1 1 J 
1*» "*" ^' "^ 5 



Ti;r + ^,+ 5;+-<^^^)' 



or the value of ^S* --r- . 



changing into irOy we get 

therefore, differentiating the logarithm of each member. 

But the coeflicieiit of d*"'* in the expansion of w cot ir9 is 

[zn ' 
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U.« 1 2^»-^S«._,7r- 



1.2.3...2n 

„ XX. 9r* 1 1 - 

Hence -h---, + -+...« — , 77-»"::7 + r: + •••- — • 



111 9r* 1 1 1 ir* 



l« + 22+32 + --e' 14 -^24 + 34 + --90" 
Obs. Calling S^^ the sum of this series, we have 
S^, ^ (2n + 1) (2n -t- 2) ^an-i 

B n* 

Now suppose w very great, then '*'*"^ = — , which 

proves the divergency of the series formed by the numbers 
of Bemouilli ; these numbers increase very rapidly, beginning 
with J?i3. 

68. To find the sum of the series 



1 1 1 

i^« ■*■ i*^ "^ 5^ 



;t«+:5; + ^-^ •••(*«»)' 



or the value of ^ST 



(2w -I- 1)*" 

.„. »..-{.-(^)]{.-(a}{-(!-:)}- 



changing Q mto — , we get 
2 



ir9 
COS 

2 



therefore, differentiating the logarithm of each member, 

:-?-¥{- (?)'r^F{-(sr 

+ &C.; 



IT 

2 

81-1 



^l'-(^T 
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and equating the coefficients of ^~' in each member, which 
are respectively, 

U) ^--V"^' * '^(i^TlT' 

Hence -t + - + — + &c. =-.-. « — . 

1* 3* 5^ 2 6 1.2 8 



9. Also ,-t: - -r: +----— -1^ ..- 

22« ^ ^2n ^ 

" !*• ■*■ 3^ ■*" 5^-,"*" ••• " i^ Vl^ "^ 2** "** 3^« ■*" "V 

(2«»'^-l)^"'Jg,,_^ 

|2n 



70. Since log (l + «"*) « e^' - :| c"^* + i^""*" - &c-, 

integrating this 292-1 times between the limits z^O^ x^ co^ 
we find for result the series just summed ; 



71. To find an approximate value of 

T{a} + 1) = 1 . 2 . 3... a?, when a? is very large. 

Making u^ = log Of in the formula 

2w, ^ f^u^-\ w, + --^d,w, - &c. 
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we find, (Art. 50,) %adding log a? to both sides, 
log {l .2.3... (a?- l)tr} = C + a? logo? - a? - ^logo? 

+ ^^ - + ... + log W 

= C + (a? + ^) log a? - a? + log (l + A), 

S S B 

putting log (1 + A) = ^- ^— + — ^ - &c., 

so that A is a quantity continually approaching to zero as w in- 
creases. 

Now to determine C, suppose a? very large so that h may 
be neglected, and change a? into 2a?, then 

log (1 . 2 . 3 ... 2a?) = C + (2a? + ^) log 2ct? - 2a? 

= C+(2a'H-^) (log a? + log 2) -2 a?, 
and log (2 . 4 . 6 ... 2a?) = log (2*. 1 . 2 . 3 ... a?) 

= a? log 2 + C + (a? + ^) log a? - a? ; 
.-. log {1.3. 5. ..(2a?- 1)} =:a?loga? + (a? + ^)log2-a?; 
2.4. 6. ..2a? ^11 1 1 « 

= C + ^ log2a? -log2; 
2. 4. 6. ..(2.2?) 

2 . 2 . 4 . 4 . 6 . 6 ... (2a? - 2) 2a? 
^^1.3.3.5.5.7... (2a?-l)(2a?- l) 

= log — , by Wallis's Theorem, 

since a? is indefinitely large; 
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.'. log(1.8.3....v) = ^log27r + (x + ^) \oga> - a> + log(l + h) 
= log vaTrar + log ( - I + log (l + *) ; 

.-. 1 .2.3...a? = y/zira>. | - J .(l+A), 

where A is to be calculated from the series 

logo +A) = — .--—.-+ — .--&c., 

and in general it will suffice to take the first term only, 
which gives 

1 . 2 . 3 ... .r = V 27r.r afe ^»'. 

Obs. The preceding series, even for large values of a?, 
becomes divergent after a certain number of terms; this will 
happen after n terms if 

(2w + 1) (2n + 2) ' a?^"+* (2n - 1) 2n o^^-^ ' 
^2n^i ^ (2n + l)(2n + 2) ^ 

or ^;; > ; r ^, 

B2n-l (2W-l)2n 

but the first member of this inequality never exceeds 

.-. (2n - 1) 2n > 47r*J7S or n > tt^. 

It can, however, be proved that an approximate value 
of the series will be obtained by taking the aggregate of 
the convergent terms only. 
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SECTION III. 



EQUATIONS OF DIFFERENCES. 



^2, We now come to the case in which the relation 
between the principal variable and any function of it is to 
be determined by means of an equation between cT, w,, and 
one or more of the successive values, or differences, of w, ; 
that is, from equations of the form 

F {w, u^, w,+i, ... z^,+.) = 0, 

or, /(a?, w,, Aw,, ... A" wj = 0, 

since, by the theorems of Arts. 18 and 23, these forms are 
convertible one into the other. An Equation of Differences 
is said to be of the n^^ order when the successive value, or 
the difference, of the highest order which it involves is the 

'JS, The complete integral of an equation of differences 
of the n}"^ order will contain n arbitrary constants. 

Let ... w„ Ug+^i ^s+2h9 ••• ^ ^ series of terms corre- 
sponding to the successive values w, x -\- h, o^ + SA, ...; and 
let 

F (a?, w„ a) = 

be the equation by which the general term is determined as 
a function of a and a, or the equation of the series, a being 
an arbitrary constant. Since this equation must hold for all 
the succeeding terms, we shall have 

F{x + h,u^^^, a) = 0. 
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Eliminating (a) between these two equations, we get an 
equation between a?, w„ and «,+»;. or, substituting w, + A w*. 
for u^^^f an equation between w, w,, and A w„ which is the 
equation of differences of the first order whose primitive 
equation is 

F(aff w,, a) = 0. 

In like manner if the equation of the general term con- 
tained two arbitrary constants a and by as 

F(a}y u^y tty b) = 0, 

we might eliminate a and b by means of the two succeeding 
equations, 

F(^+A, «,^^, ay 6) = 0, F(W'\'Qhy w,+2fc9«j6)=0, 

and thus get an equation between Wy w,, ^x+h^ ^*+2»» ^^» 
substituting 

Ug+Au, for u^^^y and w^ + 2 Aw, + A^w, for w,^^^, 

an equation between Wy u^y Aw*, A^'w,, without the constants 
a and 6, which is an equation of differences of the second 
order, having for its complete primitive the equation 

F(wy Ugy a, 6) = 0. 

Hence it appears that every equation of differences of the 
first order, or between two successive terms, will introduce 
one arbitrary constant into the equation of the series ; every 
equation of differences of the second order, or between three 
successive terms, will introduce two arbitrary constants into 
the equation of the series ; and, generally, every equation of 
differences of the n*^ order will introduce n arbitrary con- 
stants into the equation of the series. 

Linear Equation of Differences of the First Order. 

7*. The general equationr of the first order and de- 
gree is 
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Aa and B^ being functions of oo. To integrate it, assume 

••• «?*+! {w^ + A t«>,) - J, «, w', = S, ; 

and in order that this equation may resolve itself into two 
others, each of which admits of being integrated, assume 
(as we are at liberty to do, having made only one supposition 
respecting «, and w,) 

or, dividing by v^w^^ ^ 

But Alogtr, = log-^, .•. Alogr, = log-rf,, 

.-. logt?, = 21og^ = logP^,.i, (Art. 50.) 

or v, = PA,-i, the constant being unnecessary. 
The other part of the equation gives 

and «. = P^-,{S(^) + C} 
the complete integral, involving one arbitrary constant. 

• 7 
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Obs. Taking the difference of the result 

which shews that -^-j- is a factor which makes each side of 

the proposed equation integrable; and it is generally the 
most convenient way of integrating the equation to multiply 
it by this factor, 

Ex. 1. Ug^i - aUg « «?*.N 
Here J, = a, PA, = a* ; 

.-. ^ = 2 (-) « 2 (4^a0, putting i « a, 



a - 1 (a - 1)» (a - !)» 

2«+l 2 



+C;(Art.59.) 



' 1 - a (1 - o)* (1 - o)' 

Ex. 2. tt^+, - aug = cos a?d. ^ 

cos(^i7 — 1)0- a cos J70 _ 

w, « — ij 7{ + Co*. 

a^-2acos0H- 1 

Ex. 3. Two vessels which hold (a) and (6) gallons re- 
spectively are filled, the one with proof spirit, the other 
with water ; (c) gallons are taken from each and poured into 
the other; and this is repeated such a number of times as 
to make their contents of the same strength ; find the number 
of times. 
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Indirect Integrals of Equations of Differences. 

^S. Since the equation of differences of the first order, 

./(^, w^, Aw,) « 0, 

is formed by eliminating the constant (a) between the equations 

w, t^ F(ay a), u^^j^ - F(af + A, a), 

it follows that we shall arrive at the same equation of differ- 
ences, whether a be constant, or be a function 6f>r, as a,, 
provided it be such that 

F(af + A, o,+ J = F(m + A, a,). (1.) 

Now this equation is satisfied by «,+;k = fl?,j which gives 
Aa«»0, and a^^^aj a constant, and leads to the ordinary 
or direct equation to the series, 

Ug- F (/Vy a). 

Also equation (l) will be divisible by o,+,^ — «#> be- 
cause of is a value of a^^^^ whkh satisfies equation (l); 
and if dimensions of a,^^ and a, superior to the first are 
involved in it, the result of this division will be an equation 
involving a^^^ and a^; i. e., an equation of differences of 
the first order with respect to a«, the solution of which will 
give one or more values of a, in terms of a? and arbitrary 
constants ; and these being substituted for a^ in the equation 

u^^F{wy «J, 

will furnish equations of series, which are primitive equa- 
tions of 

/(af, w,, A«*J = 0, 

and each involves an arbitrary constant. 
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If equation (l) does not involve higher dimensions of a^ 
and a^^^ than the first, they will disappear from the result 
when it is divided by a^^^- o^. In this case a, will have 
only one value, viz. a, = a, and there will be only one equar- 
tion of a series corresponding to the proposed equation of 
diflferences. The mode in which the indirect solutions just 
treated of are obtained, is analogous to that in which the 
singular solutions of differential equations are obtained; but 
whereas the latter can contain no arbitrary constant, indirect 
solutions of equations of differences may contain as many 
arbitrary constants as the complete integral itself from which 
they are deduced. 

Ex. w,»a?Aw, + F(Aw,). 
Taking the difference, we find 

Aw, = Aw, + (^ + 1) A^w, + AF(Aw,), 
or = (<r + 1) A*w, -I- Ai^ (Aw,), 
which is evidently satisfied by Aw, = a, a constant; 

.-. w, = aa? + a'; 
and substituting in the proposed equation 

aa? + a' = aa? + F(a% .-. a' = F (a) ; 
.*. u^^wa+ F(a)f 
the complete integral, containing one arbitrary constant. 

For the indirect solutions we shall have 
w,= a7a, + jP(a,), 
a, being determined from the equation 

(j? + 1) <x, + i + i^K+i) = (ci? + 1) a, + iP(a,). 
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Suppose, for instance, that F (a,) « a* 



5 

,8 



.-. (<2? + 1) (0,+, - a,) + oj+i - 0^= 0, 
or, rejecting the factor a^^i - a^, 

«*+i + «, = - (^ + 1) ; 

:t = - 2(^ + 1) (- 1)' = — — (- 1)' + C, 



or a,= + C(-iy; 

1 ^ A t^ 



Linear Equations of Differences of all Orders. 

76. The linear equation of Differences of the »*^ order 
is 

all the coefficients being functions of w ; the first step towards 
its integration is to establish the following theorem. 

77- If there be n particular values 
which, when substituted for w,, satisfy the equation 
that has no term independent of u,, its complete integral is 



U„ = Oj ^U„ + Og V, + ... + ««X> 



^19 ^9««*afi being arbitrary constants. 
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For let these values be substituted in the expression 

and it becomes, (collecting the terms multiplied by the factors 

«! C«*+« + Pi '«*•+-! + ... + 1>«X) + «2 (^^M+n + l>i*«*+«-l + — 

+ p^%,) + ... + a.C««*+« + yiX+-i + .-. +1>«X)- 

Now since ^«,, *w^, ...*«*,> satisfy the proposed equation, 
each of the quantities included within brackets is equal to 
zero, therefore the whole is identically zero ; consequently the 
assumed value of Ug satisfies the proposed equation, and it 
contains n arbitrary constants, therefore it is the complete 
integral of that equation. 

78. To integrate the equation of differences, 

«**+« + PiW.+— 1 + ft «*•+-« + .- + ft«* - Qj 
all the coefficients and q being constants. 

Assume u^^v^-^-k; then by substitution we get 

^*+Ji +Pl^*+«-l +ftV*+ii-2 + ••• 

+ p„u, + A?(i +ft-hft+ ..- +Pn) -g = o. 

Let k = , then the equation becomes 

1 +j)i+j)2+ ... +ft 

Let Vg - off then a* (a* + P\a*^^ + fta"""* + ... + p«) is the 
value of the first member ; now this will vanish if a be any 
root of the equation (called the auxiliary equation), 

/(a) = a» + p,a-» +P8a'-* + .•. + p«_,a +ft = 0. 
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Hence the (n) roots of this equation a^, 029 (hi***^* 
will give n particular values of «», a% aj, al...al which 
satisfy equation (l) ; therefore its complete integral is 

Vg = CjO* + CjdJ + ... + c«aj; 

and the complete integral of the proposed equation is 

Q 



u, = CiOi + Cgo; + ... + c^a* + 



l + Pi +P2+ ... +p« 



79- If the auxiliary equation have equal roots, the 
above ceases to be the form of the complete solution ; because, 
in that case, it does not involve the due number of arbitrary 
constants ; and it must be modified as follows. Suppose two 
roots a^, Og, to be very nearly equal to one another, so that 
'lis = 01 + A, h being a very small known quantity; 

then, 

Cia' + CgttJ = (ci + C2) aj + Cg {wa'^^h + -— a'-^h^ + &c. } 

1 • % 

replacing the constants 

Ci + Cji by Ci, and — - by C« ; 

now this continues true however small h be taken, and 
therefore when AsO, when it becomes 

(Ci + C^^pX; 
.-. u, = (C*! + CuOf) al + c^al + &c. 

Similarly, if the auxiliary equation have r roots equal 
to Oi, the complete solution will be 

w# = (co+Ci«r+C2a?*+... + c,-.i^r'"'*)af + c,+ia?+, + ...+c,a;; 
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of the correctness of which, we may be assured by the fol- 
lowing reverse process. Assume t*» =«'«,, then 

and the first side of the equation becomes, when divided 
by a% 

a" (1 + A)"iJ, + pi a"""* (1 + A)**^tJ, + &c. + p^v, 

= {(o + a A)" +pi (a + aA)""^ + &c. +p.} v, 

a*A* 
«/(a + aA)t?, = {/(a) +/(a) . a A +/'(a) — -+ &c,} t), 

=/(a) . t), + ^ ./ (o) . At?, + :;^/' (a) . A'«, + &c. + a«A"t?,- 

1 1.2 

Now suppose o = ©1, and /(a) = to have r roots equal 
to ax\ this makes the terms as far as /^'"^^(o) vanish; and if 

v, = Co + CiW + ... + c,.ia7^"V 

then A'v, = 0, A'"*^*!?, = 0, &c., 

and all the remaining terms vanish; and consequently the 
equation is satisfied by 

«, = (co + Cia? + ... + c,_ia?'-*) a\. 

80. Also if the auxiliary equation have a pair of ima- 
ginary roots, 

m ± » \/^ = p (cos 9 ± \/^ sin ©), 

/ 71 

putting p = V m* + «*, tan = — , 
the corresponding terms in the value of u„ will be 
Cp' (cos d + \/^ sin e)' + Cy (cos - \/^ sin By 
" /o' (ci COS wQ + Cg sin ^d), 
changing the arbitrary constants. 
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And if there be r pairs of imaginary roots, the corres- 
ponding terms in the value of Ug will be 

+ (*o + *i^ + ••• + K-i^"^) p' (cos - \/- 1 sin 0)' ; 
or, changing the arbitrary constants, 

(co + Cia? + .•. + Cr^iOt^'^) p' cos w0 
+ (c'o + c\af + ... + c'r^iaf'^)p'mkaf0. 

Obs. The Differential Calculus being a particular case 
of that of Finite Differences, a strict analogy exists between 
the methods and results in the two subjects, as the reader 
cannot fail to have observed; indeed whenever in the latter 
a result is obtained with an indeterminate increment for the 
principal variable, it is possible to pass to the corresponding 
result in the former, by a method similar to that pursued 
in Art. 25 ; and which we shall further illustrate by the fol- 
lowing instance. 

In the equation u^+^h " ^^^*+a + (**** + ^*) ^# ™ 0> 
putting Ug^afj we find a =s (i»±«\/-l)*; 

= (m* + n*)^* |cjCos f - tan"^ —J + Cg sin (- tan*^— j > . 

Now the proposed equation is 

A*w, - 2 (w - 1) Aw, + {(m - 1)* + «*} u, a 0, 

or if we replace the known quantities m - 1 and n by other 
known quantities m^h and n^hy the equation and its scdution 
take the forms 

A*«*, At«, , - ^ 
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t., = {(1 + m,hy + n\V}^^ |c, cos \^ tan-» IT^) 

+ Cg sm r tan"* ; > . 

Now take the limits of these expressions when A = 0, and 

we find the well known results, since (l+2WiA)^* becomes c*"s 

» 

f^^ = e"i' (Cj cos »j«r + Cg sin n^w). 

81. To integrate the linear equation of differences of 
the n^ order, 

the Goeffioients being functions of of. 

If ^Ug, ^w„ ... "f^, be w particular values of tr, in the 
equation 

«•+« + Pi«»+— 1 + i>2«*+«-2 + ... + p««, = 0, (1.) 

with which the proposed coincides when its second member 
is zero, we- have (Art. 77.) 

t?, = Ci^U, + Cg ^u„ + ... + c."t^,. 

If we now divide both sides by ^tt,, and take the difference, 
we shall eliminate e,; next dividing both sides by the co- 
efiicient of Cg, which suppose '«,, and taking the difference, 
we shall eliminate Cg; again dividing by ^«j„ the coeiBcient 
of C3, and taking the difference, we shall eliminate c,; and 
proceeding in this manner till all the constants are elimsnated, 
our fiffial result will be of the form (eadi A affieeting the 
whole of the expression that follows it) 
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in which expression the coefficienlb of v,^^ i» eyidendy 



^^*+l-" X+2-' '^x+n-l ^*+» 



therefore, dividing by this coefficient, we get the expression 

Vg Vg Vg U'g 

whidi must be equivalent to the first member of equation (l). 
Therefore the same expression, only with Ug instead of «,, 
must be equivalent to the first member of the proposed 
equation, and consequently equal to JT; hence, equating these 
equals, and integrating, we get 

Ug = 'Ug^'vg'S.'vg ... 2"-'t?, 2;jzi rr2 — 1 , 

(each 2 affecting the whole of the expression which follows it) 
which is a general formula for the integration of any linear 
equation of differences whatever. 

82. As a first exemplification of this method, suppose 
the coefficients of the equation to be constant, and a^ Og, ••• a« 
to be the n roots of its auxiliary equation. 

Then «, = CjO' + Cgag-h 03©;+ ... +€^^1, 

changing, both in this, and in the similar succeeding steps, 
the arbitrary constants; 

U3/ UJ <* \Os/ *" "Us/ 
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in which exprefiAsion the coefficient of v^^^ is evidently 



therefore, dividing by this coefficient, and replacing v^ by u^^ 
we get 

or, if we choose to introduce the arbitrary constant after 
each integration, 

i#, = Cj Oj + Cg aj + ♦.. 

S3. If the auxiliary equation contain equal, or imaginary 
roots, this method is still applicable ; it is only necessary to 
assume for t^^ the value belonging to the case of equal or 
imaginary roots, as will be seen in the following instances. 

Ex. 1. To solve the equation 

Here «, = (c + c^a) — ; 
a* 

.: A«(oX) = 0, 



* For this method, and the corresponding one of solving the general linear 
differencial equation, I am indebted to Mr Gaskin. 
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in which expression, the coefficient of v,^, is a*+"; therefore, 
dividing by this quantity, and replacing v, by u^^ we get 



^A''(o'«.) = ^; 



.-. «. = l2»(^a'+«). 



Let ^-«, then r'C^a') "'"S^*'- 22'o'+' + c'» + c, 

(Art. 58.) 

^ ^ if (a^iy {a -if 

Ex. 2. To solve the equation 

Here Vg^Cip AnwO -¥ t^p* cos «0, 

where ff^nf •\' »*, tan d = -- ; 

ft} 



-f-^^Cjtan^e + Ca; 



^cos^d 

... A ( ""' \ « ^^^°^ . 

Kp'isosafd) coswOcos^w-k- 1)0^ 

.% A cos ae cos (^+ 1) d A I ^ ^' ^ I :^ 0, 
'^ V/o'cosa?^/ 

in which expression, the coefficient of v^+g is \,+2 * 
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therefore, dividing by this coeiBcient, and replacing «, by u^^ 
we get 

— -T-z. A cos w0 cos (a? + 1) A I -i — ^^ I « -ST; 

cos (a? + 1)0 "^ ' \p' cos w9f 

{ X cos (af + 1)0] 
^+2 1» 

to which we may add the terms c,p'sin^0 + C2|o*cos^0, if 
we suppose a constant to be a;dded after each integration. 

Ex. 3. To solve the equation 

where the auxiliary equation has two pairs of imaginary roots, 
so that 

«* +J>i «' + P2«^ +1^« +P4 = {(« - mf-^-n^y. 
Here t?,= (o + a'a) p'cosaO + (6 + A' .a?) p* sin a? 6, 

••• ^^ — y^^^ + a' 0}+ (b + 6'd?)tana?0, 

p' cos ^ d 

• ^ 

.-. A —— ^^ = a +6'tan^0 + (5 + fe'^) ^ -r-^, 

p'cosxQ ^ coswdcos{w-\-l}9 

••. cosi»dcos(^+ l)dA ^^=Ao'^cos(2a?+l)0 + cos0| 

^ p'coswO * ^ ^ 

+ |6' {sin (2 a? + 1) - sin^} + (6 + 6'^) sind, 

.-. A cos ^ 0COB (a? + 1) A -r; — * — -z. = — a sinS (a? + 1) d sin d 

p cos w 9 

+ 6' COS 2 (/v + 1 ) sin d + b' sin 0, 
.-. sec* (a? + l)'0 A COS 07 cos (a? + 1) A 



p*COSd?d 

= - 2 o' tan (a? + l) d sin + 26' sin 0, 
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2 a' sin* 
^ ^ ^ ^ cos(ir + l)0cos(a? + 2)0 

. /. A cos {jjo + 1) cos (^ + 2) A sec* {w + 1) Q x 

A cos a? Q cos (a? + 1) d A ^ — - = 0, 

, . , . ,, ic • ^ i« . cos (4? + 2) 
in which expression, the coemcient 01 u,+4 is j^^ ; 

therefore, dividing by this coefficient, and replacing v, by f^„ 
we get an expression equivalent to the first member of the 
proposed equation ; hence, equating it to Xj and integrating, 
we find 

t^, = p' cos «r 2 sec «r sec (a? + l) 2 cos* (a? + 1) x 

;^4 }- 

84. Having given a particular integral of the equation 

to find its complete solution. 

Let v« be a particular integral, then 

Hence, eliminating ^«, we find, putting — = «,, 

A«,+i B^v„ 

— — - « « «I7„ suppose ; 

A Xg u^^-g 

or A log ( A «,) = log M?„ 

.-. Ai8r^=CPw,.i, (Art. 50.) 
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86. If we know a particular integral of a linear equation 
of any order that has no term independent of «««, we may 
reduce it to another equation of the same kind of the order 
immediately inferior. 

Let u, s V, be a particular integral of a linear equation 
of differences of the n*^ order reduced to the form 

A*«, + 5PjA"-^«, + 5P8A"-'w,+ ... +g.w,«/(A)«,T=o (i). 

Assume u, « v^Sti;^ ; 

then A"w.« {A + (l + A)A'}"«,2ti;,, 

where A affects v^ only, and A' affects ^w^, only, (Art. 34.), 
and the proposed equation becomes 

/{A + (1 + A) A'} «,2ti;, =0. 

Now /{A + (l + A) A'} is a rational integral function 
of A and A\ which we wish to arrange according to powers 
of a'; and this may be done at once by Taylor's theorem, 
which gives 

/(A) + A7i(A)(1+A) + ^V,(A)(1 + A)» + ,-+A'-(i-hA)«. 

Hence, observing that /(A) «, = 0, since t>, substituted 
for u^ satisfies the equation, and that (l + A)^o« » v,^.,, we 
get the depressed equation 

w,/i(A)«,+i+-— Aw,/g(A)«,+,+ ... + A*-' «;,.»,+,= 0; 
1 .ss 

or, reversing the order of the terms, (since /, (A) means the 
same function of A, that d^f(a) does of a), 

«'.+«A«-*w, + («» A«,^...i + ?iV,^..i) A-^fi;, 
+ /i(A)«,+i.w, = 0, (2.) 
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a linear equation of the (n - l)*** order, of the same form 
as the original one. Similarly, if we know another particular 

value of w,, ^^j then A (^ j will be a value of w^ in equa- 
tion (2), which may be depressed to another of the same 
form of the (t* — 2)*^ order; and if we know r particular 
solutions of equation '.(l), we may in this way depress it to 
an equation of the same form of the (n — r)*^ order. As 
a linear equation of differences of the first order and degree 
can always be solved, it appears that to obtain the complete 
integral of a linear equation of the n^^ order, we must know 
7i — 1 particular solutions. 

86. Besides linear equations of differences, and such 
equations as can be reduced to that form, very little is 
known of equations of differences of the second and higher 
orders. The following are instances of equations which admit 
of reduction to linear equations with constant coefficients. 

where pi, P29 Sec. are constants, and v, is any function of 
w. Assume u, « w^Pv^^ then the equation becomes divisible 
by Pi?,+„, and is reduced to the linear equation with constant 
coefficients, 

2- «*,+« + pia'«,+._, -^Pia^Ux+n-i + ... + p^a^'u^ = 0. 

Assume ««, = u^a^^"^^"^* "^'^ then any term PrO^^'u^+n-r 
becomes 

-« « ^**(.*+2« + l) i(«-r)(n-r + l). 

therefore the equation becomes divisible by a**^'^**"''*^ and is 
reduced to the linear equation with constant coefficients, 

9 
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^' w,+s + {a + fc(-l)'} w,+i +'cw, « 0. 
Let u^ = v» y/a^ 6(-l)', 
then w,+, = v,^., \/a - fc (- 1)% 
w*+2 « t?*+« \/a + 6(-l)' ; 
••• ^M+9 + V a'* - 6* 15,+, + cv, = 0. 

4. w,+i «, + ««,+! + 6«^, + c = 0. 

Assume w, + o ■= -^ , then t«.+| + a « -^^ , 

''#4.1 /^#4.2 ^ . , (^x^l \ 

or t?,4.8 - (o - 6) t?,+i + (c - ab) v, = 0. 

The two arbitrary constants which will appear in the 
value of Vgn must be reduced to a single constant, by the 
condition of the proposed equation being satisfied. 

5. f/,+, - 2wJ+ 1 =0. 

Assume r/, « c^s v,^ then n, = ^ (c** + c~^). 



Prob. 1. To find the value of 

V S - \/2 - \/2 - ... to .T terms, 
let it s= «,, then f/,+, = v 2 - «, ; 
.-. w* + i + w, "=2; let w, as 2 cost?,, 

then 4 cos* r^^, + 2 cost?, « 2 ; 
.-. 2 cos*fj,+, = 1 - cos V, =• 2 sin*-|«,; 

*v,+i « sin ^t?, - cos f~ - iO; 



COS! 
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.% w, = 2 cos |c (- ^)' + ^1 ; 

-- + -]= V 2 = 2 cos - > .-. c = ~ , 
.-. u. = 2cos{^(-i)' + j}«2sin^{l«(-i)'i. 

Peob. 2. A swan breeds three cygnets in its second year, 
and'^four every succeeding year; ani the young ones all 
breed according to the same law; required the number at 
the end of the w^^ year. 

Let u^ = number at end of o?*^ year ; 

then in (a? + 1)*^ year, there will be bred .three each by those 
in their second year, and four each by all the others ; 



.% w' - w'^ - 3m - 1 = w (w»^ - 1) - (w + 1)' = ; 
... m=-l, and l±\/2 5 
... w, = a(-l)' + 6(l+\/2r + c(l-^/2)^ 
t^^ = = a + 6 + c, 

t^i = l=:-.a + 6(l+ \/2) + c (1 - V'^), 
i^jj=:4 = a + 6(3+ 2\/2) + c (3 - 2 \/2) ; 



''s-2i 
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'■■ «,=(-ir+i(-i=-i)(i+-/ir-i(-^+i)(i-v'i)'. 

If the law of increase be only one the second, and every 
succeeding year, 

87. The following problems relative to continued frac- 
tions, give rise to equations of differences, which can be in- 
tegrated as linear equations. 

1. To determine the value of the continued fraction 

c c c 

u, = to ^ fractional terms*. 

o + a + ... a 

.-. w,+i = — - , or «,+, (tt, + a) « c ; 
a-tu. 



assume 



Ug + a^ , then —^{ a] 



or v,+2 - « 'f^x+i - c«, = ; 
«, '=^ CiQ* -^ Cg /3*, a and /3 being roots of Ap* — oA?-e = 0; 






'"' a-^C3/y- - putting c,= - 
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C l+Cg /3 

a a + Cs)3 a 



If a and )3 be imaginary, so that 

Ji^ — ak "C^k^ — 2k p cos + p% 

, sin,r0 
then w, = - p -7-— — - . 

2. To find the value of the continued fraction 
c c c c 

c 
to w fractional terms, the last term being -« 

Proceeding as in the last example, vre find 



tt^ = c 



a' + C3/3* ' 



- c 1 + Cg 6 — a 

3. To find the value of the continued fraction to w 

fractional terms, 

1 J 1 1 ^ 

w, = £ I — 

„ 1^ 1 1 

Here w,+2 = 9 .% — a = 



1 w,+j| 6 + u* 

a + 7 

6 + w, 

or a o + w,. 

1 - aw,+8 
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Assume !-««,« , .-. 1 - a w^+a = > 

or t?,+3 - (2 + o 6) v,^i + tJ^.! = 0, 
of which the auxiliary equation is 
k'^(2 + ab)k^-^ 1 =.(A;*-\/^A;-l)(A;* + \/a6ft-l) = 0. 

Let a 9 , - a, - , be the roots of this equation, then 

a a 

«x = «' {^1 + <^ (- 0'} -i- ~ 1^3 + ^4 (- 1)'}, 

a 

where p..|±^(^. « C. ^(- ^r^; 

aMg= 7« 1 5 ;., 

.\ a6 + 1 »a'*- 1 +o~* = 7-, ••• a — «va6; 

a + ^2 a" « 

•\ Pa = 1 « C - C, .-. C = 0, C = - 1, and p, « (- 1)', 



O tt- s= 1 ^ 



a'+^ + (-l)*a"-i- 
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4. Suppose 6 to be essentially negative and « — c, so that 
the continued fraction is 



and let y/ah = \/^\ y/ac = 2\/- 1 sin 0, 

then ft^ - 2\/^ sin 0Ap -- 1 = 0, 
.'. a = cos + v-l sin ; 



sin (a?- 1)0 
sm (^ + 1)0 

according as «r is odd, or even. 



or 1 - 



cos (<r - l) 
cos(a?+ 1)0' 
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SECTION IV. 

SUMKATIOK OF SERIES. 



Integration of the General Term. 

88. One of the most direct and important applications 
of the Calculus of Finite Differences, is the general method 
which it furnishes of assigning the sum of any number of 
Terms of a Series of the general term of which we are 
able to take the integral. There will be two cases to con- 
sider, according as the general term is given explicitly in 
terms of the index, or is only given by means of an equation 
of differences. We shall begin with the former, by shewing 
that the sum of any number of terms ending with the 
general term u^^ is equal to the integral of the following 
term, plus a constant. 

Let Sg denote the sum of the w first terms of a series 
whose general term is u^^ 

then aS, = «i + Wg + ... + 1^,» 

Making ;r = 0, Sq'= = ^*U:,^i + C ; 
or, as it Js usually written, 

« 2^1 + C, 
.-. *y, = 2w,+i - 2«i. 

In general the arbitrary constant will be determined by 
the term with which we make the series commence. Thus, 
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if we make it begin with Ur instead of Ui, this amounts to 
suppose the sum of the terms preceding w^ that is, 2Wr to 
be zero; and therefore to determine the constant, we have 

*S'r_i = 2w,=, + C = 0. 

89. As the summation of series, where the general 
term is given explicitly as a function of the index, thus 
resolves itself into the cases of integration treated of in 
Section II., it will only be necessary to give a few numerical 
examples; every expression integrated in that Section gives 
the sum of a series of which it is the general term. 

1. To find the sum of the first w terms of any pro- , 
gression of figurate numbers. In the r^ order the general 
term is 

_^w(af + l)(w+ 2) ... (a? + r - 2) 
^' 1.2.3...(r-l) ' 

. . 2.u^+i = -— — + C, (Art. 45.) 

1 ,X * 3 ... T 

2^1 = + C =0 ; 

w(af + l)...(a? +r- 1) 



^- = 



1.2.3...r 



Similarly, for the series of inverse figurate numbers, ex- 
cept when r = 2, that is, for the series 1 + ^ + ^ + ^ + ... + -. 

^ . 1.2.3...(r-l) 

For we have «^, = - • 



a?(a?+ l)...(a?-Hr- 2)' 

^ ^ 1 .2.3...(r- 1) 

2w,+i = C - 7 -7 ^ ^ ^ ^ , (Art. 47.) 

(r-2)(a? + l) ...(^ + r-2) ' "^ ^ 

2wi = C - --^^ = 0, 
r -2 

10 
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{ I 1.2.3...(r-3) 1 

*^ ' " ^"^ " ' ^ t^^^Ti " (a; + 1 ) (^ + 2) . . . (^ + r - 2) r 

and iS = . 

« r-2 

2. To find the sums of the squares and cubes of the 
natural numbers, 



u,^ a^ ^ (of ^ 1) w + w 



.-. 2t.,,, = ^ +— ^ +C, . 

2«i = + C « ; 

^(a?+l)(2ay-H) 

Again, «^^ « a;* = a? (a?^ - 1) + ^ ; 

(a?-l)a?(a? + l)(ar + 2) ^(a^ + 1) 

••• 2«'+i 1 + — i +^' 

Similarly, l» + S^ + 5^^ + ... + (2a? - 1^ = -^^— . 

l3 + 38 + 5' + ... + (207 - ly = 2a?* - /P*. 

3. To sum the series 

l«-2« + 3^-4«+... ±a?*; 
«,=:(-l)-V; 

2..J-0'---^(^--^^)(-^)^^('^)"^C (Art.59.) 

-2 4 -8 
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Similarly^ 

4. To find the sum of w terms of the series 
1 1 1 o 



2«-s» 4»-S« 6«-3« 

*** * 7 — Cii — i = 7 TT ; J which falls under Art. 48. 

(2a?y-S* (2a? -3) (2^ +3) 

1 jy-H 1 GOBI'S 

'""is "" (2a; + 1) (2a? + S) " 6 ' (2a? - 1) (2a? + 1) (2a? + 3) ' 



and Sm^ — . 
18 



Similarly, 



5 8 11 3a? +4 

- + + + &c. to w terms = 2 - 



1.2.3 2.3.4 3.4.5 ' (a? + 1) (a? + 2) * 

5. 2. 4. 7 + 4. 7. 13 + 8. 13. 25 + 16.25. 49 + &C. to a? terms, 
w, = 2* (3 . 2*-' + 1) (3 . 2* + 1). 
Assume 2w,+i = -4(3.2'-* + l)(3.2'+l)(3.2'+^+l); (Art. 55.) 

.-. tt,+,«^(3.2' + l)(3.2'+' + l){3.2'+*+l-(3.2-^ + l)} 

21 
^ — A. 2'+* (3 . 2' + 1) (3 . 2'+> + 1) ; 
4 

'• "* 21' 
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2tt,+, = ^(S.2'-»+ 1)(S.2'+ l)(S.2'+'+ 1) + C; 



21 2 

4 40 
.-. *S,«— (3.2-' + 1)(3.2'+ l)(S.2'+* + l) . 



2 4 8 16 , 

6. + + &c. 

3.3 S.g 9.15 15.33 

«. ^-^ 



16 21 2 26 /2\* . 

7. . - H . - — &c. to X tenns, 

2.3.4 S.4.5S *.5.6 \SJ 

u ^6 + 5(0,-1) 

' (a;+l)(a7 + 2)(a^+3) ^ ^^ * 

Assume 2«,,. = (^^g)^^^^^) (-f)'' (A^. 56.) 

then u4 = - 3, and 
3 



*^'''i"(^ + 2)(-P + i)^"'t)'- 



Similarly, 



19 1 28 1 39 1 52 1 

7 + 7—:: — : • :: + :: — : — z- 7^ + z — ^-7; • — + 



1.2.3*4 2.3.4*8 3.4.5* I6 4.5.6' 32 
to w terms, 
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X {po + 1) (^ + 2) 



(ir'> 



8 
to w terms. 



1 1 1 . 

COS0COS20 COS20COS30 cos3dcos40 



«« = 



cos ^ d cos (a? +1)0' 



_ tan (a? + 1) ^ , ^ 

2w,+, = \ ^ + C, (Art. 53). 

sm0 ^ ^ 

sind 

tan (a? + 1) - tan Q 
smd 

9. 1 cos + 2 cos 2 + 3 cos 3 + &c., to w terms. 
Ug^i = (^ + 1) cos (a? + 1) 0, 

2t^,+i = (0^+ i)!ii^^lil^ - _!_2sin(a? + |)0 (Art. 59.) 

2 sin - 2 sin - 

2 2 

_^ (g? + 1) sin {w + ^) g cos (<g + I) g 

2 sm - I 2 sin - 

2 V 2/ 

2«,= J +- 5--,+ C = 0, 



(2 sin-) 



1 



(2 sin-) 
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„ «8iii(«+i)0 1 f. . iv/, . ^ -j^+^fll 
•• *,- g* + 5J8in(« + i)e8m--sra*-^0| 



2 sin- 28in'- 

2 2 



2 sin - \ sin - / 



2 2 

Hence also if *?,= lcos*0 + 2cos^20 + 3cos*30+ &c. 
then 2^, = l(l +COS20) +2(1 + cos40) + &c. 



4 4sind 



/sina?0\* 



10. tan-^ = + tan-* i + tan** + &c. 

1 + 1 + 1* 1 + 2+2* 1+3 + 3* 

to w terms. 

1 A* 

Ug = tan-* r = tan"^ -r , 

1 +a + ar I -{-w^w + 1) 

.-. 2 w, = tan-^ a? + C, (Art. 54.) 

.-. *?, = 2w,+, -2tti = tan-*(a?+ l) --. 

4 



Recurring Series. 

90. We next come to the case where the general term 
is not given explicitly in terms of its index, but only certain 
relations between the consecutive terms, or these and their 
indices are expressed ; of this sort of series the most remark- 
able are Recurring Series. 

A recurring series is a series in which an equation of the 
first degree with constant coefiicients, holds good between a 
Certain definite number of consecutive terms, in whatever part 
of the series they be taken. 
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For example, in the series 

3 + 5 + 9 + 17 + 33 + &c. 
we have 9 = 3.5 - 2.3, 17 = 3.9 - 2.5, &c.; and in general 
w*+2 = 3w,+i-2w,. 
The general equation of every recurring series is 

The series of coefficients which connects any term with 
the preceding ones is called the Scale of Relation. Thus 

1 +i>l+l>2+ ••• +Pn = 

is the scale of relation of the recurring series whose equation is 
«*,+« + J>iW#+«-i+ ... +l>»t^,= 0. (1) 

91. A recurring series may generally be resolved into 
two or more geometric progressions. 

For if ai, »£, Og, ... a^ be the roots of the equation 

o» + p^ a»-^ + P2 a*'^ + ... + p„ = 0, (2) 

the complete integral of the equation of the series is 
Ug = Cj a' + Cg aj + ... + c„ a^, (Art. 78.) 

Wg = Ci ttj + C2 fltg + C3 ag + ... + c„ a^. 



«*« = Ciai+C2o; + C3a;+ ... +Cnal^ 
and the series consequently is transformed into 

c, («! + aj + ... + a*) + (?8 (og + oj + ... + oj) + &c. 

+ c«(a,+ a* + ... +<). 
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92. In the particular cases in which equation (2) has 
equal, or impossible roots, the recurring series can no longer 
be resolved into geometric progressions ; for the complete in- 
tegral of the equation of the series becomes in those two 
cases, respectively, 

w, = (cq + CiW +C2af^ -¥ ... +c^«ici?*'~*)a'-h c^+ia;+i+ ...+c„a*, 

W, = (Co + Ci a? + ... + Cr-i (Jf-^) p' COS w 9 

+ (cj + c}^ + ... + cl^iOf'^) p* mi Of 6 + C^r+itttr+i + •••+C»««* 

93. Hence, to find the general term of a recurring series, 
we must integrate the equation expressing the relation be- 
tween its successive terms, and determine the arbitrary con- 
stants by making the general term Ug coincide with a sufficient 
number of given terms. When a series is known to be 
recurring, its equation may be determined by assuming it to 
be of the form (l); and then forming a sufficient number 
of equations for finding the coefficients p^, p^i &c. Pn^ by 
substituing the given terms in order. 

94. To find the sum of a? terms of a recurring series. 
First, let its general term be of the form 

Ug s c^a* + Cgaf + ... + CnKi 

»! — 1 ^2 ~ * 



Secondly, let the general term be 
^s = (co + CiW + c^o^ + ... + Cy-iO?"-^) a\ + c,+ia'+i + &c. ; 
then since 

2af a' - ^. — + -7 -r^ - &c. (Art. 59.) 

«! - 1 1 «! - 1 (ttl - U J 
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.•. iS, sa 2w,+i — 2«i, is determined. 
Thirdly, let ' 
u, = p" (co cos w9 4- Cq sin wO) + p*i€ (Cj cos mQ + c^ cos ^0) 
+ p'cf (cg cos /p0 + c/ sin a;0) + &c., 
then each term to be integrated will be of the form 

afp* cos {wB + a), 
the integral of which may be found by Art. 58, because 

2"{p'cos(^0 + a)} 
is always assignable. 

Ex. 1. To find the sum of w terms of the series 

1 + 5 + 17 + 53 + &c. 

Let the equation be w^^+a +p^«+i + 9«** = 0; 

.*. i7 + 5p + g = 0, 5S + 17y + 5gs=0, 

which give j) = — 4, g = 3 ; 

••• «^*+2 - 4w,+i + 3w, =ss 0. 

Let UgM^a*\ .-. a^- 4a +3 = 0; .-. a^S or 1; 

.•. w, = Ci3* + C2> 

1 = 3Ci + Cg, 5 « 9Ci + C2, 

.•. Ci = f , Cg - - 1, 

.-. w,=*2.3'-*- 1; 

3' 

.-. ^, = 2. -«a?+C, 0=1 + C; 

3 "" 1 

Ex. 2. 2-o-a' + 2a5-fl^-a* + 2o®-d'-&c. 

, (2a» - 2) tt 

w, = 2a*~ cos -^^ — . 

3 

11 
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Ex. 3. 1 + 4 + 18 + 80 + 356 -f &€.» to w terms. 
Here m,+2 - *«^*+i - 2 w, = ; 

.-. 2\/6 w, = (2 + \/6)' - (2 - \/6)% 
and 

^ /^o (2 + \/6r^-(2 + v/6) (2-^/6r^-(2-^/6) 

1 + y/6 1 - V 6 



Application of the Integral Calculus to the Summation of Series. 

96. As integrals are often expressed by Series, so, 
conversely, the latter may be represented by integrals; and 
it is often desirable to find the integral of which a proposed 
series is one of the developments, in order to subject it to 
the methods which we possess for calculating, at least ap- 
proximately, the value of any integral taken between assigned 
limits. We proceed therefore to notice one or two processes 
given by Euler for effecting this; they consist chiefly in 
performing certain operations on the series, by which it is 
transformed into another series which we are able to sum, 
or which is similar to the proposed one. 

96. Series which proceed according to the powers of 
some quantity ^, affected with coefficients consisting of factors 
in arithmetic progression either in the numerator or deno- 
minator, may be summed by the aid of the Integral Calculus, 
the denominators being taken away by differentiation and the 
numerators by integration. 

t 2^ 3^ 
Ex. I. Let « « - + — + — + &c. (to 00 ), 
2 3 4s 

.-. dt(fi0 = ^ + 2^ + 3^+ &c. 



-x 



V 1 ~— f 
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d,(st) 1 

= + 1 + C; 

log(l-0 . 1 



.*. 8 = 



1 -t 



Ex. 2. « = t + -r f + . 

a + )3 a + 2/3 



a + nb 

+ TT f» + ... to m terms. 

a + nfi 



Ex. 3. « = 1+ ^t+- 1, \as ^ 

a+fi (a + /3) (o + 2/3) 



fe/«*^ -a-6'^a + /3'^(a + /3)(a + 2/3) 



a (a + 6) j+2 



(a + /3) (a + 2/8) (a + 3/3) 



t" + &c. ; 
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(a + jS) (a + 2j8) 



a - 6 

hence, performing the differentiation, we find 
{jit - hf)di8 + {a-at)8^ a, 
a linear equation of the first order for finding «. 

For the sum of m terms, we shall evidently have from 
equation (1.) 



a a 



-r+1 i. l-K O --1 



^ /i L _ a(a + b)...{a + (m-Z)b} 1 

17 1.2 fy(fy+l) 

a(a + l)(a4.2)^(/3 + l)(/3 + 2) 

+ , ^ ^ 7 TTT TT ^ + &c. (to 00). 

1.2.3 7(7+1X7 + 2) ^ 

(fi'-t)cP,8+{(a + ft-^l)t-y]dtS + afi8^0. 

97* When consecutive denominators have only a single 
factor in common,* or none at all, the following is a conve- 
nient mode of summing the series. 



Ex. 



a (a + 6) (a + 26) (a + 2b) (a + 36) (a + 4fc) 

+ &c. (to 00 ). 
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g.^^^ sy 1_ 

(a + ofb) {a + (^ + 1) 6} {o + (^ + 2) ft} a + wb 

2 1 



a + (a? + 1)6 a + (a? + 2)6' 

resolving all the terms by means of this formula, by putting 
zp = 0, 1, 2, &c., we find 

a a + b a-\-2b a + 3b o + 46 

.-. 26»^=.2 / .--. 

*^g I +ar a 

™^ ^^' a (a + 6) (a + 26) (a + 36) 

"*■ (a + 3 6) (a + 4 6) (a + 5 6) (a + 6 6) "^*'^' ^*^ ^ ^* 

"e^U"^^. l + w' + a^^l' 

98. We may obtain an expression for the sum of the 
series, 

1 + ^ d,<f>(z) + ^ <^ {<P(Z)]» + j^^ «f. |0(;.)}» + &C., 

where («) is any function of x, by Lagrange's Theorem. 
For i£ yx + axp (y), then 

y - » + * (») + — rf. {0 (»)}» + &c. 

where dy<p{y) must be obtained in terms of x and «, from 
the equation y^siS'\'W<p{y). Suppose, for example, 

0(y) = i(y*-l). then dy0(jf)-y. 
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and y^x + -(f/^-l)f .-. l -ary = y/l ^2xaf + a^; 



% 



1^ '^ ^ 1.2 2» '^ ^ V'l-g^a^ + a?* 

Hence if (l -2ar^p + a?^)"i=i +Zij?+Z2a^+ ••. +Z„a?" + ..., 

then Z„ = ~V-^(^-l)*- 

The function Z„ has some remarkable properties, of which 
the following may be noticed. 

From the formula 
r /. . , . . V = -7=r log |\/a + 6^+ V-(c + ca?)[ 

it is easily shewn that 

x'c.-,„....^-..(.-^«.^) -'.1.0,(11-:), 

.-. "y^' (1 + Z^ar-i.Z^a^T'^- ...) . (l + Zi% Z^^H- ...^ 

2 / «' a* , \ • 
==_ a +--+-- + &c. . 
a V 3 5 / 

Hence the definite integral of every term of the product 
forming the first member that is not independent of r, must 
equal zero ; 

... -;(^'(Z.Z.) = 0, and -;^*'(Z.)««^. 

99* To find an expression for the sum of the reciprocals 
of the rH^ powers of the values of y^ in the equation 

«^-y + 0(y) = 0; 

where 0(^) denotes any function of y. 
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Let «-y + 0(y) = C(a,-y)(a2-y) ... (a„-y), 

•• ~r~rT • H ... H . yi*) 

ss-y + <p{y) a^-y a^ - y a^-y 

Now developing the two members of this equation in 
powers of y, the general term of the second member is 

and by Taylor's theorem the general term of the first member, 
putting i-dy0(y)=/(y), is 



fMlMdi 



M«^)- 



But if /(y) = Jo + -4i»+ •.. +-4«y"+ ... , 

smce = -+-. + ... 4- -— ^ + ... , 

z-y % sr sn^^ 

, f(y) ^ Ao + Az+ ... +Jn^ ^^ . 

•• = ... + — T y^+-..j 

% - y ss^^^ 

therefore the general term of - — - is — — y*, whatever be 

% - y »*■*■* 

the form of f{%)^ if "^rf be restricted to such terms as in- 

volve negative powers of % only ; consequently the general 

term of Ml^ ;. -^^^^.y, and therefore of 

z-y s^^^ 

hence equating coefficients of y" on both sides of equation (1), 
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Now putting for f(x) its value, the general term may be 
resolved into 

[^^^ ["^^1 " \l^' [ ^' I' 

r-i *^ I ^^ — J"~IF ^^^^ 

of which the first and last will be destroyed by the corres- 
ponding parts in the preceding and succeeding terms ; there- 
fore, changing n into n — 1, 

But the second member of this equation is what Lagrange^s 
Theorem gives for the development of y~% in the equation 
^-y + 0(y)=O; therefore «_« is equal to the sum of those 
terms of the development of y"" which involve negative powers 
ot z. 

100. By the help of the preceding proposition it may 
be shewn that Lagrange^s Theorem, applied to the solution' 
of the equation 

gives an approximation to the least root. By what precedes, 
we have 



«_,_, 1 n + r ^(x) n +r 






where each series is restricted to those terms which involve 
negative' powers of «; and the number of those terms in- 
creases with (w), and if w be supposed very great, then each 
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series may be taken ad ififinitum. Bat whatever be the value 
of w, if both series go on ad infinitum, the value of the 
second member is 

Xl=y^.sr' + r*-»^(sr)+^d.{*'->0(*)}+&c., 

as given by Lagrange^s Theorem, and which may be verified 
by actual multiplication; therefore when n is infinite^ 

limit of -^^ « y% 

(tf denoting that function of x ithicb is given for the develtfpe^- 
ment of y^ by Lagrange^s Theorem). 

But if di be the least of the roots Si, Og, Os, &c.) we 
have also (Theory of Algebraical Equations, Art* 156.) 

8 

limit of — ^ = d[ , when n is infinite ; 



.•. if=za\i and y=tti^ the leftit root* 

lOl. To find the sunns of the series 
sin a S sin da 3 sin 3a 



kcosa ApcosSa A;cos^a 



-^ kc. (to eb>i 



1« + jfc2 • 2« + A;'^ S^ + k" 



+ -^-riT+ &^- Oo od), 



* i - i. I «« ^ sin na . ^- k ctfe na 
or the values of ^S^ —, — —y 'ST —^ — — » 

Firsts we have 

(l-»*>«"^* L 

-1 — ^-i -«e-'^*+2(pcos^ + p*Cos2cJ?+f)'cos3a?+&c.)^"K 

l-SOCOStT+p' ^ 



•SpcostT+p' 
12 
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and integrating both sides from-.vsO to xsta^^ and putting 
p = l, we get 

^ = ^ -H 2c-*** '5^ -r r- - Se-'^ ^^* — r- 

k kr + n^ kr-^-n^ 

k 



+ 2'sr 



k' + n^ 



Again, integrating from ^ = to orB^ir, and putting 
pal, we find 



1 - e-«^ 



^(c-^'^+l)* + 2^*y-- .-2e-*fe-^5«T3 i- 

k kr-{-nr 1r-\-n* 

Now subtract the foregoing result from this, and call 
the two sums we are in search of, Si and 82^ then 

and changing the sign of k^ 
-Stire'^^ 1 

therefore, adding and subtracting, 

1 nsinna tt c*^'—**) - c'-*('-"«> 



iS« = ^iS~ 



Apcosna ir e^(''-'^^ + e"^^''"'^^ 1 
k^ + n^ ^2k ^^e-^ 2tf»' 



These formulas were first given by Poisson, and ma^ be 
considered as embracing the chief results which have hitherto 
been obtained relating to the summation of series of the 
sines and cosines of multiple angles. 
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102. The definite integrals required in the preceding 
investigation may be found as follows. Suppose 

(l-p*)e-*' 
1 - 2p cos w -^p^^ 

to be the equation to a curve; and let it be proposed to 

find the limiting value of its area, from a: = to cP s a, on 

the supposition that p approaches continually to unity. First, 

we observe that all the ordinates will be ultimately evanescent 

except those corresponding to d? = 0, Stt, 47r, &c., which will 

1 4.© 
be very large, because for those values becomes a factor 

of the expression for y. We will begin by supposing that 
o<27r; and therefore the only ordinates that we are con- 
cerned with are those immediately succeeding that through 
the origin ; so that, making p s 1 - p, and then supposing 
a? very small, we get successively 

y = — Ll Li ^ ^ , ultimately ; 

p« + 4(l.p)sin«J P^^ 



.-. 7;y=:2tan--, 
P 

Tyj»=i = '^j making p«0. 

Next, suppose that the area is to be found from <r«0 
toa? = 27r; then, besides the area just found, there will be 
another portion immediately preceding the point for which 
0? = 2ir ; to find this latter portion, put a? = 2ir - or', and call 
the ordinate y; 

then ^^ (^-P')^'^""^^^' =e-«fc^y, 
^ I-2pcosa?+p* 

therefore, the portion of the area immediately preceding the 
second limit = -ttc"*'^''; 
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103. As the following important theorems in Definite 
Integrals are based upon the preceding in^restigations, we 
shall here give the proofs of them, though not strictly be- 
longing to the present subject. 

"^—^^^ ^ -/(/r)H.2^5'«p-cos»(c-^) -./(o^); 



IT V ^ ^ M. O 

1 - 2p cos (c- J?) -- + J>* 

CL 

.•• "T'y,=. -"•//'/(*) + ^'S- -7/-co8«(c -w)-.f{ai). 

CL 

Now, suppose c to be less than a, then for all values 
of X between «f » — a and ^ = + a, y is evanescent, except 
when a? t= c ; if therefore we write /» = c + ^, supposing z ex- 
ceedingly small, and then integrate with respect to ar, from 
« = to jsr » «y any small finite value, and double the result, 
and make p = l, we shall obtain the value of ""X^"yj^=x. 
Put making j) » 1 - p, we have 



Zafip) pair 

~V (pay+iirgy'" 



.'• Ty 


IT 


tan" 


pa 


•• Ty,..- 


«/(c): 







Wio) - "7. /(*) + « 's^ "7, <»s « («-'«') - •/(*)» (»•) 



^a a a 



S'TT Sir 

+ cos(c-a^) — ./(ar) + cos(c-a?) — /(a?) +...?; 
a a 
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therefore, making a infinite, 

/(c)=i limit of - ""jr+'IcosCc-a-) -+co8(c-.v)— +&c.}/(^), 
TT o a a 

or, /(c) = - 7," ~'f,*'cos{c-m)z.fiw), (Integ. Cal. Art. 116), 

TT 

a theorem given by Fourier, and included, as we see, in the 
theorem (l), which was first given by Poisson, 

104. When in physical questions a definite integral arises 
whose value cannot be exhibited in finite terms, or in a form 
convenient for numerical calculation, the method of Quadra- 
tures is used as a substitute for Integration. This method 
consists in taking a series of values of the function to be in- 
tegrated, multiplying these by the differences of the cor- 
responding vQ^lues of the independent variable, and adding 
together all the results. The sum of such results approxi- 
mates to the value of the definite integral, as the intervals of 
the independent variable are diminished. The method is 
equivalent to adopting, for the area of a curve which a 
definite integral gives, the sum of the areas of a system of 
inscribed or circumscribed polygons. 

Let the equal intervals into which the independent variable 
w is divided be taken equal to unity, and let /(^) be the 
function of a? to be integrated between assigned limits a and b. 
Take a value of it /(a + w-^) at the middle of the «*^ 
of the intervals into which 6 - o is divided. If the intervals 
be very small, /(^) may be considered constant from 

«=a + n-l to .r = o + », or A'/W 

may be neglected; and the value of the definite integral is 
^f(a + n - ^), where n is to receive in succession the values 

1, 2, a, ... 6 -o. 
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For greater accuracy , suppose A* /{a) to be constant from 

x-a to w=^by or A'/C^) 
to be neglected. Then 

and the integral for the interval in question is 

fx/ia + n-^ + x), from »= - |^ to «r = |^, 

The integral, consequently, from a to 6^ is 

n receiving the same successive values as in the former case. 
In calculating the value of a definite integral by this method, 
the intervals are to be taken smaller as the variation of the 
function is more rapid. 

Convergency and Divergency of Series. 

105. A series 

Wj + ^2 + Ws + ... -h w„ + ... ( to 00 ) 

is called convergent, if the sum «„ of any number n of its 
terms approaches continually to a finite quantity 8 as its 
limit, when n is indefinitely increased ; and divergent in the 
contrary case. 

106. When a series is convergent, the sum of any 
number of consecutive terms after the n^ continually tends 
to zero as n increases. For 
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therefore, as n increases, the value of the first member con- 
tinually approaches to s - 8^ or zero, as its limit. This 
being true when m » 1, we see also that w„+i, or the general 
term u^^ continually tends to zero as n increases, or each term 
is greater than the following ; but this, although a necessary 
condition, is not sufficient to insure the convergency of a 
'series. 

Thus in the series 

111 1 

12 3 n 

fi^„ B — continually approaches to zero as n increases ; but 

1 1 1 



««+«-««= 7+ r+ .-. + 



n + l » + 2 n + m 

is evidently greater than 
m 



n + m 



> ^, if m = n ; 



and therefore the sum of. » consecutive terms after the n*, 
does not diminish indefinitely as n increases ; consequently, the 
series in question is divergent. 

107. In the geometrical progression 

a + aw + aw + ... 4- aa?*~* + ... (l.) 
1 -^ 



«. = a. 



1-^' 



1 -af+" 

«•+!« = «•—; — — > 
1 — /p 



1 -07" 



1 -ar 
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Hence if ^< l, when n is infinite a^mo, 

a 



8 = 



1 -J? 

also w«=0, ««+,-«« = 0; 

both which results shew that the series is convergent. 

But if d? > 1, then u^- aaf"^ increases indefinitely with 
n, which alone shews that the series is divergent. 

Hence the geometrical progression (l) is convergent or 
divergent, according as w is less or greater than unity ; and 
it may be used as the test of the convergency or divergency 
of other series. For if a proposed series can be shewn to 
have no term greater than the corresponding term of (l) 
when <r < 1, then that series is convergent; or if a proposed 
series can be shewn to have no term less than the correspond- 
ing term of (1) when w> l^ then that series is divergent. 

Thus in the series for e^ the base of the natural loga- 

rithmS) 

111 1 
^ +- •♦"T^ + ' ;. ^ + ••• +, + ..* 

1 1.2 1.2.3 |W~1 

the terms which follow the w*^, viz. 

1 1 1 

.— + 1 ; + , r + 8ttc. 

\n [n -f 1 [n + ^ 

are evidently less than the corresponding terms of the geome- 
trical progression 







1 

— + 
n 


1 1 

n ' n 


1 
n 


.^ + &c., 


the 


sum 


of which is 












1 
n ' 


1 

I 
1 -- 


1 


1 




n - 


r«-i' 



n 
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and continually tends to zero as n increases. Therefore the 
proposed series is convergent; and the error in taking the 
aggregate of its first n terms for its sum, is less than the 
quotient of the n^ term divided by n - i. 

108, From the measure of convergency or divergency 
which a geometrical progression furnishes, we shall now pro- 
ceed to deduce one or two other tests as given by Cauchy, 
Caurs d^Ancblyse Algibriqtie. 

109- The series Wj +«2+ ••• + «^»+ ••• is convergent, 
or will become so, if the superior limit of (u^)^ be less than 1, 
when n is infinite ; and divergent in the contrary case. 

1 
Let k denote the superior limit of (u^)^ when n is infinite ; 

and first suppose k <1; also, let a be any magnitude between 

Jc and 1, so that k<a<l; then when n is increased inde- 

finitely, (t^«)" cannot approach indefinitely near to k without 
finally becoming constantly less than a. Therefore it will be 
possible to take for n so large a valu6, that for that and all 
superior values, we may constantly have 

Consequently, the proposed series will finish by always 
having its terms less than the corresponding terms of the 
geometrical progression 

o + a* + •.. +a* + o""*'*+ ... ; 

and as this series is convergent, a being < 1, it follows that 
the proposed series will end by being, i fortiori^ convergent. 

Secondly, suppose ft > 1 ; and take, as before, a between 

1 and At, so that fc > a > 1. Then, when n is indefinitely 

1 
increased, (Uy^ cannot approach indefinitely near to k without 

finally becoming constantly >a; we shall therefore be able 

13 
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1 
to satisfy the condition {u^>a^ or u^>a^y by taking n 

sufficiently large; and consequently we shall always find in 
the series iij + w, + ... + w, + w«+i + ... , an indefinite number 
of terms greater than the corresponding terms of the geo- 
metrical progression a + a* + ... o* + 0*+^ + ..., which is di- 
vergent, a being > l ; and therefore the proposed series will 
end by being divergent. 

110. It may be shewn that if, as n increases indefinitely. 



u. 



u^ remains positive, and the ratio -— continually tends to 



w« 



become equal to a finite quantity k as its limit, the expression 

{u^ continually approaches to the same limit. Hence, the 
test of convergency or divergency in the last Article, may be 
changed into the following, which is more convenient in its 
application ; if the limit of the ratio of u^^^ to u^ when n is 
infinite be less than 1, the series is convergent ; and divergent 
in the contrary case. 

If Ap « 1, this test gives no result in either form. 

111. Suppose that the series tii + t^g + fi3+ ... consists 
of both positive and negative terms; then if Vi, t?29 ^39 &<^* 
be the numerical values of these terms, so that Wj = ±tJi, 
t^2 » ^ t?2' ^^*» ^t ^^ evident that the sum of the proposed 
series can never surpass that of the series f>i + Vg +tJs + ...; 
if therefore, the latter series be convergent, that is, if 



r^).. 



<i» 



the proposed series will be convergent ; or, if the latter series 
finish by having terms greater than any assignable magnitude, 
that is, if 






>i, 
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the same thing will happen to the proposed series, which will 
consequently be divergent. Hence, the above test is ap- 
plicable to series consisting of both positive and negative terms, 
provided we use the numerical vdues of .the terms without 
regard to signs; and it fails, as in the preceding case, when 
Ap= 1. 

112. Let 00 + 0,07 + 0^0^+ ... +a««7* 4- ••• (l.) 

be a series arranged according to positive and ascending 
powers of the variable w; the coefficients being positive or 
negative ; then, by what has been proved, this series will be 
convergent or divergent, according as ka^ (where 

when n is infinite,) is numerically less or greater than 1. 

Hence for all values of w between - - and + 7 the series 

k k 

will be convergent, and for all values of w beyond those 
limits it will be divergent. 

Ex. % + l(^Zi)«. + i(lZiHlZi)«3 + &c. 

1 1.2 1.2.3 



Here -^ « = - 1, when n is infinite; 

a. » + 1 

therefore the series is convergent or divergent, according as c^ 
lies between + 1 and •* 1, or without those limits. 

^ a a* cf ^ 

Ex. _ + _+- + &c.. 



or 
n 
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.'. «= MHM « a, when » is infinite ; 

w. » + 1 

therefore the series is convergent or divergent^ according as 
o < or > 1. 

^ 111 1 

1 [2 |s \n 

s B 0, when « IS infinite, 

t^n n + 1 

therefore the series converges* 

113. The series Wi + 1^ + ... + w« + ... is convergent, 

\ogt f^ 

or will become so, if the inferior limit of * be > 1, when 

n is infinite; and divergent in the contrary case« 

Let k denote the inferior limit of " when n is in- 

finite, and first, suppose k > 1; also let a be any number 
between k and 1, so that k > a> 1. Then when ^ is in- 

, log — 

definitely increased, " or its equal * cannot approach 

logi. ^°«" 

n 

indefinitely near to k without finally becoming constantly 
greater than a. Therefore it will be possible to take for n 
so large a value, that for that and all superior values we 
may constantly have 

log — > a log n > log n% or — > » , or u^<—* 
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Consequently the proposed series wiU finish by always having 
its terms less than the corresponding terms of the converging 
series (Ex. Art. 112). 

and therefore will itself end by being, a fartiorij convergent. 

Similarly, if Ap < 1, it may be shewn that the proposed 
series will finish by being divergent. 

114. The series t^o + t^i -f 1^2 + ^ + ^4 + % + kc. (l), 

each term of which is less than the immediately preceding 
term, 

and the series u^ + ^u^-^ 4iu^ + Su^ + iGu^ ^ &c. (2), 
are convergent or divergent at the same time. 

Suppose series (l) to be convergent and that its sum »ay 
then 

Su^ < 2U4, + 2% + 2U^ -f ^Uj 



<Uq + 2(Wj + % + «8 + &c.) < 2« - «09 

consequently series (2) is convergent. Next, suppose series 
(1) divergent, then 

2«i > «*! + M2, 
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.'. 1*0 + 2wi + 4»8 + «tt7 + 8ic. > «i^ + f#i + fi2 + «s + he 
and is therefore, divergent. 
Ex. Let series (l) be 

then series (2) is 

1 1 1 i . 

a geometric progression convergent when m > I, and di- 
vergent in the contrary case; consequently series (3) will be 
convergent if i» > 1, and divergent if m = or < 1. 

116. The series Wi - Wg + «3 - W4 + &c., is convergent, 
if the numerical value of the terms decreases without limit. 

For, by writing it in the forms 

Ui - (Wj, - Ws) - (W4 - «6) - &c., 

Wi - «jj + (t^ - W4) + &c. ; 

we see that it is >Ui'^u^ and <Ui, and therefore is con- 
vergent. 

Thus, the sum of the series 

1 1 1 1 o 

+ + &c., 

12 5 4 ' 

lies between 1 and ~; also the series 
2 

1 1 1 1 e. 

1« 2« 3« 4" ' 

is convergent for all positive values of m. 
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116. It is also shewn in the work from which these tests 
of convergency are taken, that two converging series all whose 
terms are positive, will by their addition or multiplication 
produce new converging series, whose sums result from the 
addition or multiplication of the sums of the former. 

117* An approximate value of St^, will be obtained by 
taking the aggregate of the converging terms only, in the 
series for ^u, involving J!^u«, Ug^ and the differential coefficients 
of u^; and the error will be less than the last of the con- 
vergent, or the first of the divergent terms. 

We have by Art 6S^ (omitting the index of Ug in the 
second member,) 

S S 

». - C + Ji It - -I » + j-^ d,« - Tj rf> + . . . 



J.U, 



.+(-0"-'^dr-'«+(-i)"2fi.. 0) 

where 21?, - -^^ii- e^*'u- -^!i^ d^^'u + &c. ; 



2« + 2 ' 2» + 4 



or since in general, by Art. 67, 

gg.,1 f 1 1 !_ 1 

[2n " ^ i(2w)»» "^ (4fl-)*' ■*" (eTT)" "^ / ' 

1(2^)*'+* (47r)''+* (6«-)«»+* ^ / ' 



-&c.; 
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series 



as m 



or, adding the terms vertically, and calling the resulting 
»i, »,, &c., 

df+*w df+'w df+'tt 

where t^m " ^ v-^« - ^ ^^^^ + ^ ^p,.^ - &c. 
(2Wfl-)"»+* (2w7r)*'+* (2m7r)^+* 

or «!:«. + (2m.)«t,.»^^_^, 

and integrating this by the met od of parameters. 

Art 72j ^^^ 

*^"" " (2m7rr-»-^ jC(8in2m7rjy.<g"^^u), 

fdnce, (V being an integer, the term multiplied by sin2mir^ 
disappears, and cos2mfl-^ = l; and the arbitrary constant 
28 unnecessary, being already introduced in equation (l); 

^ /•f8in2dr7r sin4a?'7r sinG^ir ^ ] «^^i 

therefore, numerically, 

2i?.< -2 f\ . .- ■>- , .- + ,, .. -t- ... >dj»+^a, 

2» ' 



rhich last quantity lies between 

■"*«-l jB«-1 j -"Sn + l J21I4.1 

-T dr u^ and T ^^—-ct'+'m; 

|2W ' |2w + 2 ' 
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if these be the last of the convergent and the first of the 
divergent terms respectively, of the series for 2w,- Conse- 
quently, the sum of all the diverging terms in the series for 
2^^^, is less than the last of the convergent or the first of the 
divergent terms. 

Interpolation. 

118. When a series of values of a quantity is ob- 
tained either by observation, or laborious calculation, it 
is of great importance to be able to insert other values 
between them, such as would have resulted from a similjir 
observation or calculation, without the labour of performing 
these. This is the object of Interpolation; and in this the 
Calculus of Finite Differences finds one of its chief uses. 
More strictly, Interpolatr^tn of Series is the inserting among 
the terms of a given sei.js, new terms subject to the same law 
as the first. In doing this, the terms of the series are con- 
sidered as particular values of the function which expresses 
its general term, corresponding to a given regular succession 
of indices; and it is the business of Interpolation to dis- 
cover that general term ; or at least to assign such a function 
of the index as shall represent the given series of values, 
and, approximately, all intermediate values. The problem 
thus requiring us to assign the analytical expression of a 
function from a limited number of its numerical values, is 
plainly indeterminate; it is the same as to form the equa- 
tion to a curve which shall pass through a limited number 
of points, whose abscissae represent the values of the inde- 
pendent variable, and the ordinates those of the function, 
without giving the species of the curve; which, as is evident, 
may be done in an infinite variety of ways. But if the 
given terms are numerous, and near to each other, the ex- 
pression for the general term, within the limits of the given 
quantities, may be found to a great degree of accuracy. 

119. There are two principal cases to be considered; 
first, when the given values of /(a?), namely, 

/(^i)> /(^i + h), f(w, + 2 A), ... f{w, + (n - 1)A}, or 

Wi, U29 ^3, ..., «„, 

14 
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as we shall write them, correspond to values of the inde- 
pendent variable, Wu j?i + A, ^, + 2A, &c., a?i + (» - l) A, in 
arithmetical progression. And secondly, when the given 
values /(jr,),/(ar2), •../(ar„), or Wi, «,,... w^ correspond to 
values a?i, djg, ... a?., of the independent variable, not obeying 
any assigned law. 

120. Having given w,, Wj, Ws, &c. «„, « values of a 
function /(^), corresponding to the n values of the indepen- 
dent variable a?!, a?i + A, ... arj + (w - 1) A, to find an ex- 
pression for any intermediate value /(o?! + A). 

Since/(a? + nh) ^f(w) + n A/(af) + ^i^Zil Ay(a?) + &c. 

(Art. 23) changing or into w^, and then replacing nh by A, 

we get 

A A (A — A) 

/(«, + ft) - u, + - A», + ^g^/ A*«. + &c. 

in which, since 

A*W, «= W8-2t«2+tti, 



. , / X (»-l)(»-2) 

A-*W|«Wn-(«-l)w,»i + f^ ^tt..,-&C....JiWi, 

1 .% 



(2). 



if we make A xs o. A, 2 A, &c., (n - l) A, the second member 
assumes the n values fi,, t^, ... t^«; and not only this, but if 
we assume for A any value whatever between and (n - l) A, 
we shall obtain the value of the function corresponding to 
that value of the independent variable. 

121. In applying the above formula, the simplest mode 
is not to calculate At^i, A^t^i, &c. by equations (2), but by 
continued subtraction of the given terms; that is, we must 



Digitized by LjOOQ IC 



107 

write down the series of given values u^, u^, ^•. u„, and sub- 
tract each from the succeeding one; next subtract each of 
these differences from the succeeding difference ; then perform 
the same operation upon the new differences ; and so on, till 
the process terminates ; the first terms of these series of dif- 
ferences are the values of Awi, A*Wi, ..- A*"'«i. Unless 
the terms of the given series by continued subtraction lead to 
a constant difference, the expression for f{af^ + k) will have 
as many terms as the given series has. 

£x. 1. Having given the values of sin 30°, sin^lS sin 32^, 
8inSS<>, to find sin (SO^ + A?'), k being between and 18C/. Here 



«,».5 


A 






t^ = .5150381 


150381 


A* 




«, « .5299193 


148812 


-1569 


A' 


«4 = .5446390 


147197 


-1615 


-46; 



.-. Awi« 0-0150381, A*Wi = - 0*0001569, A'* Wi «- 0-0000046, 
V. sm (30^ + *') = .5 + - At/, + -^^ A^t*! 

A; (A -60) (A? -120) ,, 
4- — i ^--^ ^A ««. 

If A; a 20, it will be found that sin 30^.20' a .5050299, 
which is too large only by a unit in the seventh place of 
decimals. 

Ex. 2. Having given 
log 3.14 « 0.496929, log 3.15 - 0.498310, log 3.16 = 0.499687, 
log 3.17 '^ 0.501059 ; shew that log 3.14159 » 0.497149. 

Ex. 3. Having given t^o t^, t^s, ^49 four right ascensions 
(declinations, longitudes, &c.) of the moon at intervals of 
12 hours, to find its value i hours after the time corresponding 
to the second value. Here A s 12, k ^12 + t, and if the 
required right ascension « tf^ + ^, then 
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V 2 6 / 12 2 Vl2y 6 U2/ 

when developed in powers of — . 

122. Between every two consecutive terms of a given 
series, to interpolate any number of equidistant terms. 

Let Ui^ U29 Wg, &c., Uj^ be the given series, and let w» — 1 
be the number of equidistant terms to be inserted between 
every two consecutive terms ; then the new series will be 

Wij Um+lj tC m+2 9 ... U ^m-l y U29 t£ ^m+l 9 &C. 
fit in fit fit 

If therefore v^^i denote the r + 1*^ term of this series^ 
we have 



-=;±r=/(i + 3 = 



or Vr+i «= w, + — Awi + -^ r- A*«i + &c. 

m 1.2 m* 

Hence, taking r from 1 to m - 1, we get the terms in- 
serted between Ui and Uz ; next taking r from m + l to 
2 w — 1, we get the terms between Uf^ and t^ ; and so on. 
The differences At^i, A*«i, &c., are to be computed by con- 
tinued subtraction as in Art. (121) ; and the series for v^+i 
will have as many terms as the proposed series has, unless 
those terms by continued subtraction lead to a constant dif- 
ference. 

Ex. To insert three equidistant terms between every 
two consecutive ones of the series 1, 7, 15, 28, 49, &c. 
Here wi « 4, and 

A, 6, 8, 13, 21, 

A' 2, 5, 8, 

A^ 3, 3, 
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_ 6r r(r-4) r(r-4)(r-.8) 128 + 192r-4r'-f r* 



, ,, . . 317 504 695 1113 ^ 

and the senes is 1, • — , — , 7, , &c. 

128 128 128 128 



123. The formula of Art. 120 may be presented under a 
different form by changing k into tT — a?i, which gives 

(a? — Xi)(af — tTi— h) •.. {a? — ^i - (n — 2)A} _j 



+ 



w- l.A' 



«-i 



where /(^) is a function of x which, as x assumes the n values 
^15 ^1 + A, &c., ofj + (» - l) A, successively assumes the cor- 
responding values t^j, Wg, ... w„ ; and for any other value of 
w within, or not far beyond, the limits x^ and w^ + (n — l)h, 
it gives the value of the corresponding interpolated term. If 
we put A = 1, the formula is adapted to the case where the 
increment of the principal variable is unity* 



124. In any series of consecutive equidistant values of 
a function, where one is deficient to insert that one. 

Let u^j ^25 ^5 &c., u„ be the values of the function 
corresponding to the values or^, ^^ + A, ... a?^ + (w — 1) A for w. 
Then assuming that A^'^u^^O^ or that the (»-2)*^ dif- 
ferences are constant, which will almost always be the case 
in tabulated results, we have 

(n-l)(»-2) 

A-'^Wi = t^« -• («- l)Wn-i + TIT ^«-2 "■ 8^^- =*= «i = 0, 

1 .2 
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an equation of the first degree from which any one of the 
values as u^ may be found, if the rest be known. Having 
thus completed the system of values, we may interpolate any 
intermediate term f{po^ -i- *) by the method of Art. 120. If 
two values out of n are deficient, then we must suppose 

t^-*u^ a 0, A*"*«8 = 0, or 

In - 2) (w - 3) 
tt«-i - (n - 2) w.., + — w,., - &c. ± 111 - 0, 

1 • z 

(^ - 2) (w - 3) 
«» - (w - 2) «,., + i— W«., - &c. ±«, = 0; 

1.2 

which equations will suffice to determine any two of the values 
in terms of the rest. In the same way any number of 
deficient terms may be inserted. 

Ex. 1. Given the cube roots of 121, 122, 124, 125, to 
find that of 123. 

tt, = 4.946088, «, « 4.959675, W4 = 4.986631 , ttj = 5, 

A*Wi « W5 - 4«f4 + 6««3 - 4^2 4. Uj s ; 

.-. Ws = 5 {4 (w, + W4) - t^, - u^ = 4.973190. 

Ex. 2« Given t^j, t^, f«5, t«6, to find t^, t/^. 
A*Wl = «6 - 4^4 + 6^8 — 4% + t^i B 0. 
A*W2 « «^6 - 4.W5 + 6W4 - 4«, + Wg « ; 
.-. «8 = ^ (- Sttj + IOM2 + 5^8 - 2Wg), 
Wj = ^ (- 2Wi + 5^2 + IOW5 - SWg). 
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125» We next come to the case where the given values 
/(^i)> /(^8)» &c., /K), or Wi, «2,..., «„ 

correspond to values w^^ ^29 •••^n9 T^ot obeying any assigned 
law ; and it is required to determine a rational integral 
function of n — 1 dimensions, fOc)y which shall assume the 
n given values u^, U2y..*u^9 when for w the values Xi, 
^29 ^39 &C.9 .V. are successively substituted. 

Since /(^) is of (n — 1) dimensions, we may assume 

_ /(^) .C, . C2 . ^ C, 

+ ■ ' + ••• + '^ 



H-C2(j7-^i)(ar-a?3).».(a?-a?.)+... + C,(^-a?i)(a?-a?2)...(a?-ar,_i). 

Now make x^w^^ w^y &c., a?., successively; and observ- 
ing that the corresponding values of the first member are 
«i, W8,...w„ we get 

^2 « C2 (^2 "" ^1) (j^i "" '^s) • • • (^2 •" ^a) 



+ «2 7"^- — 'V — ''x y — -^ + &C. 



(*.- 


- a?«) K 


-^s). 


.. (»i 


-^.) 


(^. 


- a^,) (a? . 


-a?,).. 


.(ar. 


■^.) 


(.»,- 


■ ^i) (<"« ■ 


-^.). 


.. («« 


-^.) 


(*. 


-jr,)(a?- 


- a'.) .. 


..(ar- 


a?,-,) 



which is Lagrange^s Theorem for interpolation. 
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Ex. To find a function of w which, when a? « 1, 3, 6, 12, 
shall assume the values 1, 7, 10, - 8. 

•^^ ^'^ 2.5.11 ' 2.3.9 

^ (^ - 1) (<r - 3) (a? - 12) {x - 1) (a? - 3) (a? - 6) 

sTi ' 11.9.3 ' 

126. To determine the maximum or minimum value of 
a function, from three pi its values near its maximum or 
minimum, and the three corresponding values of the indepen- 
dent variable. 

If Ui, %, u^ be the given values of u^ and Xu w^^ ^3 
those of 0?, we have 

_ (a?-a?2) (^-«i?3) (je-^i) (^-^3) (^-^1) (^-^2) 

(^1—^2) (^i-^s) (^2-^1) (^a-^s) (^s""^i) (^3-^2) 

Hence, putting d,w = 0, we find 

«, (a?2 - w^ (2a; - /Pg - o^g) + u^ (w^ - a?i) (2w - ^, - ^Pg) 
+ 2^ (^1 - ^2) (2^ - a?i - /Pg) = ; 
1^1 (a?| - wl) - Wg (d?i - al) + W3 (^! - a?D 



07: 



2Wi (a?g - 073) - 2t*2 (^i - ^3) + 2^3 (a?i — W2) 



the value of j? at which t^ is a maximum or minimum. This 
formula is useful in various Astronomical problems, as for 
instance, to determine the meridian altitude of a heavenly body, 
when an observation exactly on the meridian cannot be ob- 
tained. 
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